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SECTION  I 
INTRODUCTION 


The  advent  of  new  material  and  manufacturing  techniques  has 
resulted  in  the  use  of  thin,  light  weight,  and  integrally  machined 
engine  components  for  aircraft  propulsion.  These  components 
operate  under  high  centrifugal  and  thermal  load  conditions  and 
under  severe  vibratory  environment.  Surge,  stall,  and  flutter  are 
some  of  the  flow  related  vibration  excitation  mechanisms  while  the 
rotor  balance,  misalignment,  and  bearing  induced  instability  are 
mechanical  sources  of  vibrations.  If  turbine  and  compressor  blades 
and  disks  operate  under  vibratory  conditions  for  an  extended  period 
of  time  they  fail  due  to  fatigue  at  critical  locations.  Hence  the 
study  of  flexible  turbine  blades  and  bladed  disk  system  vibrations 
is  an  important  task  in  the  evolution  of  a  design  strategy  to 
ensure  structural  integrity  of  turbine  engine  structural 
components . 

Two  approaches,  analytical  and  experimental,  need  constant 
upgrading  to  meet  the  design  and  test  requirements  of  the  critical 
turbine  components.  The  geometric  intricacy,  material  and  geometry 
properties  uncertainty,  and  the  nonlinear  response  of  these 
components  limit  the  deterministic  analysis  application.  The 
upgraded  analysis  will  have  to  consider  these  effects.  Also, 
strain  gage-slip  ring  or  strain  gage-telemetry  measurement  methods 
have  pro^'ided  the  stress  and  vibration  data  at  selected  points  on 
turbine  and  compressor  blades  in  operation.  Laser  interferometry 
has  been  used  for  the  full-field  vibration  data  in  laboratory 
conditions.  However,  these  measurement  methods  suffer  from 
disadvantages  such  as  limited  data,  short  time  operation,  long 
downtime  for  sensor  repair  and  installation,  and  limited 
accessibility. 

Two  aspects  are  readily  discernable  from  this  discussion. 
First,  we  need  to  develop  reliable  and  large  vibration  data 
acquisition  systems  for  the  continuous  monitoring  of  several  blades 
of  the  many  critical  stages  of  turbine  and  compressor  rotors. 
Second,  the  new  analytical  methods  will  have  to  consider  the 
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statistical  variations  in  turbine  engine  material  and  geometric 
properties. 

1.1  PROGRAM  THRUSTS 

Recognizing  the  present  gaps  in  the  existing  turbine  engine 
design  technology,  this  research  program  identified  and  set  forth 
multifold  thrusts  in  developing  new  measurement  devices  and 
upgrading  the  analysis  techniques.  The  analytical  and  experimental 
methods  developed  were  fundamental,  laboratory  oriented,  and 
possess  potential  for  application  on  prototype  engine  components. 

The  research  and  development  thrusts  from  this  program  were: 

1.  Develop  new  noncontacting  vibration  measurement  systems 
for  rotating  bladed  disk  assemblies, 

2.  Develop  finite  element  based  deterministic  and 
probabilistic  analytical  methods  to  determine  the  dynamic 
response  of  turbine  blades  and  bladed  disk  systems, 

3.  Apply  and  evaluate  of  the  analysis  and  measurement  systems 
to  define  the  vibratory  response  of  rotating  turbine 
disks. 

1 . 2  PROGRAM  SUMMARY 

During  the  period  1  September  1985  to  31  May  1989,  the  work  on 
the  Structural  Testing  and  Analysis  Research  of  Turbine  Engine 
Components  was  done  at  the  Innovative  Structural  Research 
Laboratory  (ISRL)  of  the  Aero  Propulsion  and  Power  Laboratory, 
Wright-Patterson  Air  Force  Base,  and  also  at  the  University  of 
Dayton  Research  Institute  (UDRI) . 

A  four-channel  fiber  optic  plane-of-light  digital  data 
acquisition  and  data  display  system  was  designed,  fabricated,  and 
demonstrated  for  tip  deflection  measurements  on  rotating  blades. 
This  system  is  laboratory  suited,  equipped  to  compensate  for  shaft 
speed  changes  during  one  revolution,  and  has  a  fast  turnaround 
capability. 


2 


A  novel  two-channel  digital  Z-grid  (grid  wire)  instrumentation 
system  was  designed,  fabricated,  and  tested  to  measure  an 
instrumented  blade  tip  deflection  of  a  rotating  bladed  disk.  Both 
the  plane-of-light  and  Z-grid  wire  measurement  systems  were 
integrated  such  that  the  data  acquisition  was  achieved  through  a 
common  device. 

Suitable  data  analysis  software  programs  were  developed  to 
generate  time  and  frequency  spectra  from  the  test  data  and  to 
extract  the  blade  tip  deflection  of  all  blades  from  the  optical 
system  data  and  one  blade  from  the  grid  wire  system  data.  A 
quantitative  fringe  interpretation  method  for  the  differential 
double  pulse  interferometry  was  also  developed  and  applied  to 
define  full-field  vibration  patterns  of  a  flexible  rotating  bladed 
disk. 

A  joint  effort  with  Arnold  Engineering  Development  Center 
(AEDC)  and  United  Technologies  Research  Center  (UTRC)  to  excite  a 
rotating  bladed  disk  in  a  backward  traveling  wave  and  compare  the 
performance  of  the  AEDC/UTRC  and  ISRL/UDRI  four-channel  optical 
systems  concurrently  was  successfully  completed. 

A  broad  range  of  finite  element  based  deterministic  and 
probabilistic  analysis  methods  was  developed  and  applied  to  the 
normal  mode,  frequency  response,  and  sensitivity  calculations  of 
simple-to-intricate  turbine  components.  The  UDRI  developed  PROTEC 
computer  code  conducts  sensitivity  and  probabilistic  calculations 
for  the  turbine  structures. 

Results  from  this  research  effort  have  been  documented  in 
seven  interim  reports  and  presented  in  technical  papers  at  several 
conferences  during  the  period  of  this  program. 

1.3  ORGANIZATION  OF  THE  FINAL  REPORT 

The  work  performed  to  achieve  the  goals  of  this  research 
program  was  described  in  detail  in  22  bimonthly  and  7  interim 
technical  reports.  However,  only  brief  descriptions  of  important 
contributions  to  this  project  are  included  in  the  final  report. 
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The  following  references  and  reports  give  detailed  information  for 
specific  projects; 

1.  "Turbine  Blade  Data  Acquisition  System  -  Technical 
Reference  Manual,"  by  Michael  J.  Gutman  (1). 

2.  "Turbine  Blade  Data  Acquisition  System  --  User's  Guide,"  by 
Michael  J .  Gutman  ( 2 ) . 

3 .  "Turbine  Blade  Data  Acquisition  System  -  Software 
Reference,"  by  Michael  J.  Gutman  et  al.  (3). 

4.  "Z-Grid  Signal  Conditioning  System  -  User's  Guide  and 
Technical  Reference,"  by  J.  Michael  Aulds  et  al.  (4). 

5.  "Probabilistic  Finite  Element  Analysis  of  Dynamic 
Structural  Response,"  by  R.  A.  Brockman  et  al.  (5). 

6.  "Evaluation  of  Plane-of-Light  Noninterference  Stress 
Measurement  Systems  for  Measuring  Bladed  Disk  Vibrations," 
by  M.  Swaminadham  and  Robert  J.  Dominic  (6) . 

7.  "Data  Analysis  Guide  for  the  Noncontacting  Blade  Deflection 
Measurement  System,"  by  Thomas  W.  Held  (7). 

An  outline  of  the  contents  of  this  report  follows. 

SECTION  I  outlines  the  need  for  noninterference  stress  and 
deflection  measurement  methods  development  for  rotating  turbine 
components.  It  also  provides  documentation  format  for  the  final 
report . 

SECTION  II  describes  the  finite  element  model  generation  and 
modal  analysis  methods  for  flat,  twisted,  and  circular  plates,  and 
bladed  disk  assemblies.  Cases  of  numerical  investigation  on  the 
tuned  and  mistuned  bladed  disk  frequency  response  are  detailed.  A 
procedure  and  its  application  to  suppress  singularities  in  the 
analysis  of  advanced  compressor  blades  and  bladed  disk  assemblies 
to  save  cost  and  computational  effort  has  been  described. 

SECTION  III  deals  with  the  stochastic  variations  in  turbine 
structural  element  material  and  geometric  properties  and  develops  a 
rationale  for  deriving  the  frequency  and  mode  shape  sensitivities. 
Geometric  model  generation,  sensitivity  and  probabilistic 


4 


calculation,  and  output  presentation  are  coded  in  the  PROTEC 
computer  program.  This  section  illustrates  the  application  of 
PROTEC  for  the  deterministic  and  probabilistic  analyses  of  simple- 
to-complex  structural  components. 

SECTIONS  IV  and  V  deal  with  the  important  aspects  of  the 
design,  development,  and  test  tasks  of  the  plane-of-light  and  Z- 
grid  instrumentation  systems.  The  development  of  four  optical 
sensors  and  their  integration  with  a  10-MHz  clock  to  digitize  and 
store  the  optical  probe  signals  are  described.  The  principles 

% 

involved  in  the  digital  design  and  grid  wire  layout  of  the  Z-grid 
system  for  measuring  tangential  and  axial  displacement  components 
of  one  single  blade  of  a  rotating  disk  are  also  described  here. 
Features  common  to  both  optical  and  grid  wire  systems  are 
mentioned,  and  the  steps  to  acquire  and  display  blade  data  via  a 
common  front  end  and  control  monitor  are  explained. 

SECTION  VI  discusses  speed  compensation  and  data  analysis 
software  procedures  for  the  stationary  and  traveling  wave  analysis 
of  vibrating  bladed  disks.  The  steps  to  resolve  axial  and 
tangential  components  of  a  blade  tip  deflection  are  described. 
Fringe  interpretation  for  double  pulse  interferograms  and  the 
derivation  of  scaling  factors  for  the  three  disk  mode  shapes  are 
described  in  this  section. 

SECTION  VII  deals  with  the  hardware  design  and  disk 
instrumentation  to  excite  and  measure  response  from  two  disks  in 
stationary  and  traveling  wave  resonances.  Excitation  magnets  and 
sensor  probe  location  details  are  illustrated.  Details  of  the 
three  test  series  conducted  with  UTRC  are  provided.  Problems  of 
triple  pulse  interferometry  and  nonlinear  disk  response  are 
•  addressed.  Time  average  mode  shape  interferograms  for  static 

bladed  disks  are  furnished. 

SECTION  VIII  summarizes  the  analytical  and  experimental 
frequencies,  critical  speeds,  and  mode  shapes.  Discrepancy  in 
analytical  predictions  and  measurements  is  discussed.  Resonance 
speed,  mode  shape,  and  phase  changes  due  to  a  stationary  harmonic 
and  static  force  fields  are  explained.  The  wave  buildup  and  the 
amplitude  collapse  in  the  linear  and  nonlinear  disk  response 
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regimes  are  provided.  Tables  in  this  section  compare  the  UDRI 
analysis  predictions  with  the  optical  and  grid  wire  system 
measurements.  In  addition,  AEDC  and  UDRI  NSMS  evaluation  tests  and 
results  are  also  provided. 

Section  IX  presents  the  conclusions  based  on  the  analytical, 
experimental,  and  instrumentation  development  effort  on  vibrations 
of  rotating  turbine  engine  components.  Recommendations  to  the  Air 
Force  for  the  future  turbine  dynamics  research  is  also  provided. 
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SECTION  II 

VIBRATION  ANALYSIS  BY  DETERMINISTIC  APPROACH 


A  broad  analytical  capability  for  performing  deterministic 
dynamic  analysis  of  turbine  engine  components  has  been  developed. 
Components  ranging  from  flat  plate  to  complex  integral  bladed  disk 
assemblies  were  analytically  modeled  and  analyzed  for  normal  mode 
and  forced  response  calculations.  Wide  use  of  the  general  purpose 
finite  element  codes  to  model  turbine,  compressor  structural 
components  was  accomplished.  MAGNA  (Material  and  Geometrical  Non- 
Linear  Analysis  (8)),  SAP  (Structural  Analysis  Program  (9)),  and 
COSMIC  NASTRAN  (10)  computer  programs  have  been  implemented  on  CDC 
Cyber,  VAX  and  CRAY/XMP  computer  systems  and  the  pre-  and  post¬ 
processing  capabilities  of  these  codes  are  interchangeably  used  for 
analytical  modeling  and  results-plotting  purposes. 

Application  of  these  codes  for  the  dynamic  analysis  of  turbine 
blades  and  bladed  disks  is  described  below. 

2.1  FINITE  ELEMENT  MODEL 


Based  on  a  linear  elastic  analysis,  the  finite  element 
approximation  begins  from  the  constitutive  relationship: 


a 


ij 


^ijkl  ^kl 


(1) 


Because  of  the  symmetry  of  the  stress  tensor,  and  of  the  strain 

tensor,  the  property  of  a being  a  potential  function  of 

and  of  the  isotropy  of  the  material,  equation  (1)  reduces  to 


a 


ij 


E 

1+1/ 


1-2V 


^kl) 


(2) 


For  small  displacements,  the  linear  strain  tensor  is  given  by 


c 


ij 


1 

O 


(3) 


where  Uj^  defines  the  displacement  vector.  The  principle  of  virtual 
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work  requires  for  rhe  stresses  and  inertial  forces  the  relationship 


I  ‘'ij  “ij  =  -e  I  dV  (4) 

V  V 

which  yields  by  completion  with  the  virtual  work  done  by  the  nodal 
force  of  a  finite  element,  the  equations  of  motion 

+  k^  =  s^  (5) 

where  s^  is  the  nodal  force  vector,  is  the  modal  displacement 
vector  and  and  k^  are  the  element  mass  and  stiffness  matrice, 
respectively . 

For  a  complex  structure  comprising  j  =  1,  2,...n  finite 
elements  the  description  follows  by  adding  the  n  equations  (5)  in 
appropriate  manner  to  give 

M  0  +  k  0  =  S  (6) 

and  which  yields  by  introducing  the  boundary  conditions  and  the 

A 

assumption  of  harmonic  vibrations,  0=0  cos  «t 

p  *  *  ^ 

(-  M  +  k  )  0  =  O  (7) 

The  solution  of  the  eigenvalue  problem  gives  the  natural  frequency 
(radians)  and  the  corresponding  eigenvectors  0  which  determine  the 
mode  shapes  of  vibrating  structural  components. 

2.2  VIBRATION  ANALYSIS  OF  TWISTED  BLADES  AND  CIRCULAR  DISKS  BY 
FINITE  ELEMENT  METHODS 

A  general  purpose  in-house  finite  element  computer  program 
(MAGNA)  has  been  used  to  analyze  the  dynamic  characteristics  of  the 
basic  components  of  bladed  disk  structures.  Cantilevered  flat  and 
twisted  plates  and  centrally  clamped  circular  plates  having  the 
same  geometry  and  material  properties  as  those  tested  in  the  ISRL 
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were  analyzed  for  eigenvalues  and  eigenfunctions.  The  blade  was 
assumed  to  have  a  linear  twist  angle  about  its  centroidal  axis. 
Pre-  and  post-processing  capabilities  of  the  MAGNA  program  have 
been  used  to  evaluate  several  natural  frequencies  and  mode  shapes 
of  the  test  objects.  Results  from  the  flat  and  pretwisted 
cantilevered  plate  and  centrally  clamped  circular  plate  analyses 
are  presented  in  Tables  1  and  2  respectively.  Mode  shapes  of 
twisted  plates  using  MAGNA' s  contour  plotting  features  depicting 
constant  amplitude  vibrations  are  presented  in  Figure  1.  Circular 
plates  dynamic  deformed  shapes  as  shown  in  Figure  2  were  generated 
by  the  MAGNA's  geometry  plot  (G-plot)  capability.  These  display 
features  enhanced  our  understanding  of  the  blade-alone  and  disk- 
alone  eigenfunctions. 


2.3  BLADED  DISK  VIBRATIONS 

Finite  element  analysis  method  has  been  applied  to  integrate 
circular  disk  and  blades  into  one  assembly  and  to  evaluate  the 
natural  frequencies  and  mode  shapes  of  a  12-in. -diameter ,  12- 
bladed  disk  assembly.  SAP  finite  element  code  has  been 
implemented  on  the  CDC  computer,  and  several  annular  disks  were 
modeled  using  thin  shell  plate  elements.  Several  frequencies 
obtained  through  SAP  agreed  well  with  MAGNA' s  results  and 
experimental  values.  However,  the  SAP  IV  program  provided  output 
in  written  data  only  and  did  not  have  plotting  capability.  To 
interpret  the  SAP  results  qualitatively,  a  computer  program  has 
been  written  to  convert  SAP-output  into  a  format  accepted  by  the 
plotting  program  used  by  the  MAGNA  code.  This  feature  allowed  us 
to  plot  model  geometry  and  vibration  contours  of  the  bladed  disk. 

A  comparison  of  analytical  and  experimental  frequencies  of  a 
bladed  disk  is  shown  in  Table  3.  Finite  element  model  and  the 
natural  modes  of  a  12-in. -diameter  bladed  disk  analyzed  by  SAP  IV 
and  plotted  using  MAGNA' s  plotting  capabilities  are  shown  in 
Figures  3a  and  3b. 
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TABLE  1 

NATURAL  FREQUENCIES  (HZ)  OF  FLAT  AND  PRETWISTED  PLATES 


Geometry  and  Material  Properties 

Length  x  width  x  thickness:  2  x  2  x  5/64  in. 


Maximum  twist  angle  at  the  tip:  45“ 


Young's 

Poisson' 

Density 

Modulus : 

s  ratio  : 

• 

• 

30  X  10^  psi 

0.3 

0.000735  Ib-m/in 

Mode 

Flat  Plate 

Twisted  Plate 

IB 

735 

774 

IT 

1641 

3446 

2B 

4351 

4630 

3B 

5784 

6725 

2T 

6129 

8342 

4B 

9894 

10338 

B  - 

■  Bending,  T  -  Torsion 

TABLE 

2 

NATURAL  FREQUENCIES  (HZ) 

OF  CIRCULAR  PLATES 

Geometry  and  material  orooerties 

Outer  radius: 

6  in.  and  4  in. 

Inner  radius: 

0.5  in. 

Thickness  : 

5/64  in. 

Young's  Modulus: 

30  X  10^  psi 

Poisson's  ratio 

:  0.3 

Density  : 

0.000735  Ib-m/in 

MAGNA 

LEISSA  (11)  BLEVINS  (12) 

Mode  (n,  s) 

4”  Dia. 

4"  Dia. 

1  (1,  0) 

199.6 

168.5 

168.1 

2  (0,  0) 

238.4 

215.5 

215.4 

3  (2,  0) 

281.1 

281.0 

280.8 

4  (3,  0) 

- 

594.3 

594.5 

5  (4,  0) 

•• 

1305.2 

1304.8 

n 

-  nodal  dia,  s  - 

nodal  circle 

10 
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TABLE  3 

NATURAL  FREQUENCIES  (HZ)  OF  A  12-IN. -DIAMETER,  12-BLADED  DISK 


Mode 

(n,s) 

Analytical  (SAP) 

Experimental 

(Holography) 

(1,0) 

73.3/74.1 

68.4/68.8 

(0,0) 

91.9 

84.5 

(2,0) 

106.3/107.7 

104.9/108.0 

(3,0) 

199.8/203.8 

193.1/193.5 

(4,0) 

292.1/299.0 

269.6/271.1 

(0,1) 

367 

n  -  nodal  diameter,  s  -  nodal  circle 


2.3.1  Centrifugal  Force  Effect  on  Disk  Modal  Characteristics 

Rotational  speeds  have  considerable  effect  in 
increasing  resonance  frequencies  of  rotating  blades  and  bladed 
disks  of  turbine  engines.  Solutions  of  the  equations  of  motion  of 
rotating  bladed  disks  enabled  us  to  study  the  centrifugal  force 
effect  on  natural  frequencies  and  mode  shapes.  Some  details  of  the 
analytical  basis  are  described  here. 


are 


The  equations  of  motion  of  a  rotating  blade  disk  (13) 


[M]  (u)  +  [B]  (0)  +  [[K]  -  (u)  =  (P) 


(8) 


where 


(u) 


U 

y 


(9) 
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Inner  Die  =  1/2  inch 
Outer  Dia  =  8  indi 
Blade  Length  ^  2  indi 
Overall  Disk  Dia  »  12  inch 

Thickness  =  0.078  inch 


gure  3a.  Finite  Element  Model  of  a  12 -Inch -Diameter  Disk. 


3  D  Sin  0,  199.8  Hz  3D  Cos  0,  203.8  Hz 

Figure  3b(cont'd).  Mode  Shapes  of  12-Inch-Di2uneter  Disk  (SAP) 


[M]  =  =  [Tgg]’’ 

-  m 

o 

a  0 

°] 

f^GB^ 

(10) 

o 

m 

[B]  = 

1 - 1 

0  0  o 

o 

o 

m 

-2] 

o 

f^GBl 

(11) 

[K]  =  [K] 


Global 


and  CM^J  =  =  CTsbI’' 


0  0  0-1 
o  m  o 
o  o  m 


(12) 

(13) 


A  special  feature  of  [M^]  matrix,  associated  with  0  is 
that  it  contributes  to  stiffness  in  the  plane  orthogonal  to  the 
rotational  axis.  Elements  of  the  stiffness  matrix  (Equation  12), 
were  generated  from  four  node  quadrilateral  plate  elements  with 
five  d.o.f.  per  node  (three  in  plane  translations  and  two 
rotations) .  Appropriate  boundary  conditions  were, 

at  the  inner  radius,  {u}={u^}=0  (14) 

and  at  the  free  ends,  (u^^)  =  *  0.  (15) 

When  nonstructural  concentrated  mass  (M)  was  used  at  a 
grid  point  to  mistune  a  bladed  disk,  a  contribution  to  the  kinetic 
energy  was  added.  The  expression  for  T  was  appropriately  changed. 

The  form  of  mass  matrix  for  concentrated  mass  (M)  is 


[M]  = 


M 


M 


Sym 


M 


ZM  -YM 


-ZM 


XM 


YM  “XM 


^xx  ”^xy  ^xz 
^yy  "^yz 


zz 


(16) 
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2.3.2 


The  equation  of  notion  of  a  bladed  disk  systen  with  no 

danping  is 


[M]  {u}  +  CK]  {u)  =  {P}  (17) 

Here,  (u)  are  the  displacements  and  (P)  are  the  loads  acting  on  the 
system.  For  normal  mode  analysis,  the  above  equation  gets 
uncoupled,  and  each  mode  is  represented  by  one  equation  of  motion. 


Let  u.  =  • 

1  j=l  13  D 

(18) 

or  (u)  =  [0]  ii) 

(19) 

Equation  (17)  now  becomes 

[M]  [0]  (h  +  [K]  [0]  (f)  = 

{P) 

(20) 

t^]^[M]  [0]  (C)  +  [0]’’(K)  [0]  {()  = 

[^]j{P) 

(21) 

and  letting  =  1  (0^}'^[K]  (0.) 

(22) 

and  /ij  =  {0j) 

(23) 

Equation  (21)  becomes 

{fj]  +  (fj)  * 

(24) 

where  j  denotes  the  mode  number 

i  denotes  the  displacement  coordinate 

To  obtain  displacement  in  frequency  domain,  Laplace 
transformation  with  S  =  S/St  is  used  to  get 


(u)  =  [0]  - \ [0]’’{P.)  (25) 

M.  (s'"  +  w‘)  J 
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Analogous  to  experimental  technique,  wherein  a 
sinusoidal  forcing  function  sequentially  sweeps  through  all  the 
modes  of  a  bladed  disk,  a  numerical  simulation  of  a  sinusoidal 
forcing  function  of  the  following  form  was  implemented. 


P(u)  =  A 

P(«)  =  Aue^*^ 

P(w)  =  A[cosu+isinti)]  (26) 

2.3.3  Effects  of  Sneed  and  Concentrated  Masses-A  Numerical 

Study 

A  12 -in. -diameter,  12-bladed  disk  was  considered  for 
analysis.  A  reasonable  number  of  4  node  quadrilateral  plate 
elements  with  five  degrees  of  freedom  per  node  was  used  in  the 
analysis  solution  format  (COSMIC  NASTRAN  Rigid  Format  13)  to 
evaluate  twin  modal  frequencies  of  the  bladed  disk. 


Natural  frequencies  and  mode  shapes  of  a  bladed  disk 
with  the  following  dimensions  were  determined  for  several  speeds. 


Inner  diameter  of  the  disk 
Outer  diameter  of  the  disk 
Blade  length 

Overall  outer  diameter  of  bladed  disk 
Blade  width 

Angle  of  twist  at  the  tip  of  blade 
Number  of  blades 
Thickness  (constant) 

Young's  Modulus 
Poisson's  ratio 
Density 


1  in. 

8  in. 

2  in. 

12  in. 

2  in. 
.o 


=  45 
=  12 
=  5/64 


in. 


=  30x10  psi 
=  0.3 

=  0.000735  Ib-m/in' 


Concentrated  mass,  M^  =  3.90  grams,  M2  =  7.8  grams. 

From  Table  4,  it  is  obvious  that  the  centrifugal  force 
effect  on  a  tuned  system  increased  its  natural  frequencies. 
Repeated  (double  and  equal)  frequencies  were  present  for  every 
diametral  (circumferential)  mode  except  for  zero  circular  (radial) 
mode.  Mode  shape  corresponding  to  the  second  of  each  pair  of 
frequencies  was  rotated  by  an  angle  governed  by  ir/2r,  where  n  is 
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the  number  of  diametral  mode.  For  example,  in  ID  mode,  nodal 
diameters  were  orthogonal  to  each  other,  while  In  2D  mode,  nodal 
diameters  are  rotated  by  45  degrees  and  so  on.  Contour  lines 
corresponding  to  normal  displacements  were  symmetrical  with  respect 
to  x-y  axes.  Mode  shapes  for  a  tuned  system  are  presented  In 
Figure  4. 


TABLE  4 

ROTATIONAL  SPEEDS  VS.  NATURAL  FREQUENCIES  (HZ) 
OF  A  TUNED  BLADED  DISK 


Speed 

from) 

ID 

OC 

Modes 

2D 

3D 

4D 

0 

76.5 

92.1 

110.2 

200.3 

276.9 

76.5 

110.2 

200.3 

276.9 

1000 

78.6 

93.5 

112.5 

202.0 

278.3 

78.6 

112.9 

203.6 

281.3 

1700 

82.2 

95.8 

116.7 

205.1 

280.8 

82.2 

117.2 

206.8 

283.1 

2000 

84.3 

97.2 

119.2 

206.9 

282.3 

84.3 

119.6 

208.6 

285.3 

4000 

104.0 

110.4 

142.8 

225.8 

298.0 

104.0 

143.3 

227.3 

300.0 

5000 

116.6 

119.1 

^58.2 

238.8 

309.2 

116.6 

^8.6 

240.3 

312.0 

6000 

130.1 

128.7 

175.3 

253.9 

322.3 

_ 130.1 

_ 175.6 

4 _ 

_ 325.0 _ 

Mlstunlng  due  to  concentrated  masses  at  specific 
locations  on  the  disk  was  studied.  Changes  in  frequencies  are 
presented  in  Table  5.  With  mlstunlng  of  this  nature,  repeated 
roots  will  degenerate  Into  two  modes.  The  dependency  and  extent  of 
this  degeneration  depends  on  the  magnitude  and  the  relative 
positions  of  mlstunlng  masses.  When  only  one  mass  at  the  outermost 
radius  was  used  to  deliberately  mlstune  the  disk,  one  of  the  roots 
lowered  In  frequency  while  the  other  remained  essentially  the  same. 
This  effect  was  seen  In  all  diametral  modes.  Rotational  speeds 
stiffen  the  disk  and  Increase  the  numerical  value  of  the  resonant 
frequencies. 
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Fiqure  4.  12-lnch-Diameter  Tuned  Disk  Mode  Shapes. 


TABLE  5 

MISTUNED  DISK  SPEED  VS.  NATURAL  FREQUENCIES  (HZ) 

CONCENTRATED  MASS  M.  OF  7.8  GRAMS 
AT  GRID  POINT  257 


Speed  Modes 


3D 

4D 

0 

54.2 

86.6 

100.0 

163.6 

236.2 

76.6 

110.5 

200.3 

279.3 

1000 

57.3 

88.2 

102.2 

164.4 

78.6 

112.9 

202.1 

1700 

62.7 

91.1 

106.3 

169.4 

82.2 

117.4 

205.5 

2000 

65.6 

92.7 

108.7 

172.2 

244.6 

84.3 

119.9 

207.4 

284.9 

3000 

77.1 

99.4 

119.1 

184.1 

254.5 

93.0 

131.0 

216.0 

291.9 

5000 

104.4 

117.0 

148.1 

215.3 

282.3 

116.4 

160.8 

241.8 

312.9 

6000 

118.6 

130.2 

165.3 

233.2 

299.0 

256.9 

326.5 

with  a 

Single  mass 

of  7.8  grams 

,  equivalent 

of  one 

blade  mass 

(same  trend  was  evidenced  for  half  the  mass) 

,  and  for 

the  case  of 

ID  mode 

,  mode  that 

corresponds 

to  lower  frequency  had 

distorted  shape.  The  mode  corresponding  to  the  unaltered  frequency 
remained  symmetrical  with  the  nodal  line  passing  through  the 
mistuning  mass.  The  nodal  diameter  was  attached  to  concentrated 
mass,  and  the  mass  does  not  participate  in  mode  formation.  This 
trend  is  true  for  all  diametral  modes  with  a  single  concentrated 
mass.  Circular  modes  were  also  disfigured  in  the  sense  that  the 
disk  half  containing  the  concentrated  mass  had  less  displacement 
compared  to  the  other  half.  Also,  the  center  of  the  circular 
contour  lines  moved  away  from  the  geometric  center  of  the  disk.  In 
2D,  3D,  etc.,  modes,  one  of  the  split  modes  was  altered  while  the 
other  remained  unaltered.  In  the  unaltered  modal  configuration, 
nodal  diameters  pass  through  the  mistuning  concentrated  mass.  The 
results  for  a  mistuned  mass  of  3.9  grams  are  presented  in  Table  6. 


with  the  increase  of  speed,  mode  shapes  corresponding 
to  the  lower  and  unaltered  frequencies  behaved  as  before  except 
that  both  frequencies  raised  in  proportion  to  the  speed  of  rotation 
and  the  displacements  were  less.  Mode  shapes  for  a  case  of  » 

7.8  grams  and  N  =  2,000  rpm  are  presented  in  Figure  5.  The 
centrifugal  effect  of  the  concentrated  mass  reduces  the  resonant 
displacements  in  one  of  the  split  modes. 


A  study  with  two  equal  concentrated  masses  at  180°  and 
also  at  90°  to  each  other  was  conducted,  and  the  effect  of  these 
two  masses  was  to  lower  both  frequencies  of  the  diametral  and 
circular  modes  as  shown  in  Table  7.  The  mode  shapes  remained 
symmetrical  in  the  first  case,  but  in  the  latter  case,  nodal  lines 
did  not  attach  themselves  to  concentrated  masses. 


TABLE  6 

MISTUNED  BLADED  DISK  SPEED 
VS.  NATURAL  FREQUENCIES  (HZ) 

CONCENTRATED  MASS  M  OF  3.9  GRAMS 
AT  GRID  POINT  257 


Speed 

(rpm) 

ID 

OC 

Modes 

2D 

3D 

4D 

0 

64.5 

88.1 

101.8 

170.7 

241.1 

76.6 

110.5 

200.3 

279.3 

1000 

67.0 

89.6 

104.1 

173.0 

243.0 

78.6 

112.9 

202.1 

280.7 

1700 

71.4 

92.3 

108.2 

177.1 

246.5 

82.2 

117.3 

205.3 

283.3 

2000 

76.6 

95.5 

113.4 

182.1 

250.8 

86.7 

122.6 

209.4 

286.6 

3000 

84.0 

100.2 

121.0 

189.6 

257.3 

93.0 

130.4 

215.5 

291.6 

5000 

109.3 

117.7 

150.0 

218.0 

282.8 

116.6 

159.7 

240.0 

312.2 

6000 

122.8 

130.3 

167.3 

235.1 

298.7 

_ 128.6 

177.0 

255.4 

_ 325.5 
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vibration-amplitudes  for  the  grid  points  257,  188, 
and  170  were  evaluated  as  a  function  of  excitation  frequency  in  a 
simulated  sine  dwell  test.  Variations  in  amplitudes  for  small 
increments  of  frequencies  in  close  proximity  to  resonances,  for 
the  ID,  OC,  2D  and  3D  modes,  are  presented  in  Figure  6.  Modal 
damping  values  were  estimated  from  these  graphs,  by  the  halfpower 
point  bandwidth  method.  Resonances  were  widely  separated  even  for 
split  modes,  and  no  response  contribution  from  other  modes  was 
evident.  Phase  differences  in  responses  for  the  grid  points  257 
and  170  which  were  spatially  offset  by  90°  are  seen  in  diametral 
modes.  Modal  damping  values  are  quite  small  and  correspond  to 
hysteresis  damping  of  the  disk  material.  The  results  from  a  study 
on  bladed  disk  mistuning,  rotational  effects  and  frequency  response 
can  be  found  in  reference  (14). 

TABLE  7 

MISTUNED  BLADED  DISK  DYNAMICS 

CONCENTRATED  MASS  M-  AND  M, ,  EACH  OF  3 . 9  GRAMS 
AT  GRID  POINTS  257  AND  230  AT  90* 


Speed 

(rpm) 

ID 

OC 

Modes 

2D 

3D 

4D 

0 

64.1 

86.8 

92.6 

159.2 

229.4 

64.8 

110.5 

177.9 

266.5 

1000 

66.7 

88.2 

95.1 

161.6 

231.3 

67.2 

113.0 

180.1 

268.6 

1700 

71.3 

90.9 

99.7 

166.0 

231.8 

71.4 

117.4 

184.0 

272.4 

2000 

76.5 

94.1 

105.2 

171.5 

239.1 

76.7 

122.9 

188.9 

277.1 

3000 

83.5 

98.8 

113.2 

179.5 

245.7 

84.3 

130.9 

196.2 

284.2 

5000 

108.0 

116.9 

143.1 

209.7 

271.6 

110.5 

160.8 

224.1 

311.3 

6000 

125.0 

120.1 

160.6 

227.6 

287.8 

_ 128.6 

_ 178.5 

_ 241.0 

_ UJLxl _ 
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2.4  ADVANCED  COMPRESSOR  BLADE  VIBRATIONS 

Two  finite-element  programs  were  used  to  study  the  resonant 
characteristics  of  an  advanced  compressor  blade.  SAP  IV  with 
quadrilateral  elements  allowed  tuning  of  the  blade  model  to 
represent  its  half-plate  and  half-shell  geometric  appearance.  The 
COSMIC  NASTRAN  program  with  triangular  plate  elements  was  also  used 
to  analyze  blade  vibrations.  An  added  feature  of  NASTRAN  was  its 
accommodation  of  rotational  effects.  SAP  analysis  provided  a  close 
agreement  with  the  experimental  determination  of  several  lower 
order  modes. 

The  titanium  blades  of  an  axial  flow  compressor  disk  had  a 
transonic  airfoil  cross  section,  were  integrally  machined  with  the 
disk,  and  had  an  aspect  ratio  of  1.255.  These  blades  had  a 
meridional  chord  curved  at  the  base  which  gradually  leaned  out 
at  the  tip.  Leading  edge  sweep  of  the  blade  tip  was  obtained  by 
increasing  the  chord  length  in  the  midspan  region.  Nonradial 
stacking  of  the  centroids  of  the  cross  sections  towards  the  tip  and 
reflex  curvature  in  the  meridional  plane  of  the  leading  edge  were 
used  to  limit  excessive  stresses  on  one  side  of  the  airfoil 
surface.  All  these  features  of  the  blade  profile  were  obtained 
from  the  geometry  file  supplied  by  the  PW  and  also  from  the 
reference  (15) .  Variations  in  chord  length,  thickness,  stagger 
angle  and  the  stacking  centroids  can  easily  be  observed  in  Figure 

7. 

2.4.1  SAP  Analysis 

The  SAP  IV  program  was  employed  with  plate  and  shell 
elements  to  analyze  the  blade  for  natural  frequencies  and  mode 
shapes.  The  model  for  the  UDRI  SAP  IV  analysis  is  shown  in  Figure 

8.  The  SAP  IV  analytical  model  had  143  nodal  points  and  120 
quadrilateral  plate/shell  elements,  with  grid  points  131  through 
143  completely  restrained  to  simulate  the  fixed  root  condition  of 
the  blade.  Each  nodal  (grid)  point  had  six  degrees  of  freedom  and 
the  analysis  method  can  selectively  free  or  restrain  the  rotations 
X',  Y',  and  Z'  independently  to  simulate  the  plate  or  shell  element 
characteristics.  The  bending  properties  of  the  plate  and  shell 
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Figure  7.  Variable  Geometry  Compressor  Blade 
Cross  Section. 
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elements  accurately  simulated  the  dynamic  properties  of  the  blade. 
Actual  blade  geometry  data  was  modified  slightly  by  assuming  that 
each  blade  element  had  a  constant  thickness  equivalent  to  the 
average  thickness  of  the  local  area. 

For  the  purpose  of  understanding  the  effect  of 
rotational  restraints  on  the  frequencies,  a  finite  element  model 
with  the  grid  points  having  the  degree  of  freedom  per  node  (0-free 
or  l-fixed)  was  studied.  Several  cases  were  studied  and  the 
results  from  four  cases  are  presented  in  Table  8.  Relaxing  the 
restraint  in  rotation  about  the  Z  axis  for  all  grid  points  up  to 
the  mid-span  and  modeling  a  thick  fillet  at  the  base  brought  the 
analysis  fundamental  frequency  much  closer  to  that  of  the 
experiment.  However,  NASTRAN  runs  performed  by  Pratt  &  Whitney 
Aircraft  Company  (PWA)  and  by  Puterbaugh  (15)  of  AFWAL  did  not 
incorporate  the  base  fillet  and  predicted  a  fundamental  frequency 
which  was  lower  than  the  experimental  value.  The  analysis  by  UDRI 
considered  the  effect  of  a  thick  fillet  at  the  trailing  edge  and 
our  fundamental  frequency  prediction  was  very  close  to  the 
experimental  value.  There  was,  however,  a  variation  in  the 
frequencies  of  higher  modes. 

Table  9  compares  the  results  of  the  three  analyses. 

The  PWA  analyses  illustrated  the  effect  of  the  centrifugal  force 
field  due  to  rotation  on  the  modal  frequencies.  The  flapwise, 
edgewise,  and  torsional  modes  were  affected,  with  the  flapwise 
fundamental  frequency  being  drastically  increased  from  530-540  Hz 
at  zero  speed  to  760  Hz  at  21,500  rpm.  Modal  frequency  changes  due 
to  speed  and  the  comparison  of  frequencies  and  modes  at  zero  speed 
can  also  be  seen  in  the  Campbell  Diagram  shown  in  Figure  9. 

Contour  plots  of  modal  displacements  are  provided  in  Figure  10.  It 
was  rather  simple  to  identify  flapwise,  edgewise,  and  torsional 
modes  from  these  plots.  There  was  good  correlation  of  the  mode 
shapes  furnished  by  the  PWA  with  our  analytical  results.  We  should 
note  that  the  laser  holographic  experiments  did  not  record  the 
edgewise  bending  mode.  This  mode  was  essentially  in  the  plane  of 
the  disk,  and  the  recording  technique  was  not  sensitive  to  small 
inplane  displacements  of  the  edgewise  mode.  The  analytical  effort 


TABLE  8 

SAP  FINITE  ELEMENT  ANALYSIS  OF 
LOW  ASPECT  RATIO  FAN  (LARF)  BLADE 


Case  Degrees  of  Freedom  eigenfrequencies, 

No.  Node  Numbers  X  Y  Z  X"  Y'  Z'  (Hz) 

1  1  through  130  0000  0  1  523 

and  1308 

131  through  143  1111  1  1  1500 

2433 

2628 

3276 

_ 2321 

2  1  through  130  00000  1  662 

and  1348 

131  through  143  1111  1  1  1769 

2480 

2814 

3358 

_ 3767 

3  1  through  77  000001  589 

1320 

78  through  130  0000  0  0  1557 

and  2460 

131  through  143  1111  1  1  2681 

3319 

_ 3687 

4  1  through  77  000000  598 

1325 

78  through  130  000000  1594 

119  0001  1  1  2482 

118  2701 

and  1111  1  1  3333 

_ 131  through  143 _ 1111  1  1 _ 3712 

Key  to  the  Degree  of  Freedom:  0  -  free 

1  -  fixed 
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TABLE  9 

COMPARATIVE  TABLE  OF  THE  NATURAL  FREQUENCIES  (HZ) 
OF  A  LOW  ASPECT  RATIO  BLADE 


Mode 

UDRI 

N=0 

ISRL* 

N=0 

APL  (15) 
N=21500 

PWA-GPD 

N=0 

•  (16) 
N=21500 

SAP  IV 

Case  1 

SAP  IV 

Case  4 

NASTRAN 

(CTRIA  2) 

FROM 

GRAPH 

1 

523 

598 

549 

589 

765 

534 

756 

2 

1308 

1325 

1363 

- 

1501 

- 

1744 

3 

1500 

1594 

1695 

1728 

1851 

1700 

1889 

4 

2433 

2482 

2487 

2049 

2533 

1980 

2649 

5 

2688 

2701 

- 

2718/ 

2763 

- 

2875 

- 

N  -  Speed  in  rpm 

*  Laser  holographic  interferometry  experimental  values 
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Figure  10.  Contour  Plots  for  Mode  Shapes  of  Compressor  Blades 


conducted  by  UDRI  is  based  on  the  information  provided  by  other 
sources.  A  more  ideal  approach  would  be  to  generate  geometrical 
data  of  the  blades  from  a  numerical  control  machine  tape. 

Three  aspects  of  the  analysis  to  define  dynamics  of  fan 
blades  were  investigated.  These  were: 

a.  use  of  triangular  elements  to  investigate  the 
rotational  effects, 

b.  development  of  a  procedure  to  suppress  the 
signularities  in  the  finite-element  analysis,  and 

c.  use  of  quadrilateral  elements  available  with 

NASTRAN. 

2.4.2  TRIA2  Elements 

Blade  geometry  was  modeled  with  triangular  plate 
elements.  Blade  geometry,  mesh  nature,  and  the  axis  of  rotation 
are  illustrated  in  Figure  11.  The  blade  root  fillet  was  modeled 
with  several  elements.  Care  was  exercised  to  avoid  biasing  of 
stiffness  properties  of  triangular  elements  by  orienting  the  longer 
sides  of  elements  to  the  left  and  right  alternatively.  Triangular 
elements  have  both  membrane  and  bending  properties.  Membrane 
properties  of  the  elements  simulated  well  the  in-plane  centrifugal 
force  field.  The  triangular  and  rectangular  elements  with  their 
nodal  degrees  of  freedom  are  shown  in  Figure  12. 

Rigid  format  13  of  COSMIC  NASTRAN  was  used  to  perform 
the  normal  mode  analysis  of  a  rotating  compressor  blade.  Table  10 
furnishes  the  natural  frequencies  for  several  modes  of  a  blade  at 
various  rotational  speeds. 

A  comparison  with  the  results  of  PWA  and  General 
Electric  Company  (GE)  was  made  for  the  fundamental  mode  only.  GE 
used  ANSYS  and  their  fundamental  frequency  at  operating  speed  was 
776  Hz.  PWA  used  MSC  NASTRAN  and  their  frequency  was  673  Hz.  Our 
use  of  COSMIC  NASTRAN  yielded  a  frequency  770  Hz.  Differences 
between  PWA  and  our  results  were  primarily  due  to  small  variations 
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Figure  11.  NAST 


TRIA2  Element  with  15  d.o.f.  (T  =  Translation 
R  =  Rotation) 


Figure  12.  Triangular  and  Rectangular  Plate  Element  Degree 
of  Freedom. 


4 


TABLE  10 

ROTATIONAL  EFFECTS  ON  NATURAL  FREQUENCIES  OF 
LARF  BLADE  WITH  TRIA2  ELEMENT 


Frequency .  Hz 


Hfi. 

(rpro) 

*1 

^2 

^3 

^4 

^5 

1 

0 

550.3 

1,363.9 

1,696.1 

2.492.5 

2,667.7 

2 

6,000 

573.0 

1,379.4 

1,709.2 

2,497.8 

2,671.9 

3 

12,000 

635.9 

1,420.6 

1,749.3 

2,511.9 

— 

4 

18,000 

726.9 

1,478.3 

1,818.5 

2,530.4 

— 

5 

20,500 

770.8 

1,504.7 

1,856.3 

2,537.8 

2,703.3 

in  blade  geometry.  Comparison  of  results  at  higher  modes  or  at 
other  speeds  could  not  be  made  as  their  results  were  not  provided. 

2.4.3  Sinaularities 

Triangular  plate  elements,  as  mentioned  before,  have 
five  degrees  of  freedom  per  node  with  the  out-of-plane  rotational 
constraint  at  the  node  point  of  the  structural  model  confined  to  a 
single  plane.  This  out-of-plane  rotational  constraint  was  easily 
incorporated  in  the  input  as  single  point  constraints  for  all  of 
the  elements.  But  the  compressor  blade  was  severly  twisted  and 
curved,  and  the  assumption  of  out-of-plane  rotational  restraint 
with  respect  to  global  coordinate  systems  was  not  justified. 

Either  the  releasing  or  the  constraining  in  full  of  this  degree  of 
freedom  causes  under  or  over  estimation  of  the  frequency. 

Puterbaugh  (15)  in  his  analysis  earlier  had  underestimated 
frequency  and  also  had  several  singularity  points.  The  grid  points 
that  had  singularities  were  marked  with  an  asterisk,  as  shown  in 
Figure  13.  When  singularities  appeared,  results  were  not  accurate. 
Therefore,  our  approach  was  to  devise  a  scheme  and  apply  a 
multipoint  constraint  method  to  suppress  these  singularity  points. 
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Figure  13.  Grid 
Mode 


The  multipoint  constraint  method  required  the 
determination  of  normal  vectors  to  each  of  the  element  grid  points, 
and  the  direction  cosine  to  the  normal  was  given  as  a  proportional 
and  related  constraint  at  the  points  where  singularities  appeared. 

The  procedure  to  determine  normal  vectors  and  direction 
cosines  was  as  follows: 

Let  A,  B,  and  C  be  the  nodal  points  of  the  TRIA2 
element.  Their  coordinates  w.r.t  global  coordinate  systems  were 

^a'^a'  ^a'^b'  ^b'  ^b  ^c'  ^c'  ^c*  Figure  14  for  element  and 

coordinate  details. 

Vectors  A  and  B  were  given  by 

a  =  (Xj^-x^)i  +  (yjj-y^)j  +  (Zj^-Zg)k  (27) 

and  vector  joining  B  and  C  was  given  by 

^  (28) 

c  =  a  X  b  (29) 


or 


e  = 


b  a 


b  ■'z  b  a 


x_-  y^-  y.  z_  z. 


‘c  "b 


'c  -^b  “c  “b 


(30) 


and  the  unit  vector  along  c  was  given  by 


^  ‘cl 


(31) 
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where  |c|  was  the  magnitude  of  vector  c.  The  direction  cosines  of 

the  unit  or  normal  vector  were  the  coefficients  a,  fi,  -f  of  the  unit 

vector  u  .  a ,  ,  y  were  shown  in  Figure  15. 

o 

Multi-point  constraints  (MFC)  require  that  the 

condition 

U^+  fljUj  +  =  O  (32) 

be  satisfied.  In  the  equation  the  coefficients  a^,  a^,  were  the 
direction  cosines  and  u. ,  u. ,  u.  were  the  d.o.f.  at  each  nodal 

1  ]  K 

point. 

A  computer  program  (VECTOR)  was  written  to  directly 
read  the  geometric  coordinates  of  all  the  grid  points  from  a 
NASTRAN  file  for  TRIA2  elements  and  to  estimate  direction  cosines 
of  unit  vectors  for  all  grid  points.  Relevant  direction  cosine 
components  for  the  singular  points  were  input  through  MFC  cards  as 
described  in  the  NASTRAN  User's  Manual.  With  the  suppression  of 
singularities,  the  analysis  was  rerun  and  the  confidence  in  the 
results  was  restored.  The  method  and  its  application  to  curved 
blade  vibrations  was  presented  in  reference  (17) . 

2.4.4  Quadrilateral  (OUAD2i  Elements 

Superposition  of  four  TRIA2  elements  with  bending  and 
stretching  properties  would  generate  QUAD2  elements  in  NASTRAN.  A 
QUAD2  element  was  illustrated  in  Figure  12b.  It  has  20  d.o.f.  in 
all  with  the  out-of-plane  rotational  restraint  at  nodal  corners. 

The  model  compressor  blade  was  divided  into  120  QUAD2  elements  and 
the  input  file  was  the  same  as  the  one  used  for  SAF  IV.  The  model 
was  evaluated  during  disk  rotation  by  releasing  all  d.o.f.  for  all 
grid  points  up  to  the  midspan  of  the  blade  where  curvature  was 
appreciable.  Rotational  restraint  normal  to  the  surface  was 
retained  for  the  top  half  of  the  blade.  The  results  obtained  from 
this  analysis  were  similar  to  those  of  SAF  IV  and  were  closer  to 
the  experimental  values.  Table  11  provides  a  comparison  of  results 
of  COSMIC  NASTRAN  analysis,  with  TRIA2  and  QUAD2  elements,  at  zero 
and  at  operating  speeds  of  the  compressor.  Experimental  and 
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analytical  results  comparison  from  the  available  sources  is  also 
provided  in  this  table. 


TABLE  11 

COMPARISON  OF  THE  MODAL  FREQUENCIES  (HZ)  OF  A  LARF  BLADE  AS 
DETERMINED  BY  VARIOUS  METHODS 


UDRI 

SAP  IV  COSMIC  NASTRAN  ISRL  PWA 


Mode 

No^ 

Speed  Model  Model 

(roro)  _ 1 _  2  TRIA2 

Model 

1 

Model 

2 

Experi¬ 

ments 

-  MSC 
NASTRAN 

S£ 

ANSYS 

1 

0 

523 

598 

550 

500 

571 

589 

534 

580 

2 

0 

1,308 

1,315 

1,364 

1,178 

1,415 

- 

- 

- 

3 

0 

1,500 

1,594 

1,696 

1,611 

1,657 

1,718 

1,700 

- 

4 

0 

2,433 

2,482 

2,492 

2,321 

2,617 

2,049 

1,980 

- 

5 

0 

2,688 

2,701 

2,668 

2,771 

2,920 

2,718/ 

2,763 

2,875 

- 

1 

20 

,500 

- 

- 

770 

- 

788 

- 

756 

776 

2 

20 

,500 

- 

- 

1,508 

- 

1,511 

- 

1,744 

- 

3 

20 

,500 

- 

- 

1,856 

- 

1,860 

- 

1,889 

- 

4 

20 

,500 

- 

- 

2,538 

- 

2,639 

- 

2,649 

- 

5 

20 

,500 

- 

- 

2,703 

- 

3,035 

- 

- 

- 

NOTES:  NASTRAN  Model  2  has  out-of-plane  rotational  restraint  for 
grid  points  1-77  and  all  d.o.f.  for  grid  points  78-120. 

Models  1  and  2  of  SAP  IV  were  compatible  with  Models  1  and 
2  of  COSMIC  NASTRAN. 

Model  1  has  6  d.o.f.  for  all  grids  except  the  fixed 
boundary . 

2.5.  INTEGRAL  BLADED  DISK  MODELING  AND  ANALYSIS 

The  analysis,  thus  far,  has  considered  the  blades  to  be 
mounted  on  an  infinitely  rigid  disk.  This  assumption  does  not 
accurately 
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represent  the  blade  and  disk  interface  compliance.  Therefore,  our 
analysis  considered  flexibility  of  the  disk  in  analyzing  the  blade 
dynamics . 

When  the  combined  flexibility  of  the  disk  and  blades  was 
considered  as  a  whole,  the  magnitude  of  such  an  analytical  problem 
was  large  in  terms  of  computer  core  memory  and  computational  time. 
Modeling  effort  and  analysis  time  can  be  drastically  reduced  if  the 
cyclic  symmetry  of  the  structure  was  used  in  the  analysis.  Only 
one  blade  and  its  associated  section  of  the  disk  need  be  considered 
for  analysis.  The  photo  in  Figure  16  shows  details  of  the 
axisymmetric  integral  bladed  disk  assembly.  The  disk  has  20  blades 
in  all,  and  the  analysis  was  performed  on  an  18-degree  segment  of 
the  disk.  The  rotational  cyclic  symmetry  feature  available  with 
the  COSMIC  NASTRAN  finite  element  program  was  used. 


The  scheme  to  analyze  the  problem  consisted  of  two  tasks: 
first,  use  of  a  preprocessor  program  (PATRAN)  which  automated  the 
generation  of  the  finite  element  model;  and  second,  operation  of 
the  NASTRAN  finite  element  program  on  the  PATRAN  model  to  define 
resonant  characteristics  of  the  bladed  disk.  Some  details  about 
the  generation  of  the  geometric  and  finite  element  model  and  the 
subsequent  analysis  of  the  disk  sector  are  provided  in  the 
following  paragraphs. 

2.5.1  PATRAN-G  Modeling  of  Bladed-Disk  System 

PATRAN  was  an  interactive  computer  software  program 
(18)  for  analytical  solid  modeling  and  finite  element  pre-  and 
post-processing.  The  procedure  for  bladed-disk  modeling  using 
PATRAN  was  as  follows;  A  set  of  grid  points  called  the  control 
points  was  selected  that  define  the  outlines  of  the  structure  to  be 
modeled.  These  grid  points  were  input  to  the  PATRAN  program  in 
terms  of  their  identifying  numbers  and  their  three  spatial 
coordinates.  The  object  contour  was  generated  from  these  point 
coordinates.  The  material  properties  associated  with  these 
geometrical  entities  were  also  defined  at  this  stage.  This 
completed  analytical  solid  modeling  phase  of  the  disk.  The  solid 
and  thin  lines  in  Figure  17  show  the  analytical  geometry  modeled 


48 


Size  on  Diameter 


All  Dimensions  in  mm 


Size  on  Diameter 


DESIGN  OUTLINE 

ANALYTICAL  MODEL 
OUTLINE 


Figure  17.  Analytical  Model  of  Bladed  Disk  Assembly 


and  the  finite  element  mesh  pattern  to  follow.  Starting  with  the 
control  grid  points  and  connecting  lines  and  surfaces  in  the 
meridional  plane  of  the  blade  and  disk  sections,  PATRAN  generated 
the  three-dimensional  solid  body  shown  in  Figure  18.  Figure  19 
shows  a  hidden  line  plot  of  the  body  as  seen  from  the  inlet  side. 
Since  only  a  curved  surface  representation  was  needed  for  finite 
element  modeling  of  the  blade,  no  solid  volume  was  generated. 
Although  only  an  18-degree  sector  of  the  disk  was  used  in  the 
analysis,  PATRAN  can  be  used  to  develop  the  model  for  the  complete 
20  bladed-disk  system  with  equal  ease. 

The  next  step  in  PATRAN  modeling  was  the  finite  element 
discretization  of  the  analytical  solid  body  created  earlier.  This 
was  done  simply  by  specifying  the  nodal  subdivision  and  connecting 
element  types  for  the  disk  geometrical  entity.  The  surface 
representing  the  blade  was  discretized  using  the  four-node  plate 
element  QUAD.  The  three-dimensional  regions  representing  the  disk 
were  discretized  using  eight-node  HEX  elements.  Six-node  wedge 
elements  were  also  used  where  needed  by  the  geometry.  The  varying 
blade  thickness  was  specified  using  the  data  patch  utility  of 
PATRAN.  A  data  patch  was  defined  by  specifying  the  thickness  at  16 
Gauss  points  of  the  surface  representing  the  blade.  Blade 
thickness  at  element  centroids  was  interpolated  from  this  data 
patch.  Finally,  the  model  was  completed  by  equivalencing  the  nodes 
lying  on  common  boundaries  of  adjacent  connecting  regions  and  then 
resequencing  the  node  numbers  so  as  to  give  an  optimum  bandwidth. 
These  operations  were  also  performed  with  PATRAN.  Figure  20  shows 
the  completed  finite  element  model  of  the  18-degree  sector  of  the 
bladed-disk  system. 

The  above  procedure  created  an  output  data  file  which 
was  not  specific  to  any  particular  finite  element  analysis  program. 
The  data  translator  program  PATCOS  was  then  used  to  convert  the 
model -data  file  to  the  required  COSMIC  NASTRAN  format.  Further 
modifications,  such  as  the  specification  of  multipoint  constraints 
to  enforce  compatibility  between  the  QUAD  and  HEX  elements,  were 
performed  separately.  The  rotational  displacements  of  the  blade  at 
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Figure  19.  Hidden  Line  Plot  of  Bladed  Disk. 


53 


the  blade-disk  interface  were  represented  in  terms  of  the 
translational  displacements  of  the  disk  elements  as: 

®x  ^^p  ■  /  AZ  -  (Wp  -  Wq)  /  AY  =  0  (33) 

By  +  (Wp  -  Wg)  /  AZ  -  (Up  -  Uq)  /  AZ  =  0  (34) 

(^p  ■  /  ay  -  (Vp  -  Vg)  /  AX  =  0  (35) 

where  0  ,  0  ,  and  0  were  the  rotational  displacements  of  the  nodes 

X  ^  z 

at  the  interface;  U,  V,  and  W  were  the  translational  displacements 
of  the  nodes  p  and  q  on  the  disk  in  the  immediate  vicinity  of  the 
interface;  and  AX,  AY,  and  AZ  were  the  distances  between  the  nodes 
p  and  q  in  three  coordinate  directions.  Figures  21a  and  21b  show 
the  rotational  and  translational  degrees  of  freedom  relationship  at 
the  interfacial  grid  points  on  the  blade  and  disk  sections.  The 
assumptions  made  and  various  steps  involved  in  deriving  the  above 
equations  and  their  application  for  bladed  disk  vibration  analysis 
were  presented  in  reference  (19). 

2.5.2  NASTRAN  Finite  Element  Analysis  of  Axisvmmetric  Disks 

Consider  a  circular  disk  of  radius,  a,  and  variable 
thickness,  h(r) ,  and  subjected  to  transverse  shearforce  and  bending 
moment  around  the  radius,  a.  The  governing  differential  equation 
of  motion  of  a  distributed  model  of  the  disk  system  was 

V'‘w(r,0  +  w(r,e)  =  0  (36) 

w(r,f)  was  the  transverse  deflection, 
p  was  the  density  of  disk  material, 

D(r)  was  the  flexural  rigidity, 

^  was  the  angular  coordinate. 

Finite  element  method  discretizes  a  structural  system 
and  the  equation  of  motion  of  a  discretized  bladed  disk  sector  was 

[K]  -  x2[M]  (x)  =  0  (37) 


where 


and 
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where  [K]  was  an  assembled  stiffened  matrix  of  disk  and  blades. 


[M]  was  mass  matrix, 

2 

X  was  frequency  parameter  and  equals  to 


and  (x)  was  a  matrix  of  generalized  coordinates.  From  the  given 
geometry  information  of  grid  points  and  elements  and  the  material 
properties,  the  computer  program  generated  stiffness  and  mass 
matrices,  imposes  multi-point  constraints  as  specified  before  and 
solved  for  eigenvalues  (X  )  of  the  system.  Associated  eigenvectors 
(X)  were  evaluated  and  the  deformed  mode  shapes  were  plotted  using 
PATRAN  post-processing  capabilities. 

The  bladed  part  of  the  segment  was  modeled  with  QUAD2 
elements  with  four  nodal  points,  each  with  six  degrees  of  freedom. 
The  disk  part  of  the  segment  was  modeled  with  eight-node  HEX  solid 
elements  with  three  d.o.f.  per  node.  At  the  interface  of  the  blade 
and  disk,  rotational  degrees  of  freedom  of  QUAD  elements  were 
retained  through  multipoint  constraints  as  specified  in  equations 
(33),  (34),  and  (35). 

The  basics  of  the  cyclic  symmetry  principle  (20)  for 
dynamic  analysis  of  an  18-degree  segment  of  a  bladed  disk  were  as 
follows:  rotational  transformation  between  a  symmetrical 

structural  segment  was  governed  by  the  equation, 

n  -0  ^  vKL  -Kc  ,  - .  „ 

u  =  u  +  2  u  cos(n-l)  Ka 

K=1 

+  u^®  sin(n-l)Ka  +  (-1) (39) 
where  n=l,  2,  . N  =  Number  of  segments. 

K  =  cyclic  index 

Kj^<2  was  the  limit  of  K 

/ 
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-  a  variable  like  displacement,  stress,  etc.  in  the 


nth  segment. 

u  ,  u  ,  u  =  symmetrical  components  (cyclic 
coefficients)  used  in  the  solution  which  defines  the  entire 
structure ,  and 

2'7r 

a  =  —  =  circumferential  angle  for  each  segment. 

In  the  rotational  symmetry  analysis,  the  dependent 
displacements  (side  2  of  the  model,  the  node  points  as  specified 
through  CYJ0IN2  card  in  the  Bulk  Data)  were  related  to  the 
independent  displacements  (Side  1,  as  specified  in  CYJOINl  Card)  by 
the  following  equations  of  constraint  for  the  boundaries. 


u^2  =  u. 


U2^  =  u^*^  cos  Ka  +  sin  Ka 


u^^  =  u^^  sin  Ka  +  u^®  cos  Ka 


ON/2  ^  -N/2 


(40) 

(41) 

(42) 

(43) 


where  K  =  1,  2,  ...  Kj^. 

The  only  symmetrical  components  coupled  by  the  compatibility 
constraints  were  1C  and  IS,  2C  and  2S,  etc.  From  the  compatibility 
constraint  relations,  we  know  that  the  selection  of  each  K  index 
predetermines  the  vibration  modes  which  may  occur  in  the  entire 
disk  model.  In  our  case,  K  =  10  was  identical  to  K  =  0,  since  Ka  = 
25r  and  maximum  K  that  could  be  used  was  10. 

Table  12  presents  the  results  obtained  from  NASTRAN 
analysis  using  various  constraint  conditions.  The  list  of 
frequencies  shows  that  the  fundamental  frequency  of  the  coupled 
system  for  K=0  was  712  Hz  while  the  experimental  value  was  832  Hz. 
Several  blade-alone  frequencies  were  also  evaluated  by  this 
analysis.  Analytic_al  mode  shapes  of  the  system  corresponding  to  the 


displacements  in  the  axial  direction  (parallel  to  rotational  axis) 
are  illustrated  in  Figure  22.  Similar  contour  plots  corresponding 
to  the  displacements  in  the  tangential  direction  were  also 
generated.  Holographic  experiments  have  recorded  the  mode  shapes  of 
the  bladed  disk  and  experimental  blade  mode  shapes  are  shown  in 
Figure  23. 


TABLE  12 

NATURAL  FREQUENCIES  OF  COMPRESSOR  BLADED  DISK 
BY  CYCLIC  SYMMETRY  ANALYSIS 


Mode 

Type 

Boundary 

Conditions 

in 

Analytical 

Model 

Analytical 

Frequency 

(Hz) 

Experimental 

Frequency 

(Hz) 

Umbrella 

6  nodes  on  Flange-bottom 

392 

mode ,  OC 

and  3  nodes  on  disk-face 

883 

fixed.  7.3%  over  estimate 

1279 

in  frequency  of  OC 

1781 

Umbrella 

6  nodes  on  flange-bottom 

391 

588 

mode ,  OC 

fixed.  6.3%  under  estimate 

771 

823 

in  frequency  of  OC 

1275 

- 

1778 

1728 

Blade-alone 

Disk-blade  Interface 

549 

modes 

fixed 

1364 

1696 

2493 

2667 

59 


Figure  22.  Analytical  Mode  Shapes  of  Bladed  Disk- 
Cyclic  Symmetry  Analysis. 
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Mode  4 


d) 


Mode  5 


Mode  6 


Analytical  Mode  Shapes  of  Bladed  Disk- 
Cyclic  Symmetry  Analysis. 


6 


3,627  Hz  3,766  Hz 

Figure  23(cont'd).  Experimental  Mode  Shapes  of  Bladed  Disk. 
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SECTION  III 

PROBABILISTIC  ANALYSIS 


This  chapter  outlines  the  general  approach  used  to  develop  an 
analysis  for  estimating  the  variance  of  structural  response 
variables,  given  the  mean  values  and  variances  of  system  properties 
which  are  probabilistic  in  nature.  The  statistical  analysis 
adopted  was  rudimentary,  to  be  sure;  however,  the  treatment  is 
consistent  with  the  level  of  information  which  is  readily  available 
to  the  engineer  and  lends  itself  to  the  analysis  of  relatively 
large  and  complex  systems.  The  paragraphs  below  discuss  the 
philosophy  of  the  approach,  the  statistical  parameters  of  interest, 
and  the  mathematical  relationships  needed  for  computing  variances 
of  response  guantities  such  as  displacement,  stress,  and  natural 
frequency.  The  details  of  this  development  are  contained  in 
reference  (5) . 

3.1  STRUCTURAL  STATISTICAL  ANALYSIS 


The  notion  of  a  probabilistic  analysis  encompasses  numerous 
possible  analytical  techniques.  Given  that  certain  properties  or 
dimensions  of  a  system  are  subject  to  uncertainty,  the  proper 
choice  of  analysis  method  depends  strongly  upon  the  difficulty  and 
cost  of  a  single  simulation,  and  upon  one's  knowledge  of  the 
statistical  parameters  of  interest.  Some  of  the  possible 
approaches  are  listed  below. 

Stochastic  analysis  involves  stating  the  differential  system 
of  interest  in  terms  of  stochastic  quantities,  and  solving  directly 
for  the  response  in  statistical  terms.  This  approach  is  an  active 
research  area  in  applied  mathematics  (21) .  Solutions  of  one¬ 
dimensional  problems  still  represent  a  formidable  challenge  with 
this  class  of  methods  (22),  and  consideration  of  very  complex 
systems  was  not  feasible  at  this  time. 

Random  field  simulation  is  a  relatively  new  approach  developed 
by  Liu  and  co-workers  (23)  .  In  addition  to  discretizing  the 
deterministic  system  of  interest,  new  unknowns  are  introduced  in  a 
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finite  element  (or  other  numerical)  model  which  describe  higher 
statistical  moments  of  the  response.  This  augmented  problem  is 
solved  in  a  single  step  for  both  the  mean  values  (deterministic 
response)  and  the  additional  statistical  variables.  This  method  is 
capable  of  considering  detailed  autocorrelations  for  the 
statistical  variables,  but  is  best  suited  to  moderately  sized 
systems . 

Monte  Carlo  simulation  is  appropriate  if  the  statistical 
nature  of  the  (few)  independent  variables  is  well-known,  and  the 
cost  of  a  single  simulation  is  small.  Known  information  about  the 
statistical  parameters  is  used  to  generate  a  large  number  of 
samples  with  representative  values,  and  a  deterministic  analysis  is 
performed  for  each  sample.  The  result  is  a  sample  of  the  response 
from  which  statistical  information  can  be  derived  by  standard 
methods  (24) . 

In  the  present  analysis,  we  viewed  probabilistic  properties  of 
a  system  as  discrete  random  variables.  The  elastic  modulus  of  a 
turbine  blade,  for  instance,  might  vary  from  point  to  point  in  a 
different  fashion  for  every  blade  manufactured;  we  chose  to 
characterize  this  modulus  by  a  mean  value  and  a  single  value  of  the 
variance.  The  information  needed  to  perform  a  meaningful 
probabilistic  analysis  with  this  approach  is  usually  available  or 
can  be  estimated  with  a  fair  degree  of  accuracy.  For  example,  if  a 
modulus  value  is  quoted  as  being  "E+AE,"  we  normally  interpret  the 
quantity  AE  as  representative  three  standard  deviations;  the  range 
E+AE  therefore  includes  approximately  99.7  percent  of  all  samples. 

The  discrete  random  variable  approach  required  a  similar  level 
of  information  about  all  statistical  variables.  Therefore,  routine 
quality  control  data  or  manufacturer's  tolerances  are  the  only 
additional  information  needed  beyond  that  used  to  construct  a 
deterministic  finite  element  model.  Together  with  the  relatively 
low  cost  associated  with  solving  relatively  large  models,  this 
simplicity  makes  the  present  method  attractive  for  routine  analysis 
work. 
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Note  that  the  results  of  an  analysis  based  on  the  discrete 
random  variable  approach  are  not  related  in  a  simple  way  to  results 
from  the  alternative  methods  mentioned  previously.  However,  the 
basic  trends  predicted  by  either  method  will  agree;  that  is,  if  the 
dispersion  in  a  particular  variable  is  large,  and  if  the  structural 
response  varies  significantly  with  the  variable  in  question,  then 
we  expect  a  large  variance  in  the  response.  In  some  cases,  it  is 
possible  to  show  that  the  variances  predicted  using  the  present 
method  are  conservative  (overestimated) .  For  example,  the  variance 
in  natural  frequencies  predicted  when  a  physical  property  is 
assumed  to  vary  with  position  is  generally  less  than  that  computed 
when  the  property  is  constant  throughout  the  model,  but  subject  to 
the  same  variation  in  magnitude. 

3.2  STOCHASTIC  ANALYSIS  OF  BLADED  DISK  SYSTEMS 

The  efforts  on  the  probabilistic  finite  element  analysis 
development  were  concentrated  in  three  areas: 

□  steady-state  forced  response  sensitivity; 

□  graphical  output  of  probabilistic  results;  and 

□  example  problem  data  pieparation. 

Activity  in  each  of  these  areas  is  summarized  briefly  below. 

We  completed  a  solution  procedure  for  the  probabilistic 
analysis.  This  path  deals  with  steady-state  forced  vibration  and 
the  associated  sensitivity  analysis.  The  steady-state  analysis 
still  works  with  the  dynamic  equations  of  motion: 

MU  +  KU  =  F(t)  (44) 

However,  the  nodal  forces  (which  result  from  both  concentrated  and 
distributed  forces)  are  assumed  to  vary  harmonically,  that  is: 

F(t)  =  f  sin(«t)  (45) 
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If  the  damping  forces  are  small,  the  response  U(t)  at  large  't'  is 
sinusoidal  as  well,  and  oscillates  at  the  forcing  frequency,  w. 
Therefore  we  assume 


D(t)  =  u  sin(ut)  (46) 

in  which  u  represents  the  steady-state  amplitude  of  the  nodal 
displacements.  The  resulting  system  to  be  solved  is 

(K-w^M)u  =  f  (47) 

Given  the  amplitudes  of  nodal  forces  and  distributed  loads  which 
define  f,  and  the  forcing  frequency,  w,  the  basic  steady-state 
vibration  solution  determines  the  corresponding  nodal  displacement 
amplitudes,  u,  and  the  stress  amplitudes 

a  =  DBu  (48) 

We  were  interested  in  a  range  of  frequencies,  either  in  the 
neighborhood  of  a  resonance  or  a  specified  number  of  excitations 
per  revolution.  The  steady-state  harmonic  solution  accepts  a  list 
of  forcing  frequencies,  and  generates  a  separate  solution  for  each 
one  in  a  single  computer  run. 

Notice  that  for  some  systems  the  response  near  resonance  is  of 
primary  interest,  and  damping  forces  must  be  included  in  the 
vibration  analysis.  This  is  common  in  noise  control  problems,  in 
which  broad-band  excitation  is  the  rule.  For  bladed  disk  systems 
the  effect  of  structural  damping  is  less  pronounced,  since  the 
system  cannot  tolerate  prolonged  excitation  close  to  a  resonant 
frequency. 

The  steady-state  forced  response  sensitivity  solution  does  an 
additional  analysis  at  each  forcing  frequency,  for  each  statistical 
parameter  in  the  system  (modulus,  density,  thickness,  or  geometric 
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variable) .  For  the  parameters  now  being  considered,  it  is 
reasonable  to  expect  that  the  forcing  frequency  is  independent  of 
the  statistical  parameters;  therefore,  the  sensitivities  of  the 
displacement  amplitudes  with  respect  to  any  parameter  'p'  can  be 
obtained  from 


(K-w 


(49) 


The  corresponding  stress  amplitude  sensitivities  are  simply 


ap 


|Obu 

3p 


„3B 


(50) 


in  which  both  u  and  du/dp  are  dependent  upon  forcing  frequency. 


The  steady-state  forced  vibration  sensitivity  branch  has  been 
coded,  run  successfully,  and  validated  using  examples  for  which 
analytical  results  could  be  obtained.  As  with  the  static  and 
natural  frequency  sensitivity  calculations,  this  branch  obtains  a 
solution  quite  economically  because  of  the  ability  to  reuse  the 
already-factored  coefficient  matrix  in  the  sensitivity  solution. 


The  actual  probabilistic  calculations  were  performed  as  post¬ 
processing  operations  because  of  the  wide  variety  of  operations 
which  are  possible.  Recall  that,  for  a  response  quantity  which 
depends  upon  a  number  of  parameters. 


^  ~  ^  (Pj^ r P2 ' '  •  • ' (51) 

the  variance  of  T  is  related  to  the  variances  of  the  statistical 
parameters  by 


Var(7) 


(52) 
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One  interpretation  of  the  probabilistic  results  which  was  both 
simple  and  informative  was  in  terms  of  "percentile"  values.  The 
mean  value,  as  computed  from  the  basic  finite  element  solution,  was 

by  definition  a  50^*^  percentile  value;  that  is,  we  might  expect  the 
actual  value  (stress,  for  example)  to  be  greater  about  half  the 

time,  and  less  about  half  the  time.  A  percentile  value  was 

such  that  the  true  value  of  the  variable  was  less  than  or  equal  to 
Tq  in  Q  percent  of  all  cases.  Let  represent  the  mean  value  of 

our  variable  T,  and  the  standard  deviation  (the  square  root  of 

the  variance) .  If  t  has  a  normal  (or  Gaussian)  probability 
distribution,  the  interval  of  values  represents  about 

68  percent  of  all  possible  values  of  T.  The  interval  (r  ,T  +t  ) 

H  fi  o 

therefore  includes  approximately  34  percent  of  all  possible  values 

of  T ;  this  means  that  the  value  of  T  will  be  less  than  T  +t  84 

M  o 

percent  of  the  time.  That  is,  the  value  t+t  was  the  84^^ 

fM  V/ 

percentile  value  of  t .  The  percentile  value  also  can  be  viewed  as 
a  figure  of  reliability  or  confidence  level.  The  changes  between 
percentile  values  provide  a  direct  indication  of  the  relative 
uncertainty  in  a  particular  response  quantity. 

3.3  SENSITIVITY  ANALYSIS  OF  BLADED  DISK  SYSTEMS 


The  calculations  necessary  for  computing  both  a  fundamental 
solution  and  sensitivity  data  for  the  natural  frequency  problem 
were  implemented  and  verified  for  a  number  of  examples.  In  the 
following,  we  summarize  the  primary  computational  procedures  and 
present  a  couple  of  typical  examples. 

3.3.1  Element  Stiffness  Derivatives 

Consider  a  problem  in  which  the  stiffness  and  mass 
properties  of  a  structure  depend  on  a  collection  of  independent 
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statistical  parameters.  For  the  purpose  of  computing  response 
derivatives,  we  denote  a  typical  one  of  these  parameters  by  p,  and 
adopt  the  following  notation  for  parametric  derivatives: 

(.)'  =  3(.)/5p  (53) 

While  it  was  possible  to  compute  K'  directly  for  an 
element  and  then  obtain  the  product  K'U,  this  approach  was 
unnecessarily  time  consuming.  We  prefer  to  form  the  product  K'U 
directly  in  vector  form,  on  an  element-by-element  basis. 

Let  the  stiffness  matrix  for  an  element  be  given  by 
K  =  f  A^b'^DBA  IjI  dll,  (54) 

h, 

in  which  A  was  a  transformation  from  local  to  global  coordinates, 
u=AD,  B  was  a  strain-displacement  matrix,  and  D  was  the  elasticity 
matrix.  The  region  was  the  domain  of  the  element  in  parametric 

coordinates.  The  transformation  matrix  A  may  vary  from  point  to 
point  for  curvilinear  elements,  but  was  constant  over  an  element  in 
many  simpler  elements. 

3. 3. 1.1  Intrinsic  Parameters 

The  product  K'U  was  simplest  when  parameter  p 
corresponds  to  an  intrinsic  property,  such  as  the  modulus  or 
thickness,  since  only  the  elasticity  matrix  was  affected.  Noting 
that  AU=u,  the  local  displacement  vector,  we  can  compute 

a  =  D'Bu  (55) 

and 

K'U  =  f  dn^  (56) 
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The  product  in  equation  (55)  was  analogous  to  the  usual  process  of 
stress  recovery,  so  that  the  calculation  of  K'U  for  an  element 
resembles  an  evaluation  of  the  internal  nodal  forces. 

3 . 3 . 1 . 2  Geometric  Parameters 

When  the  parameter  of  interest  affects  the 
nodal  positions,  nonzero  derivatives  may  occur  for  the  element  of 
area  1j|,  the  strain-displacement  matrix  B,  and  possibly  for  the 
coordinate  transformation  matrix  A.  We  assumed  that  the 
derivatives  of  the  nodal  coordinates  were  known,  and  represented 
these  by  X*.=dX.  /dp,  in  which  i  ranges  from  one  to  three,  and  K 

Xil\  Xll\ 

from  one  to  the  number  of  nodes  in  the  element.  The  calculation  of 
B'  and  |j|'  depends  primarily  upon  the  computation  of  sensitivities 
for  the  shape  function  derivatives,  3 (N^  i)/dp.  We  can  show  that 

K,  1 


dp  dp 


(57) 


Notice  that  x^^  denote  the  coordinates  in  the  local  (element) 

system.  Summation  was  implied  over  the  appropriate  range  of  all 
repeated  indices.  Since  the  derivatives  of  the  coordinate  trans¬ 
formation  will  be  considered  explicitly  later,  we  can  hold  A 
constant  to  give 


mJ 


=  A  X' 
mn  nJ 


(58) 


It  was  straightforward  to  show  that  the  product  in  equation  (57) 
simplifiec.  to 


N  N  .x' 
K,m  J, i  mJ 


(59) 


The  derivative  B'  may  be  formed  directly  from  N'  . ,  in  the  same  way 

K,  1 

the  original  B  matrix  was  obtained  from  the  shape  function 
derivatives  N„  . . 

K,  i 
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Similarly,  since; 


^  =  Ul  N 


K,i 


(60) 


we  obtain  for  the  parameter  derivative  of  the  Jacobian 

|j1'  =  |j1  (61) 

The  derivative  of  the  transformation  A  can  be  obtained  as  a  by¬ 
product  of  computing  A  itself,  with  a  minimal  amount  of  additional 
computation. 


For  geometric  control  parameters,  then,  the 
product  K'D  is,  in  general 


+  a'^b'^C  (a+a)  I  j|  +  ajj]']  }  dO^ 

in  which 


a  = 

DBu 

(63a) 

a  = 

DBu 

(63b) 

a  = 

DB'U 

(63c) 

u  = 

AU 

(63d) 

u  = 

A'D 

(63e) 

We  have  verified  the  computational  algorithms  used  to  implement 
(62)  and  (63)  with  a  number  of  examples  for  which  the  results  can 
be  evaluated  in  closed  form  and  for  additional  cases  using  finite 
difference  approximations  for  the  derivatives.  Isoparametric 
finite  element  geometric  sensitivity  analysis  is  described  in 
reference  (25) . 
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3.3.2 


Suppose  that  a  solution  has  been  performed  for  several 
of  the  dominant  modes  of  a  system 

KD  -  W^MD  =  0  (64) 

Differentiating  (64)  with  respect  to  the  parameter  of  interest 
leads  to  the  frequency  sensitivity  expression 

dT(K'-w?M')D. 

u?  =  — ^ -  (i  not  summed)  (65) 

2u^Dj^llD^ 

for  the  i^^  mode  of  vibration.  Equation  (65)  remains  valid  when 
repeated  roots  were  present,  and  for  any  method  of  normalizing  the 
eigenvectors  The  denominator  was  a  scalar  multiple  of  the 

generalized  mass  for  mode  i,  which  we  choose  to  evaluate  at  the 

T 

system  level.  The  product  may  be  computed  element  by 

element,  using  the  procedure  outlined  previously.  We  discuss  the 
evaluation  of  the  vector  M'D^  below. 

Letting  v,V  denote  a  particular  component  of  the 
element  displacement  vector  in  local  and  global  axes,  respectively, 
we  may  write  the  contribution  to  the  mass  matrix  for  component  V 
as 


=  f  |j|  dO  (66) 

■'fig 

in  which  fi  was  a  function  of  the  element  density  and  thickness. 

The  best  procedure  for  the  sensitivity  calculation  in  this  case  was 

element  dependent.  However,  the  fact  that  h'^AV  =  v(x),  the 
pointwise  value  of  v,  can  always  be  exploited.  Similarly,  the 
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product  N  A'D  resembles  a  point  displacement  value,  but  without  the 
same  physical  interpretation.  Again,  the  basic  sensitivity 
calculations  were  limited  to  A'  and  | j| . 

3.3.3  Solution  Procedure 


This  section  outlines  computational  procedures  used  in 
evaluating  response  sensitivitiei,  which  were  needed  for  later  use 
in  the  probabilistic  calculations.  We  assume  that  element-level 
routines  were  available  for  evaluating  the  vectors  K'U  and  M'U,  and 

T  T 

the  scalar  products  D  K'D  and  D  M'D  as  required.  To  implement 
these  methods,  we  perform  element  calculations  for  a  number  of  load 
cases  or  modes  and  a  number  of  sensitivity  parameters,  all  in 
parallel.  Sensitivity  parameters  may  include  the  material  modulus 
or  density,  element  thickness,  and  any  geometric  control  parameter 
defined  in  terms  of  derivatives  of  the  global  Cartesian  coordinates 
at  selected  nodes  with  respect  to  the  parameter.  We  first  review 
briefly  a  solution  procedure  for  static  problems  to  help  fix 
ideas . 


3 . 3 . 3 . 1  Static  Response  Sensitivity 

In  static  analysis,  we  first  factor  the 
original  stiffness  and  solve  for  the  nodal  displacements: 

K  =  LDl"^  (67a) 

LDl'^D  =  F  (67b) 

For  the  first  pass  of  sensitivity  calculation,  form  the  right-hand 
side  and  solve  for  displacement  sensitivities: 

^el 

R  =  F'  -  X  (K'U)  (68a) 

e=l 


U)l7xJ*  =  R 


(68b) 
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The  second  pass  of  element  sensitivity  calculations  yields  the 
element  stress  sensitivities: 


a'  =  (DBA) 'D  -f  DBAU'  (69) 

Which  of  the  matrices  D,  B,  A  possesses  nonzero  derivatives  was  a 
function  of  parameter  type. 

Notice  that  if  a  particular  sensitivity 
parameter  does  no;  affect  a  given  element  directly,  the  calculation 
of  K'D  may  be  skipped,  and  the  stress  sensitivity  reduces  to  a'  = 
DBAU'.  In  practice  it  was  convenient  to  maintain  a  list  of 
switches  for  each  element,  indicating  the  status  (active  or 
inactive)  of  all  parameters.  The  selection  of  parameters  such  as 
modulus,  density,  and  thickness  may  be  tied  to  material  or  property 
set  numbers,  making  it  easy  to  determine  whether  a  specific  element 
was  affected.  If  geometric  parameters  were  defined  in  terms  of 
nodal  coordinate  derivatives,  the  parameter  was  inactive  for  a 
given  element  only  if  all  derivatives  for  each  node  connected  to 
the  element  were  zero. 

3. 3. 3. 2  Frequency  and  Mode  Shape  Sensitivity 

For  eigenvalue  problems,  we  first  solve  the 
eigensystem  and  compute  a  generalized  mass  for  each  mode: 

KO^  -  u?MD^  =  0  (70a) 

(i  not  summed) 

ro^  =  dTmd^  (70b) 

For  each  parameter  and  mode,  the  frequency  sensitivity  equation 
(65)  may  be  summed  element  by  element: 

^el 

"i  =  2^;;^  ^  t  ^  summed)  (71) 

i  i  6=1 
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In  computing  sensitivities  of  the 


eigenvectors,  we  adopt  a  modal  representation.  For  the  mode, 
let  the  eigenvector  derivative  be 

Uf  =  (72) 

in  which 

♦  =  [  D^,  U2,  ...,  ]  (73) 

is  the  modal  matrix,  and  was  a  vector  of  modal  participation 

factors.  Introducing  (72)  into  the  derivative  of  the  original 

T  th 

eigenvalue  equation,  and  preraultiplying  by  ♦  gives,  for  the  i 
mode 


(k-w^m)^.  = 


-♦'^(K'-W?M')U^  + 


2W MD. 
11  1 


(i  not  summed)  (74) 


T  T 

Here  k  =  #"K*  and  m  =  ♦  M#  were  the  diagonal  generalized  stiffness 
and  mass  matrices.  Note  that  only  the  i^^  component  of  the  product 
•  MU^,  which  was  a  column  of  m,  was  nonzero.  The  element  of 
corresponding  to  mode  'n'  was  therefore 


(k  -w?m  ) 
'  nn  1  nn' 


(i-^n; 


i,n  not  summed) 


(75) 


Let  J  be  the  degree  of  freedom  which  attains  the  largest  value  for 
mode  i;  that  is 


(Ui)j  =  sup  (U.)^  (76) 

n 

We  force  the  normalizing  basis  for  mode  i  to  remain  constant  by 
requiring  (U?)j  =  0.  This  condition  was  sufficient  to  determine 

the  remaining  element  of 
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1 


(77) 


(B.)  .  = 

vk3_/  1 


(Hi), 


n=i 

n^i 


in  which  N  was  the  number  of  modes  retained  for  the  sensitivity 
solution.  The  necessary  products  besides  the  system  generalized 

.  T  T 

stiffness  and  mass  consist  of  D^K'D.  and  D  M'D. ,  which  may  be 

n  1  n  1 

evaluated  on  an  element  basis. 

3.3.4  Numerical  Examples 

Some  typical  validation  problems  used  in  verifying  both 
the  static  and  natural  frequency  sensitivity  solutions  were 
summarized  briefly  below.  The  validation  problems  solved  to  date 
include  simple  problems  with  analytical  solutions,  as  well  as 
complex  cases  in  which  computed  sensitivities  were  compared  with 
approximate  parameter  derivatives  obtained  with  finite 
differences  (26) . 

Example  1.  Bending  of  a  Cantilever  Beam 

Figure  24  shows  a  cantilever  beam  subjected  to  a 
transverse  force  at  the  tip.  Again,  an  analytical  solution  was 
possible  both  for  the  displacement  and  rotation,  and  for  the 
sensitivities  with  respect  to  modulus,  E,  thickness,  t,  and  width, 
b.  Five  bilinear  Mindlin  plate  elements  (27,28)  were  used  to  model 
the  beam;  this  number  was  sufficient  for  good  accuracy  but,  since 
the  elements  have  only  linear  displacement  and  rotation  fields, 
does  not  reproduce  the  exact  solution.  Table  13  summarizes  the 
computed  results  for  the  displacements  and  rotations.  Note  that 
the  displacement  sensitivities  were  no  less  accurate  than  the 
displacements  themselves  (all  were  approximately  1  percent  in 
error) ,  and  that  the  rotational  results  were  exact.  Table  14  shows 
moment  and  shear  results,  and  the  force  sensitivities  which  were 
nonzero.  In  all  cases,  the  moment  and  shear  sensitivities  were 
exact,  despite  small  errors  in  the  displacement  solution. 
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TABLE  13 

DISPLACEMENT  SENSITIVITY  DATA  FOR  CANTILEVER  BEAM^ 


Ouantitv 

Exact 

Exact  fx=L) 

Comouted 

w 

2P(3Lx^-x^)/Ebt^ 

-0.4 

-0.39602 

dw/dE 

-2P(3Lx^-x^)/E^bt^ 

-4.0X10“® 

-3.96x10"® 

dw/dt 

-6P(3Lx^-x^)/Ebt^ 

-12.0 

-11.88 

dw/dh 

-2P(3Lx^-x^)/Eb^t^ 

• 

o 

1 

-0.39602 

e 

6P(2LX-X^)/Ebt^ 

0.06 

0.06 

dd/dE 

-6P(2LX-X^)/E^bt^ 

-6.0X10"^ 

-6.0x10"® 

d$/dt 

-18P(2Lx-x^)/Ebt^ 

-1.8 

-1.8 

de/dh 

-6P(2LX-X^)/Eb^t^ 

VO 

o 

• 

o 

1 

VO 

o 

• 

o 

1 

^E=lxlO^;  i/=0;  t=0.1;  b=l;  L=10;  P=1 


TABLE  14 

FORCE  SENSITIVITY  DATA  FOR  CANTILEVER  BEAM 


t 


Element  Centers 


Ouantitv 

X=1 

x=3 

■Hi 

X=7 

X=9 

M  Exact 

-9.0 

-7.0 

-5.0 

-3.0 

-1.0 

Comp. 

-9.0 

-7.0 

-5.0 

-3.0 

-1.0 

dM/db  Exact 

9.0 

7.0 

5.0 

3.0 

1.0 

Comp. 

9.0 

7.0 

5.0 

3.0 

1.0 

Q  Exact 

1.0 

1.0 

1.0 

1.0 

1.0 

Comp. 

1.0 

1.0 

1.0 

1.0 

1.0 

dQ/db  Exact 

-1.0 

-1.0 

-1.0 

-1.0 

-1.0 

Comp. 

-1.0 

-1.0 

-1.0 

-1.0 

-1.0 

t  7 

'E=1xio  ;  »/=0;  t=0.1; 

b=l; 

L=10;  P=1 

Only  nonzero  values 

shown ; 

in  in 

dE  at 

aE  at 

=  0  : 

identically 
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For  the  cantilever  beam  (Figure  24) ,  an  analytical 
solution  for  parameter  sensitivities  of  the  natural  frequencies  was 
quite  simple.  Defining 


^  = 


y 


y  12pL 


The  bending  frequencies  were  where  or^  was  independent  of 

the  geometry  and  properties  of  the  beam.  In  particular,  the  first 
three  frequencies  have  a  =  3.52,  22.0,  and  61.7,  respectively. 

Table  15  summarizes  the  results  obtained  for  a 
particular  case,  using  three  different  meshes.  The  sensitivity 
parameters  were  modulus,  E,  density,  p,  thickness,  t,  width,  b,  and 
length,  L.  It  was  instructive  to  study  the  results  from  a 
relatively  coarse  model  (five  bilinear  elements)  first;  this  model 
was  labeled  Mesh  1  in  the  table.  All  computed  results  for  the 
first  mode  were  quite  good  for  Mesh  1,  with  the  error  in  frequency 
sensitivities  being  similar  in  magnitude  to  the  frequency  error 
itself.  For  the  next  two  modes,  the  sensitivities  for  intrinsic 
parameters  E,  p,  and  t  were  at  least  equal  in  accuracy  to  the 
frequencies,  for  reasons  explained  in  the  last  example.  The  length 
sensitivity  in  Mesh  1  has  been  defined  by  attributing  coordinate 
sensitivities  only  to  the  end  nodes,  however,  and  was  rather  poor: 
this  "local"  geometry  parameter  does  not  enter  the  finite  element 
frequency  equation  in  the  same  manner  as  in  the  analytical 
solution. 


Meshes  2  and  3  represent  the  two  obvious  solutions  to 
the  poor  accuracy  of  Mesh  1  for  the  length  parameter  in  higher 
modes.  Mesh  2  was  a  refined  model,  in  which  ten  elements  were 
used,  and  the  coordinate  sensitivities  were  defined  for  the  end 
nodes  only,  as  in  Mesh  1.  All  results  were  much  improved,  as 
expected;  but  the  du/dh  sensitivity  was  still  less  accurate  (15.9 


percent  error  for  the  third  mode)  than  the  frequency  (7.2  percent 
error) . 


Mesh  3  was  different  from  Mesh  2  only  in  the 
specification  of  the  nodal  coordinate  sensitivities,  which  now  were 
specified  so  that  all  nodes  move  proportionally  when  the  length 
changes.  Mesh  3  produces  results  which  were  essentially  exact  for 
Mode  1  and  reduces  the  error  in  du/dL  by  half  for  the  higher  modes. 
For  Mesh  3,  the  error  in  all  of  the  sensitivity  results  was 
generally  no  larger  than  the  frequency  error  for  each  case 
considered. 

Example  3.  Twisted  Plate  Vibration 

Figure  25  shows  a  twisted  cantilever  plate  which  has 
been  used  extensively  for  the  comparison  of  natural  frequency 
predictions  by  MacBain,  Kielb,  and  Leissa  (29) .  We  wish  to  compare 
natural  frequency  sensitivities  obtained  with  the  present  analysis 
to  those  derived  from  finite  differencing.  For  the  case 

.  7 

considered,  we  take  E=10  ,  »/=0.30,  p=0. 00026;  the  dimensions  were 
length  a=3,  width  b=l,  and  thickness  h=0.050,  all  constants  being 
in  inch  units. 

For  the  present  analysis,  we  employ  a  6x6  mesh  of 
bilinear  elements,  which  was  adequate  for  the  first  few  modes.  As 
evidence  of  this.  Table  16  summarizes  the  first  several  modes  pre¬ 
dicted  for  a  twist  angle  of  9=30*.  Nondimensional  frequencies 
obtained  from  NASTRAN  using  a  mesh  of  128  TRIA2  elements  were  tab¬ 
ulated  as  well,  and  the  two  solutions  were  in  reasonable  agreement. 

Table  17  shows  the  computed  frequency  sensitivities  for 
modes  1-4,  for  a  twist  angle  of  32*.  Finite  difference  estimates 
for  the  sensitivities  were  obtained  from  separate  eigenvalue 
solutions  performed  for  twist  angles  of  31.9*  and  32.1*.  The 
agreement  of  the  predictions  was  reasonably  good  and  was  quite 
accurate  where  the  sensitivity  was  large  in  magnitude. 
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TABLE  15 

FREQUENCY  SENSITIVITIES  FOR  CANTILEVER  BEAM^ 


Mode 

(i) 

du/dt 

8u/ah 

du/3L 

1 

Exact 

201.6 

l.Olxlo"^ 

-3.97X10^ 

2016. 

0. 

-40.32 

Mesh-1 

202.3 

l.Olxlo"^ 

-3.98X10^ 

2023. 

3.9xl0“^° 

-40.67 

Mesh-2 

201.6 

l.Olxlo"^ 

-3.97x10^ 

2016. 

4.9xl0"^° 

-40.37 

Mesh- 3 

201.6 

1.01X10“^ 

-3.97x10^ 

2016. 

4.9xl0“^° 

-40.32 

2 

Exact 

1260.1 

6.30X10”^ 

-2.48X10^ 

12601. 

0. 

-252.03 

Mesh-1 

1391.9 

6.96xio“^ 

-2.74X10® 

13904. 

3.8xl0“^° 

-310.24 

Mesh-2 

1292.9 

6.46X10“^ 

-2.55x10® 

12917. 

4.6xl0“^^ 

-265.25 

Mesh-3 

1292.9 

6.46X10“^ 

-2.55x10® 

12917. 

4.6xl0”^^ 

-258.46 

3 

Exact 

3534.1 

1.77X10"'^ 

-6.96X10® 

35341. 

0. 

-706.82 

Mesh-1 

4727.9 

2.36X10""^ 

-9.31X10® 

47120. 

1.5xl0"^° 

-1260.49 

Mesh-2 

3790.1 

1.90X10"^ 

-7.46x10® 

37809. 

-11 

5.8X10  ^ 

-818.98 

Mesh-3 

3790.1 

1.90X10’^ 

-7.46X10® 

37809. 

-11 

5.8x10 

-757.10 

=1x10^;  v=0;  t=0,l;  p=0. 000254;  L=10;  b=l 

TABI.E  16 

FREQUENCY  PARAMETER  COMPARISON  FOR  3 O’ 

TWISTED  PLATE ^ 

Mode 

Tvoe 

Present 

1 

1-B 

3.42 

3.20 

2 

2-B 

19.10 

19.08 

3 

1-T 

26.04 

26.52 

4 

3-B 

60.15 

61.62 

5 

1-EB 

73.00 

74.19 

6 

2-T 

78.50 

82.93 

t  2  - 

'Normalized  frequencies  were  X  =  ua  Jph/D 
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TABLE  17 

FREQUENCY  SENSITIVITIES  FOR  32“  TWISTED  PLATE 


Mode 

Tvoe 

ur31.9“^ 

0^32*^ 

wr32.1“) 

d<i)/d$  ^ 

Lu>/L0 

1 

1-B 

921.16 

920.23 

919.29 

-467.3 

-535.1 

2 

2-B 

5354.91 

5346.21 

5337.53 

-4668.9 

-4979.0 

3 

1-T 

9006.48 

9017.26 

9028.02 

4927.2 

6170.8 

4 

3-B 

17141.7 

17119.2 

17096.8 

-12791.5 

-12862.9 

^Note  that  the  sensitivity  variable  0  was  defined  in  radians. 


3.4  PROBABILISTIC  FINITE  ELEMENT  CODE  AND  ITS  APPLICATION  TO 

BLADED  DISK  SENSITIVITY  AND  FREQUENCY  RESPONSE  ANALYSIS 

The  finite  element  sensitivity  and  probabilistic  analysis 
techniques  were  implemented  in  a  computer  program  called  PROTEC 
(Probabilistic  Response  of  Turbine  Engine  Components) .  PROTEC  was 
written  in  ANSI  FORTRAN  77  and  has  been  executed  on  CDC  Cyber, 
CRAY/XMP,  and  DEC  VAX  machines.  The  input  to  PROTEC  was  arranged 
in  a  series  of  input  blocks  such  as  boundary  conditions,  nodal 
coordiantes,  coordinate  derivative  data,  element  connections,  nodal 
force-moments,  material  properties,  element  geometry  properties, 
etc.  Data  translation  to  and  from  PATRAN  was  performed  by  two 
intereface  programs  -  PATPRO  (PATRAN-TO-PROTEC)  and  PROPAT  (PROTEC- 
TO-PATRAN) .  PATPRO  converted  a  finite  element  neutral  file  from 
the  geometric  modeling  program  PATRAN  into  a  standard  input  file 
for  finite  element  analysis  by  PROTEC.  PROPAT  transformed  a  PROTEC 
results  file  into  a  PATRAN  results  file  for  post-processing. 

The  modeling-analysis-post  processing  cycle  began  in  PATRAN 
where  the  finite  element  geometric  model  was  generated.  PROTEC 
performed  vibration  and  sensitivity  analysis  from  the  information 
generated  and  received  from  PATPRO  code.  The  results  file  POSFIL 
was  then  processed  by  PROPAT  to  make  compatible  plotting  data  for 
PATRAN.  PRODIS  (PRTEC-TO-DISPLA)  was  an  output  processor  for 
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PRTOEC  that  performed  probabilistic  computations  and  presented 
graphics  using  the  DISSPLA  library.  Two  modes  of  presentation  were 
included  in  PRODIS  for  display  of  probabilistic  data.  Nodal  or 
elemental  variance  data  from  static  or  natural  vibration  solution 
could  be  displayed  in  histogram  form  to  compare  different 
parametric  effects  from  different  analysis  cases.  Also,  the 
results  from  steady-state  harmonic  analysis  with  forcing  frequency 
as  an  independent  variable  were  typically  presented  as  x-y  plots  or 
3-D  surfaces. 

A  concise  general  procedure  flowchart  for  the  stochastic 
analysis  was  shown  in  Figure  26.  A  38-bladed  disk  designed  for  the 
United  Technology  Research  Center  NSMS  demonstration  test  was 
chosen  for  an  analysis  model.  The  finite  element  model  of  the 
tuned  system  was  created  with  PATRAN  and  was  shown  in  Figure  27. 
Through  PATPRO,  the  model  data  were  input  to  PROTEC  wherein  the 
stress,  deformation,  and  sensitivity  information  at  the  natural 
frequencies  of  the  model  were  calculated. 

Since  the  conventional  stress  and  deformation  results  were  not 
the  focal  point  in  this  effort,  the  emphasis  was  placed  on  the 
sensitivity  analysis.  The  natural  frequencies  of  the  tuned  disk 
and  the  frequency  variations  caused  by  mistuning  characteristics 
were  tabulated  in  Tables  18  and  19,  respectively.  The  variations 
were  plotted  in  Figure  28. 

The  values  in  Table  19  were  obtained  from  a  probabilistic 
equation  that  has  the  following  form: 

Var.  [w]  =  Var.(P).  (79) 

where  P  was  the  sensitivity  parameter  (i.e.,  modulus,  thickness, 
length)  and  a  was  the  natural  frequency.  Moreover,  the  variance  of 
the  sensitivity  parameter  can  be  obtained  by  a  statistical 
computation  (e.g.,  Monte  Carlo  simulation)  or  an  analytical 
approach  (e.g.,  micromechanics). 
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Figure  26.  PROTEC  Solution  Flow  Chart. 
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Figure  27.  PATRAN  Generated  38  Bladed  Disk  FE  Model. 


To  compare  the  sensitivity  parameters  with  each  other,  they 
can  be  represented  by: 

S.V.  =  (|^)  P  (80) 

where  P  was  the  sensitivity  parameter  and  S.V.  was  the  sensitivity 
value,  an  index  for  comparison. 

TABLE  18 

NATURAL  FREQUENCIES  OF  38-BLADED  DISK  -  PROTEC  ANALYSIS 
Disk  Geometry 

Overall  Diameter:  11  inch 
Clamped  Inner  Diameter:  2  inch 
Blade  Length:  1  1/2  inch 
Disk  Thickness:  1/16  inch 
Material:  Steel 


MODE 

NASTRAN 

PROTEC 

HOLOGRAPHIC 

EXPERIMENT 

ID 

72.2/72.2 

77.8/77.8 

75 

OC 

85.8 

82.8 

86 

2D 

103,3/103.3 

102.8/102.8 

106 

3D 

195.0/195.0 

184.1/184.1 

186 

4D 

294.1/294.1 

— 

287 

TABLE  19 

SENSITIVITY  VALUES  FOR  BLADE  MODULUS  (E)  , 

BLADE  THICKNESS  (T) ,  AND  BLADE  LENGTH  (L)  VARIATIONS 


Mode 

svrEi 

SV(T) 

SV^L) 

1 

0.030 

5.48 

9.0 

2 

0.062 

3.30 

6.0 

3 

0.242 

8.19 

10.5 
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Based  upon  the  results  shown  in  Table  19  and  Figure  28,  the 
variation  of  the  blade's  length  had  the  largest  "influence”  in  the 
natural  frequency  sensitivity.  Blade  thickness  had  the  next 
largest  effect  while  modulus  variations  had  only  minor  effect.  In 
other  words,  the  research  focus  should  be  placed  more  upon  the 
shape  imperfections  than  on  the  modulus  variations  of  a  mistuned 
bladed  disk. 

To  understand  the  harmonic  response  characteristics  of  a  38- 
bladed  disk  system  (UTRC  model) ,  three  harmonic  loading  cases  were 
analyzed.  As  shown  in  Figure  29,  case  A  was  a  single  point  load, 
case  B  was  a  double  point  load  with  equal  magnitude  but  opposite 
sign,  and  case  C  was  a  four  point  load  with  ecpial  pairs  of 
magnitude  and  opposite  sign.  The  objective  of  the  harmonic 
analysis  was  to  determine  the  response  behavior  of  the  disk  as  it 
transited  through  resonance. 

The  first  six  natural  frequencies  of  the  bladed  disk  system, 
presented  in  Table  18  and  shown  in  Figure  30,  were  calculated  with 
PROTEC.  The  chosen  frequency  range  of  these  six  resonances  was 
from  zero  to  300  Hz.  The  harmonic  response  calculations  were 
performed  within  this  frequency  range.  Finer  frequency  increments 
were  used  in  the  neighborhood  of  the  resonances.  Case  A  with  the 
single  point  applied  load  had  the  response  characteristics  of  the 
free  vibration  resonances,  as  we  might  expect.  Case  B  had  three 
resonance  responses  within  the  frequency  range  stated  before,  while 
case  C  had  only  two  resonance  responses.  Cases  B  and  C  could  not 
excite  the  OC  (Umbrella  Mode)  response  because  of  their  bipolar 
loading  conditions  and  could  excite  only  the  odd  or  even  modes 
because  of  their  geometries. 

Figures  31  and  32  show  the  odd  and  even  node  stress 
distributions  respectively  on  the  bladed  disk  system  due  to  the 
three  different  harmonic  loads.  The  stress  comparison  between  the 
single  point  load  system  and  the  double  point  load  system  was 
depicted  in  Figure  32.  Since  the  double  point  load  system  only  had 
three  resonance  responses  within  the  frequency  range  analyzed,  the 
stress  comparison  with  the  single  point  load  system  was  limited  to 
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Figure  29.  Harmonic  Loading  on  Bladed  Disk. 


three  corresponding  natural  frequencies.  A  similar  stress 
comparison  between  the  single  point  load  case  and  the  four  point 
load  case  was  shown  in  Figure  32. 
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SECTION  IV 

PLANE-OF-LIGHT  PROBE  SYSTEM  DEVELOPMENT 


In  the  past,  measurement  methods  involving  strain  gage-slip 
ring,  strain  gage-telemetry,  and  laser  interferometry  systems  have 
been  used  to  monitor  blade  vibrations.  These  methods  have 
disadvantages  like  (a)  limited  survivability  in  high  temperature 
and  high  speed  environments,  (b)  required  disassembly  of  parts  for 
sensor  installation  and  repair,  (c)  location  of  sensors  and 
transmitting  elements  on  the  rotor,  and  (d)  nonaccessibility  to 
viewing  optics  except  for  the  front  stage  of  the  compressor.  The 
drawbacks  associated  with  the  strain  gage  and  interferometry 
methods  have  led  to  the  development  of  noncontacting  deflection 
measurement  systems  for  rotating  turbine/compressor  blades.  This 
section  describes  the  concept  and  the  optical  and  electronic 
development  of  a  four-channel  fiber  optic  plane-of-light  sensor 
method,  also  referred  as  the  noninterference  stress  measurement 
system  (NSMS)  in  this  report,  to  monitor  blade  tip  vibrations  of 
rotating  turbine  disks. 

4.1  NONINTRUSIVE  VIBRATION  MEASUREMENT  METHODS  -  RELEVANCE  AND 
PAST  DEVELOPMENT 


The  idea  of  measuring  compressor  blade  vibrations  by  a  remote 
optical  sensor  was  proposed  and  implemented  by  Hohenberg  (30)  and 
Zablotskiy,  et  al.,  (31,32).  Nieberding  and  Pollack  (33)  of  NASA 
Lewis  and  McCarty  and  Thompson  of  AEDC  developed  (34)  and  measured 
(35)  the  compressor  blade  tip  vibrations  by  a  single  spot  sensor. 
Roth  (36)  and  Endoh  (37)  developed  multiprobe  spot  sensor  methods 
to  measure  rotating  blade  vibrations.  Spot  sensors  lack  the 
ability  to  always  measure  the  same  location  on  the  blade  tip  cross 
section.  This  leads  to  gross  measurement  errors.  The  use  of 
multiple  plane-of-light  sensors  to  look  at  blade  tip  corners  was 
proposed  by  Camarata  et.al.,  (38).  Chi  et.al.,  (39)  developed 
software  programs  to  process  the  digital  data  from  the  plane-of- 
light  sensors.  Jones  adapted  (40)  the  multichannel  NSMS  system  to 
measure  the  integral  order  vibrations  (41)  of  compressor  blades. 
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United  Technologies  Research  Center  (UTRC)  assembled  a  test  model 
of  this  system  and  demonstrated  it  for  integral  order  vibration 
measurements  for  AEDC  (42) .  Concurrently,  UDRI  designed  and 
assembled,  for  ISRL,  a  NSMS  using  the  same  technological  basis  but 
implemented  in  innovatively  different  hardware  and  software 
systems . 

4.2  PRINCIPLE  OF  OPERATION  OF  PLANE-OF-LIGHT  PROBES 

The  method  measures  the  time  of  arrival  of  the  rotating  blade 
tips  and  the  differential  arrival  time  for  a  vibrating  blade  as 
compared  to  an  unexcited  blade  is  a  direct  indication  of  the  blade 
tip  deflection. 

To  express  the  idea  in  basic  mathematics,  let  a  blade  tip 
travel  with  constant  velocity,  v(=«R) ,  where  w  is  the  angular 
velocity  and  R  is  the  radius  of  blade  tip.  The  distance,  y, 
traveled  by  blade  tip  in  time,  t,  is 

y  »  Vt  (81) 

Introducing  variations  in  variables  y  and  t,  we  obtain 

Ay  =  V(At) .  (82) 

Obviously  from  this  equation,  the  differential  time 
measurements  (At)  will  yield  changes  in  positions  or  deflections 
(Ay)  of  the  blade  tip. 

The  schematic  diagram  in  Figure  33  explains  the  concept 
application  to  a  rotating  bladed  disk.  Let  a  plane-of-light  sensor 
be  mounted  in  a  stationary  case,  surrounding  a  constant  speed,  four 
bladed  rotor.  The  sensor  receives  evenly  spaced  reflected  light 
pulses  off  the  blade  tips  when  blades  are  in  the  undeformed  state. 
Time  interval  (T^)  is  digitally  measured.  If  the  blades  are 
deflected  to  the  left,  the  light  pulses  occur  at  a  later  time  (Tj^) 
and  Tj^  >  Tq.  Similarly,  for  the  blades  deflected  right,  T^  <  T^. 
The  time  differential  in  either  case,  T^^  -  T^  or  -  T^  is  the 
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direct  measurement  of  the  deflection  of  blade  tips  of  a  rotating 
bladed  disk. 

The  above  simple  explanation  can  easily  be  extended  to 
determine  the  average  blade  tip  deflection  (6)  over  several 
rotations  of  the  disk. 

Differentiation  of  equation  (81)  gives 


Ay  =  Avt  +  vAt 


(83) 


Neglecting  small  speed  changes  and  averaging  the  differential 
deflections  and  time  measurements  over  several  rotations  at  a  non¬ 
vibrating  reference  speed  gives 
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Similarly  the  average  deflections  at  a  resonating  speed  are 
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The  net  blade  tip  displacement  over  N  disk  rotations  is 


(85) 


6  =  ^2  -  (86) 

Our  analysis  software  considers  and  corrects  for  shaft  speed 
variations  (Av)  in  equation  (83)  for  each  of  the  recorded 
rotations. 


4.3  PLANE-OF-LIGHT  INSTRUMENTATION  OPTICS  AND  ELECTRONICS 

4.3.1  One  Pulse-Per-Revolution  Electronics  Probe 

The  Applied  Physics  Division  of  UDRI  designed, 
constructed,  and  implemented  an  electro-optic  system  to  provide 
accurate  position  information  about  the  angular  orientation  of  a 
turbine  disk.  The  system  provided  a  fast-rise  time,  one  pulse-per 
revolution  (1  PPR)  signal  to  provide  synchronous  clocking  pulses. 

99 


The  1  PPR  encoder  provided  pulse  rise  time  of  0.5  microsecond 
results  in  an  angular  disk  rotation  of  0.003  degree  for  the 
rotational  speed  of  1,000  RPM  and  0.018  degree  for  a  rotational 
speed  of  6,000  RPM.  These  angles  represent  0.0003  inch  and  0.0018 
inch  of  blade  tip  motion  for  a  12-inch-diameter  fan.  The  signal 
was  processed  to  provide  up  to  1  part  in  100  resolution  during  the 
pulse  rise  and  the  device,  therefore,  provided  up  to  0.00003  degree 
of  angular  resolution  at  1,000  RPM  and  0.00018  degree  of  resolution 
at  6,000  RPM. 

The  UDRI  1  PPR  encoder  relied  on  a  long  lever  arm  to 
sweep  the  retro-reflected  laser  beam  across  a  detector  surface  at  a 
high  velocity.  This  was  accomplished  by  a  20x  microscope  and  a  3- 
meter  lever  arm.  The  optical  arrangement,  including  the  top  and 
side  views  of  the  1  PPR  encoder  mounted  to  the  chamber  with  the 
test  disk  is  shown  in  Figure  34. 

The  single  pulse  per  revolution  encoder  contained  three 
modules:  (1)  a  Helium-Neon  laser  light  source,  (2)  microscope 

assembly,  and  (3)  a  PIN  photodiode  detector.  The  associated  optics 
consisted  of  one  beam  splitter,  an  access  window  to  the  chamber, 
four  aluminum  flat  mirrors,  and  mounting  hardware.  The  mirrors  and 
beam  splitter  directed  the  laser  beam  to  the  microscope  module  and 
back  to  the  detector  module. 

The  Helium-Neon  laser  was  a  CJ  Laser,  Series  1,000, 
with  a  2 -milliwatt  linearly  polarized  output.  The  laser  operated 
in  the  TEM^q  mode  with  a  beam  divergence  of  1.1  milliradians  and 
6 3 2. 8 -nanometer  wavelength.  The  laser  was  self-contained  with 
power  supply,  laser  tube,  and  shutter  housed  in  a  3-  x  12-  x  2 -inch 
housing.  The  laser  was  mounted  on  a  1/2-inch-thick  aluminum  plate 
bolted  to  the  chamber,  and  provided  a  very  stable,  rugged,  light 
source . 

The  microscope  assembly  constructed  by  UDRI,  consisted 
of  a  2 Ox  microscope  objective  mounted  on  a  linear  translation 
stage,  and  an  aluminum  turning  mirror  mounted  on  a  three-degree-of- 
motion  mount.  This  provided  motion  in  the  vertical  direction  as 
well  as  pitch  and  roll  motion  to  allow  the  He-Ne  beam  to  travel 
through  the  optical  axis  of  the  microscope  objective.  The  linear 
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translation  of  the  objective  allowed  fine  focusing  of  the 
microscope.  The  assembly  was  mounted  on  a  1.5-inch-thick  aluminum 
block  which  was  clamped  to  the  aluminum  table  on  which  the  test  fan 
drive  motor  was  mounted. 

The  detector  was  a  Hewlett-Packard  PIN  photodiode 
detector.  The  signal  was  conditioned  by  UDRI-designed  electronics. 
The  electronic  circuit,  shown  in  Figure  35,  provided  a  TTL 
(transistor-transistor  logic)  signal  compatible  with  the  data 
acquisition  system. 

The  He-Ne  beam  path,  shown  in  Figure  36,  was  as 
follows:  The  He-Ne  laser  beam  was  transmitted  through  a  beam 

splitter.  The  reflection  portion  was  collected  by  a  beam  dump  and 
the  transmitted  portion  was  directed  into  the  chamber,  by  an 
aluminum  mirror,  through  a  1-inch-diameter  glass  window  in  the 
chamber  cover.  As  the  beam  reached  the  microscope  module,  a 
turning  mirror  directed  the  beam  through  the  microscope  objective. 
This  focused  the  light  onto  a  mirror  on  the  drive  motor  shaft  which 
reflected  the  beam  back  through  the  same  optics.  The  reflected 
beam  returned  to  the  beam  splitter  which  directed  a  portion  of  it 
to  the  detector.  The  horizontal  turning  mirror  on  top  of  the 
chamber  lengthens  the  optical  path  causing  the  beam  to  sweep 
through  a  longer  arc,  and  this  long  optical  path  shortened  signal 
rise  time  for  the  detector. 

Pulse  shape  with  the  optical  1  PPR  system  was  shown  in 
Figure  37.  A  2  microsecond  actual  rise-time  for  the  signal  shown 
was  achieved  which  proved  sufficient  to  ensure  accuracy  of  half-a- 
mil  in  blade  tip  deflection  measurements. 

4.3.2  Plane  of  Light  Sensor  Development 

The  UDRI  Optics  Group  designed  a  fiber-optic  plane  of 
light  sensor  that  provided  fast  rise-time  and  easy  operation  for 
the  turbine  blade  deflection  measurements.  The  sensor  geometry  and 
the  emitting  and  sensing  end  details  are  shown  in  Figures  38  and 
39  respectively.  The  linearly  arrayed,  bifurcated  sensor  provided 
light  input  and  collection  within  a  single  assembly.  Input  and 
output  bundles  were  separate  and  the  fibers  had  minimal  cross-talk. 
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(See  Below) 


Figure  38.  Plane-of-Light  Sensor  Geometry. 


About  75  individual  fibers  were  in  both  the  input  and  output 
portions  of  the  sensor.  The  active  areas  of  illumination  and 
collection  were  each  0.5  x  0.002  inch  with  a  separation  of  0.002 
inch  between  arrays.  This  configuration  provided  a  compact  source 
of  illumination  and  plane  of  visual  acceptance. 

Some  details  of  the  fibers  and  the  probe  assembly  were 
provided  below. 

Fiber:  200  microns  core  diameter 

230  microns  ±  5  microns  clad  diameter 

Line  dimensions:  0.5*'  x  0.01”  nominal  (both) 

Note  #1:  The  thickness  of  the  line  will  be  de¬ 
pendent  on  the  diameter  of  the  fiber. 

Note  #2:  Lines  to  be  separated  by  a  shim 

(0.002”  to  0.005”) 

Aluminum  Housing  Dimensions:  2”  x  0.75”  x  0.2” 

Sheathing:  PVC/Monocoil 

Overall  Cable  Length:  3  1/2' 

The  sensor  made  to  the  specifications  was  supplied  by 
the  Applied  Fiber  Optics,  Southridge  MA.  All  probes  had  well 
polished  ends  and  good  acceptance  angles. 

We  used  the  sensors  to  detect  the  passing  of  the  tips 
of  simulated  turbine  fan  blades.  When  the  blade  passed  the  fiber 
optic  linear  array,  each  individual  fiber  collected  the  reflected 
light  and  signified  the  actual  position  of  the  blade.  The  line 
sensor  output  was  coupled  to  a  linear  detector  array  conveniently. 

The  light  from  each  of  the  output  fibers  was  integrated 
to  result  in  a  single  output  pulse  signifying  the  passage  of  each 
blade.  The  two  approaches  used  to  accomplish  this  integration 
were  (1)  collecting  the  output  light  with  a  diffusing  plate  after 
the  fiber  optic,  and  (2)  collecting  the  output  light  with  an 
integrating  sphere.  The  overall  pulse  width  was  increased,  and  the 
pulse  shape  was  other  than  a  pure  square  wave,  but  the  front  edge 
rise  time  was  very  short.  The  electronic  thresholding  was  adjusted 
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to  the  sharp  transition  that  corresponded  to  the  first  Individual 
optic  fiber  sensing  in  the  bundle. 

4.3.3  NSMS  Instrumentation  Electronics 

A  software  program  was  developed  and  implemented  to 
interface  a  nine  track  tape  drive  with  the  microprocessor  unit 
(MPU) .  The  MPU  and  parallel  interface  were  tested  successfully  to 
verify  reading  and  writing  of  the  various  ports. 

Three  input/ output  (I/O)  cards  and  three  MPU  were 
acguired,  assembled  and  tested  in  parallel  with  the  three  fiber 
optic  plane-of-llght  probes.  An  additional  I/O  card  and  MPU  were 
designed  and  integrated  with  probe  data  acquisition  system  to  work 
as  system  controller.  The  menu  for  the  entire  data  acquisition, 
display  and  write  was  driven  from  this  controller. 

The  prototype  counter  and  latch  cards  were  fabricated 
and  tested  for  arbitrary  short  cable  run  lengths,  but  later 
redesigned  to  allow  the  maximum  counter  reading  to  be  latched.  The 
latch  card  was  redesigned  using  differential  line  drivers  to  drive 
the  longer  cables.  Because  the  drivers  were  differential,  a 
differential  line  receiver  card  was  designed  to  interface  to  the 
I/O  card.  These  cards  were  fabricated  and  tested  for  satisfactory 
functioning.  Because  of  the  additional  circuitry,  a  larger  power 
supply  for  the  extra  load  was  required. 

4.3.4  Plane-of-Liqht  Sensor  System  Operation 

A  brief  description  of  the  interaction  between  four 
optical  probes  and  the  process  electronics  was  provided  here.  In 
the  systems  data  acquisition  mode,  each  processor  was  designed  to 
take  a  time  reading  as  each  blade  passed  the  sensor.  The  source 
illumination  through  the  fiber  optic  bundle  reached  the  turbine 
blade  tip  and  the  reflected  light  reached  the  photomult ipler  tube 
(PMT) .  The  PMT  output  signal  proportional  to  the  light  level  was 
amplified  and  compared  with  the  threshold.  If  the  light  level  was 
above  the  threshold,  then  the  time  of  arrival  was  recorded  by  the 
processor.  The  electronic  block  diagram  shown  in  Figure  40 
explains  the  system  operation. 
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Figure  40.  Data  Acquisition  Operation  by  the  Optical 
Probe  System. 


The  time  for  the  four  sensors  was  derived  from  a  20-MHz 
crystal  oscillator.  This  clock  time  was  divided  by  a  2  to  generate 
a  buffered  complementary  10  MHz,  50-percent  duty  cycle  square  wave. 
The  clock  output  was  used  to  clock  a  20-bit  synchronous  counter 
from  five  4 -bit  counters  with  external  look  ahead  carry  to  speed  up 
their  operation.  The  outputs  from  each  counter  card  were  buffered 
to  drive  a  latch  card.  The  counter  was  reset  by  the  one  pulse  per 
revolution  encoder  on  the  turbine  blade  shaft.  The  encoder  was 
asynchronous  with  the  master  clock  but  the  counter  reset  must  by 
synchronized  with  the  clock  which  was  accomplished  by  using  three  D 
flip-flops  (bistable  toggle  switches} .  In  these,  preset  and  clear 
were  overriding  controls.  The  first  was  set  high  by  the  rising 
edge  of  the  reset  pulse  from  the  encoder.  This  enables  the  second 
D  flip-flop  which  was  set  by  the  next  leading  edge  of  the  10-MHz 
clock.  The  counter  was  reset  by  this  signal  which  also  enables  the 
third  D  flip-flop  which  was  set  by  the  following  clock  pulse 
leading  edge.  The  output  of  the  third  D  flip-flop  cleared  the 
first  and  second  D  flip-flops.  When  the  D  flip-flop  was  cleared, 
the  next  leading  edge  of  the  clock  clears  the  third  D  flip-flop. 

When  the  light  level  on  a  PMT  exceeds  the  threshold 
value,  the  current  counter  value  was  set  in  the  latch  card  for  that 
channel.  The  PMT  output  was  asynchronous  with  the  10-MHz  clock  and 
was  synchronized  with  the  clock  by  the  same  type  of  circuit  as  that 
used  to  synchronize  the  reset  for  the  counter.  The  latch  signal 
was  also  used  to  trigger  the  computer  to  read  the  latched  values. 
The  latches  had  a  high  current  output  stage  to  drive  the  cables 
connecting  the  front  end  electronics  to  the  computers. 

Each  data  channel  computer  recorded  the  arrival  time  of 
each  blade  as  it  reached  the  sensor  of  that  channel.  The  four  data 
channel  computers  record  the  four  stations.  After  the  data  had 
acquired,  the  data  channel  computers  transferred  the  data  to  the 
system  computer  for  recording  on  a  nine-track  magnetic  tape.  There 
were  separate  records  from  the  four  data  channel  computers.  Each 
record  might  start  with  a  different  blade  number  occurring  as  the 
first  blade  in  that  record,  depending  on  the  blade  and  sensor 


spacing.  The  system  computer  also  controlled  the  setup  and 
collection  of  data  for  the  four  channels. 

The  design  of  the  four-channel  system  Including  counter  cards 
and  latch  cards  was  described  In  full  In  References  (1,2). 

4.4  TURBINE  BLADE  DATA  ACQUISITION  SYSTEM 

The  UDRI  designed  turbine  blade  data  acquisition  system 
consisted  of  five  major  components  -  the  1  PPR  detector,  1  PPR 
amplifier,  the  detector/electronic  chassis,  the  computer  system 
rack  and  the  video  graphics  terminal.  A  printer  was  attached  to 
the  system  for  obtaining  hard  copy  outputs  of  the  terminal  screen. 
These  components  were  shown  In  the  block  diagram.  Figure  41. 

The  1  PPR  detector  together  with  Its  amplifier  provided  a  TTL 
compatible  signal  to  synchronize  the  system  operation.  The 
detector/electronics  chassis  received  the  signals  from  the  1  PPR 
and  four  probes  and  converted  them  to  time  data  which  were 
transferred  to  the  computer  for  storage.  A  10-MHz  counter  In  the 
detector/electronics  chassis  drove  five  latch  cards,  one  for  the  1 
PPR  and  four  for  four  probes.  When  the  system  received  a  1  PPR 
pulse,  the  value  of  the  counter  was  latched  In  to  the  1  PPR  latch 
card,  and  then  the  counter  was  reset.  When  a  blade  passed  a  probe, 
a  pulse  was  generated  and  the  counter  value  was  latched  Into  the 
card  corresponding  to  that  station.  The  card  then  sent  a  valid 
data  signal  to  the  computer  which  read  the  data  and  reset  the  card. 
Thus  all  data  were  referenced  In  time  to  the  1  PPR  and  the  blade 
arrival  time  data  extracted. 

A  nine-track  tape  unit  was  provided  to  allow  the  data  to  be 
stored.  This  was  a  wrlte-only  tape  system,  and  the  data  was  to  be 
read  and  analyzed  on  another  computer. 

The  menu  driven  operating  checklist  provided  for  Item  entries 
such  as  test  date  and  time,  test  specimen  details,  data 
acquisition,  data  display  and  data  write  and  data  discard.  The 
system  presents  two  modes  of  display  -  one,  histogram  for  each  of 
the  four  probe  data  versus  number  of  blades,  and  two,  a  real-time 
display  of  blade  arrival  times  for  up  to  five  blades  at  a 
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particular  station.  The  user  can  sequentially  obtain  and  display 
data  for  all  the  probes. 

A  screen  display  of  control  monitor  menu  scheme  was  shown  In 
Figure  42.  Data  values  (number  of  counts)  versus  the  number  of 
blades  of  the  rotating  test  disk  were  displayed  In  a  typical 
histogram  In  Figure  43.  This  data  display  could  be  updated  for  a 
pair  of  revolutions  of  user's  choice.  The  gradual  step  Increase  In 
data  values  was  due  to  speed  change  of  the  test  disk  In  each 
revolution.  Figure  44  shows  a  typical  real-time  display  of  the 
data  for  five  blades  from  a  particular  probe.  The  data  values  were 
normalized  with  reference  to  the  1  PPR  count  values.  Similar  data 
displays  were  provided  by  each  of  the  three  plane-of-llght  probes. 
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NSMS  Data  Acquisition  System 


Revolution  204  of  1000 

Current  Minimum  Maximum 
Speed:  625  625  625  RPM 

C  =  Change  Date  /  Time  /  Specimen  Data 
D  =  Discard  Data 
H  =  Histogram 
R  =  Real  Time  Display 
S  =  Start  Taking  Data 
W  =  Write  Data  To  Tape 


Figure  42,  Menu  Display  of  the  NSMS  System 
Operation. 
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Figure  44.  Real  Time  Display  of  Blade  Deflection  from  NSMS 


SECTION  V 

Z-GRID  INSTRUMENTATION  DEVELOPMENT 


In  view  of  the  objective  of  this  technical  program,  we  also 
developed  a  nonintrusive,  frequency  modulated  (FM)  vibration 
measurement  system  and  evaluated  its  performance  at  ISRL.  The 
basic  concept  involved  in  this  technique  was  that  a  miniature 
permanent  magnet  mounted  at  the  tip  of  a  rotating  and  vibrating 
turbine  or  compressor  blade  would  modulate  a  carrier  frequency  of  a 
case  mounted  stationary  grid,  and  the  subsequent  demodulation  and 
integration  provided  a  signal  which  was  proportional  to  the 
vibratory  motion  of  the  blade.  The  proposed  technique  had  a  good 
potential  for  application  in  real  engine  systems  in  operating 
conditions.  Therefore,  UDRT  pursued  this  technique  and  designed  an 
FM  technique  to  demonstrate  the  process  under  laboratory 
conditions.  The  following  paragraphs  describe  various  phases  of 
grid  wire  development  and  its  application  for  measuring  blade 
vibrations. 

The  first  phase  involved  the  tasks  of  selection  of  a  suitable 
permanent  magnet,  the  design,  fabrication  and  evaluation  of  a 
single  channel  FM-grid  to  generate  acceptable  signal  level.  The 
next  phase  involved  the  design  and  testing  of  a  two-channel  FM 
system  for  measuring  two  components  of  vibratory  motion  of  an 
instrumented  blade  and  comparing  the  results  with  those  of  an  array 
of  fiber  optic  plane-of-light  sensors. 

5.1  PRINCIPLE  OF  OPERATION  OF  THE  GRIDWIRE  SYSTEM 

The  basic  principle  of  operation  of  the  proposed  frequency 
modulated  system  and  its  implementation  to  measure  the  vibratory 
motion  of  a  blade  in  the  laboratory  conditions  are  explained  here. 

A  conductive  wire  grid  was  firmly  embedded  in  the  casing  wall,  and 
the  grid  was  laid  out  in  the  form  of  an  evenly  pitched  square  wave. 
A  permanent  magnet  such  as  Samarium  Cobalt  buried  in  the  tip  of  a 
rotating  blade  generated  an  electromagnetic  force  (emf)  in  the 
stationary  grid  circuit  each  time  it  crossed  a  bar  of  the  grid. 
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Figure  46.  Two  Square-Grid  Wire  Layout. 
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Several  runs  were  made  for  each  grid,  with  different  air  gap 
separations  between  the  grid  and  the  magnet.  In  all  these 
experiments,  the  aluminum  disk  carrying  the  magnet  was  positioned 
parallel  to  the  stationary  wooden  disk  with  grids.  The  magnet 
swept  past  the  grid  wires  with  constant  velocity.  The  disk  was  run 
at  several  speeds  and  the  output  terminals  were  directly  connected 
to  an  oscilloscope.  Amplitude  and  frequency  of  the  resulting 
signals  were  read  off  the  oscilloscope  screen.  Results  are 
presented  in  Tables  20  through  23  and  also  in  graphical  form  in 
Figures  47  and  48.  These  graphs  show  the  dependence  of  the  output 
signals  on  the  variables  of  speed  of  rotation  (velocity)  and  gap 
between  the  magnet  and  grid. 

The  observed  output  signals  were  clean  at  all  speeds.  Flats 
between  pulses  were  more  obvious  for  the  Grid  1  whose 
circumferential  spacing  was  large  (0.51  inch)  compared  to  the 
magnet  size  (0.25  inch).  Figure  49  shows  a  gradual  increase  in 
quality  output  signal  level  with  the  increase  in  rotational  speed. 
Output  signal  levels  of  2.5  mV  when  the  magnet  was  separated  by  a 
distance  of  10  mm  from  the  grid  wires  was  considered  adequate  for 
further  signal  processing. 

We  have  also  studied  the  grid  emf  signals  from  a  partial  grid 
loop  installed  in  a  12.4-inch-inner-diameter  instrumentation  ring. 
Two  1/4  X  1/8 -inch-square  super  magnets  were  installed  on  opposite 
blades  of  the  0.065-inch-thick,  38-bladed,  11-inch-diameter  test 
disk  which  was  installed  in  the  spin  test  setup.  We  mounted  the 
12.4-inch  ID  grid  support  ring  eccentrically  to  the  magnet  spin 
circle  with  various  clearances  to  the  12  grid  wires  as  shown  in 
Figure  50.  The  test  disk  then  was  spun  at  speeds  up  to  4,200  rpm. 
The  voltage  signals  generated  in  the  test  grid  at  speeds  of  3,600 
and  4,200  rpm  for  the  passage  of  the  strongest  of  the  two  magnets 
past  the  grid  are  shown  in  Figure  51.  The  minimum  pulse  amplitudes 
of  35  to  45  millivolts  and  the  pulse  shapes  certainly  were  amenable 
to  recording  of  the  pulse  peak  occurrence  times  on  the  NSMS 
instrumentation  system. 
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TABLE  20 

GRID  WIRE  SIGNAL  OUTPUT  FROM  GRID  1  WITH  GAP  1 


Grid  Details:  Inner  Radius  *  4.5”;  Outer  Radius  «  5.0” 

Number  of  Grid  Wires  =  58;  Gage  No.  26 

Magnet  Details:  Dia.  =  1/4”;  Length  -  1/8”;  Path  Radius  =  4.75” 


Gap  between  magnet  tip 
and  grid 


Smm 


Speed  (rpm) 

Signal  Output  (mV) 

Comments 

200 

10.0 

Triangular  pulses  with 
flats  in  between 

400 

15.0 

Triangular  pulses  with 
flats  in  between 

600 

17.5 

Triangular  pulses  with 
flats  in  between 

800 

20.0 

Triangular  pulses  with 
flats  in  between 

1000 

25.0 

Triangular  pulses  with 
flats  in  between 

1200 

26.0 

Triangular  pulses  with 
flats  in  between 

1400 

30.0 

Beating  of  signals 

1600 

35.0 

* 

2000 

40.0 

Noisy 

2400 

47.5 

- 

3000 

50.0 

Beating  of  Signals 
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TABLE  21 

GRID  WIRE  SIGNAL  OUTPUT  FROM  GRID  1  WITH  GAP  2 

Grid  Details:  Inner  Radius  =  4.5";  Outer  Radius  «  5.0" 

Number  of  Grid  Wires  =  58;  Gage  No.  26 

Magnet  Details:  Dla.  =  1/4"; 


Gap  between  magnet  tip  3^ 

and  grid  Left 

Window 


8  nan 


Speed  (rpm) 

Signal  Output  (mV) 

Comments 

200 

7.5 

Triangular  pulses  with 
flats  In  between 

400 

10.0 

Triangular  pulses  with 
flats  in  between 

800 

12.0 

Triangular  pulses  with 
flats  in  between 

1200 

15.0 

Beating  of  Signals 

1600 

17.0 

- 

2000 

20.0 

- 
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TABLE  22 

SIGNAL  OUTPUT  FROM  GRID  1  WITH  GAP  3 


Grid  Details:  Inner  Radius  ==  4.5”;  Outer  Radius  =  5.0” 

Number  of  Grid  Wires  »  58;  Gage  No.  26 


Speed  (rpm)  Signal  Output  (mV)  Comments 


200 

2.5 

Less  Triangularity  in 
pulses  and  with  less 
flats  between  them. 

400 

7.5 

Less  Triangularity  in 
pulses  and  with  less 
flats  between  them. 

800 

10.0 

Less  Triangularity  in 
pulses  and  with  less 
flats  between  them. 

1200 

12.0 

Less  Triangularity  in 
pulses  and  with  less 
flats  between  them. 

1600 

15.0 

Less  Triangularity  in 
pulses  and  with  less 
flats  between  them. 

2000 


16.0 


Less  Triangularity  in 
pulses  and  with  less 
flats  between  them. 


Grid  Details 


TABLE  23 

SIGNAL  OUTPUT  FROM  GRID  2 


Inner  Radius  *  3.5”;  Outer  Radius  =  4.0” 
Number  of  Grid  Wires  =  58;  Gage  No.  26 

Magnet  Details:  Dia. 


Gap  between  magnet  tip 
and  grid 


Speed  (rpm) 

Signal  Output  (mV) 

Comments 

200 

7.5 

Triangular  pulses  with 
less  flats  in  between 

400 

8.0 

Triangular  pulses  with 
less  flats  in  between 

800 

10.0 

Triangular  pulses  with 
less  flats  in  between 

1200 

15.0 

Triangular  pulses  with 
less  flats  in  between 

1600 

17.0 

Triangular  pulses  with 
less  flats  in  between 

2000 

20.0 

Triangular  pulses  with 
less  flats  in  between 
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Figure  48.  Grid  Wire  Signal  Growth  with  Speed  and  Gap. 
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800  RPM 


Figure  49.  Growth  of  Output  Signal  with  Speed. 


129 


12.4"  ID  Grid  Support 
Ring  Mounted  Eccentrically 
To  Axis  of  Rotation  of 
11.0"  OD  Bladed  Disk 


Two  Rare  Earth  Supermagnets 
1/8"  Squqre  x  1/4"  Long 
Bonded  to  Blades  10  and  29 
To  Extend  1/10"  from  Blade 
Tips 


Figure  50.  Grid  Wire  with  Varying  Magnet  Gaps. 


?Omy  >1"'^ 


?CCmV 


3,600  rpm 


4,200  rpm 


Figure  51.  Oscilloscope  Display  of  Eccentric  Grid 
Wire  Signals. 
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A  ring  with  36  grid  wires  described  in  the  next  section  was 
supported  symmetrically  around  the  ISRL  test  shaft  spin  axis  by 
bolting  it  to  the  optical  probe  supporting.  A  source  magnet  and 
counter  weight  were  mounted  to  the  tips  of  two  diametically 
opposite  blades  on  an  11-inch-diameter  disk  as  shown  in  Figure  52. 
The  test  disk  was  installed  on  the  drive  shaft  in  axial  alignment 
with  the  center  of  the  grid  system  ring.  The  two  grid  system 
signal  leads  then  wired  to  a  dual  channel  oscilloscope. 

The  resulting  pulse  trains  from  the  two  grids  as  displayed  on 
the  oscilloscope  screen  were  shown  in  Figure  53.  The  regularity  of 
the  pulse  amplitude  and  pulse  intervals  for  a  spin  speed  of  4,noo 
rpm  with  a  magnet -grid  gap  of  0.060  inch  was  extremely  good. 

Signal  amplitude  was  adequate,  and  there  was  no  mutual  interference 
contamination  of  the  two  grid  signals.  The  pulse  arrival  times 
were  resolved  by  the  ISRL  NSMS  since  the  signal  trace  widths  were 
due  to  trace  blooming  in  the  screen  phosphors,  as  expected,  and  not 
due  to  high  frequency  noise. 

5.3  DESIGN  OF  THE  DUAL  Z-GRID  SYSTEM 

5.3.1  Z-Grid  Support  Rina 

As  a  result  of  preliminary  investigation  on  grid 
signals,  a  grid  support  ring  was  designed  for  use  with  an  11-inch- 
diameter  disk.  Two  orthogonally  crossed  grid  wires  (30  AWG)  at  45° 
to  the  rotational  axis  of  the  disk  were  used  to  resolve  both 
horizontal  and  axial  components  of  an  instrumented  blade  vibration. 
These  dual  grid  wires  overlaid  in  the  letter  form  Z  (hence,  Z-Grid 
system)  were  insulated  and  wound  around  the  inner  diametral  surface 
of  a  support  ring.  The  geometry  and  the  layout  of  the  dual  grid 
wire  system  and  its  support  ring  are  shown  in  Figure  54.  The  grid 
support  ring  with  thirty-six  30-mil-deep  grooves  to  accommodate  the 
conducting  grid  wires  fabricated  from  polypropylene  material  was 
coaxially  installed  with  the  test  disk.  Thirty-six  grids  for  each 
of  the  two  grid  loops  were  used  to  resolve  at  least  nine  modes  of  a 
rotating  disk  with  sufficient  accuracy. 
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Figure  52.  Source  Magnet  on  One  of  38  Blades  on  Concentric 
Disk  with  Wire  Grid. 


Figure  53.  EMF  Pulse  Train  for  Concentric  Wire  Grids. 
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5.3.2  Digital  Design  of  the  Grid  Wire  System 

The  optical  pulse  train  due  to  the  reflections  off  the 
blade  tips  in  the  optical  system  and  the  emf  bipolar  pulse  train 
due  to  a  source  magnet  in  the  Z-grid  system  had  similar  digital 
characteristics.  Both  signals  after  referring  to  the  1  PPR  signal 
will  yield  differential  time  measurements  for  a  specific  blade. 
This  differential  time  measurement  was  a  direct  measure  of  blade 
tip  deflection  in  the  optical  and  Z-Grid  systems.  This 
measurement-similarity  in  both  systems  prompted  us  to  design  a 
digital  Z-Grid  system  for  measuring  blade  tip  deflections  of  a 
rotating  disk.  Instead  of  analog  demodulation  of  the  grid  signals 
as  was  done  by  the  Rolls  Royce  Company,  UDRI  designed  a  digital 
system  for  the  grid  signal  demodulation.  This  step  allowed  us  to 
use  the  existing  optical  (NSMS)  data  acquisition  and  data  storage 
hardware.  UDRI  developed  speed  compensation  algorithm,  digital 
time-series,  and  Fast  Fourier  Transform  analysis  also  were  used 
effectively  with  minimum  cost  and  time. 

However,  the  following  modifications  had  to  be 
incorporated  into  the  Z-Grid  system  electronics  to  make  it 
compatible  to  the  optical  system.  First,  bipolar  grid  signals  had 
been  converted  to  monopolar  pulse  train.  These  signals  were  then 
amplified  and  filtered  by  a  bank  of  filters,  by  an  electronic 
signal  conditioning  arrangement  as  shown  in  Figure  55. 

5.4  Z-GRID  SYSTEM  ELECTRONICS 


The  Z-Grid  signal  processing  system  was  designed  to  process 
two  frequency  modulated  waveforms.  These  waveforms  were  generated 
by  the  motion  of  a  magnet  mounted  on  one  rotating  turbine  blade, 
past  two  wire  grid  rings.  The  individual  conductors  of  the  two 
grids  were  insulated  from  each  other  and  crossed  orthogonally.  The 
rotation  of  the  blade  produced  variable  phase  differences  between 
the  two  FM  signals,  dependent  on  the  location  of  the  magnet  as  it 
passes  under  the  series  of  paired  conductors  of  the  grid.  The 
signal  processing  system  processed  these  waveforms  and  generated 
pulses  corresponding  to  the  peaks  of  the  input  waveforms.  These 
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pulses  drive  a  light  emitting  diode  (LED)  for  each  of  the  two 
channels.  The  LEDs  replace  two  of  the  fiber-optic  bundles  on  the 
NSMS  Data  Acquisition  System.  The  data  were  then  acquired  similar 
to  that  of  plane-of-light  sensor  method. 

5.4.1  System  Connections 

The  Z-Grid  signal  processing  system  hardware  consisted 
of  two  combination  power  supply  -  signal  conditioning  units,  one 
for  each  of  the  two  grid  wire  circuits.  After  mounting  these  units 
in  the  vibration-instrumentation  rack  of  the  ISRL  rotating  test 
rig,  two  grid  wire  circuits  and  the  digital  hardware  were 
connected,  because  this  system  was  designed  to  be  used  in 
conjunction  with  the  NSMS  data  acquisition  system.  A  BNC  cable 
from  the  grids  to  the  preamplifiers  for  the  signal  conditioning 
system  was  connected.  The  data  acquisition  block  diagram  for  the 
digital  Z-Grid  system  is  illustrated  in  Figure  56. 

5.4.2  System  Operation 

The  operation  of  the  Z-Grid  signal  processing  system 
was  designed  to  be  transparent  to  the  user  of  the  NSMS  System.  All 
procedures  for  operating  the  system  were  identical  to  those  of  NSMS 
except  setting  the  nximber  of  blades  equal  to  36,  where  the  number 
36  was  equal  to  the  number  of  grid  wires  in  the  dual  grid  rings. 

The  output  display  for  both  channels  can  be  accessed  in  the  same 
way  as  in  NSMS.  A  typical  histogram  for  one  of  the  Z-Grid  channels 
is  shown  in  Figure  57.  This  histogram  resembled  the  optical  system 
histogram  except  for  the  number  of  blades  on  the  abscissa. 

Two  adjustments  on  the  Z-Grid  signal  processing  system 
were  the  preamp  level  and  the  time  delay  adjustment.  These  were 
adjusted  when  a  new  speed  was  used  or  the  configuration  of  the 
system  was  changed. 

The  preamp  level  was  adjusted  while  the  system  was 
running  for  an  output  level  between  one-half  and  one  volt,  and  this 
was  achieved  by  using  a  specific  potentiometer  on  the  preamp  board. 
At  present  the  Z-Grld  signal  conditioning  system  was  set  to  process 
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Figure  56.  Data  Acquisition  Block  Diagram  for  the 
Digital  Z-Grid  System. 


station 


Figure  57,  Z-Grid  Histogram  Display 


data  for  the  speed  range  of  1,000-6,000  rpm.  It  was  possible, 
however,  to  change  these  settings  for  other  test  speeds. 

The  first  method  of  setting  the  delay  was  to  set  the 
switches  based  on  a  given  speed.  Using  this  method,  we  located  the 
speed  to  run  the  turbine  in  the  RPM  column  and  set  the  dip  switches 
on  the  output  card  using  the  binary  switch  column  of  the  chart. 

The  second  method  was  to  set  the  switches  for  a  given  delay  time. 
The  relationship  between  the  grid  input  and  LED  output  signals  is 
shown  in  Figure  58.  The  delay  time  that  triggers  the  generation  of 
LED  input  pulses  of  the  grid  wire  emf  signal  peak  had  to  be  set  so 
that  the  output  pulse  fell  within  the  period  labeled  AT.  We  needed 
to  ensure  that  the  delay  time  was  not  so  short  that  the  output 
pulse  occurred  within  the  same  input  waveform,  nor  so  long  that  the 
output  pulse  occurred  during  the  next  input  waveform. 

The  Z-Grid  system  design,  electronics,  operation,  and 
diagnostic  checks  were  fully  described  in  Reference  (4) . 
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Figure  58 .  Grid  Wire  EMF  and  LED  Output  Signal  Relationship. 


SECTION  VI 

NSMS  BLADE  DEFLECTION  DATA  REDUCTION  AND  ANALYSIS 

UDRI  developed  a  set  of  data  analysis  programs  to  extract  the 
deflection  and  frec[uency  Information  from  the  data  generated  by  the 
two  noncontacting  blade  deflection  measurement  systems  discussed 
earlier.  The  four-channel  plane-of-llght  system  provided  vibratory 
data  for  all  blades  of  a  rotating  bladed  disk  at  a  specific  sensor 
location  while  the  dual  Z-Grld  system  provided  similar  dynamic  data 
for  one  blade  over  the  entire  revolution.  Both  data  sets  are 
processed  by  the  UDRI  algorithms.  Both  systems  acquire  the 
vibration  data  through  a  common  data  acquisition  system  and  stores 
It  on  a  nine-track  magnetic  tape.  The  NSMS  data  acquisition  system 
Is  not  capable  of  analyzing  the  data  In  real-time;  therefore,  the 
data  analysis  had  to  be  performed  In  a  separate  computer  system.  A 
set  of  UDRI  developed  analysis  programs  performed  the  data 
reduction  and  analysis  tasks  on  a  DEC  VAX  computer  at  the  ASD 
Computer  Center.  The  following  subsections  describe  raw  data 
acquisition,  data  correction  for  speed  variation,  dynamic  analysis 
by  optical  and  Z-Grid  methods,  and  the  FPT  spectra  determination  of 
blade  tip  deflection. 

6.1  DATA  REDUCTION  AND  ANALYSIS  STEPS 

The  rotating  blade  deflection  and  frequency  derivation  cycle 
consists  of  four  major  steps  (Figure  59) : 

1)  Make  data  acquisition  runs  in  the  laboratory  and  store  the 
data  on  a  magnetic  tape. 

2)  Store  the  data  tape  In  the  ASD  Computer  Center  tape 
library  and  transfer  these  data  on  command  to  ASD  VAX 
disk. 

3)  Run  the  appropriate  software  programs  using  the  test  tape 
data  as  Input  and  plot  blade  deflections  In  time  and 
frequency  domains. 

4)  Archive  the  data  onto  the  ASD  central  file  system. 
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6.2.  TEST  SHAFT  SPEED  CORRECTION  EFFORTS 


Early  data  runs  showed  Instability  in  the  speed  of  the  test 
shaft  drive  motor  as  well  as  Intermittent  inconsistencies.  To 
provide  information  on  the  uniformity  of  the  motor  speed  within  a 
single  revolution,  a  30  faceted  polygon  mirror  was  fastened  to  the 
motor  shaft.  Data  were  taken  on  this  polygon  with  the  fiber  optic 
plane-of-light  sensor.  A  He-Ne  laser  served  as  the  light  source 
and  the  fiber  optic  plane-of-light  sensor  was  set  at  approximately 
one  meter  from  the  polygon  to  produce  a  2-microsecond  rise-time. 
Figure  60  shows  the  installation  of  the  polygon  on  the  test  shaft, 
and  Figure  61  shows  the  setup  for  making  the  measurements. 

Sets  of  data  were  taken  on  the  turbine  fan  with  and  without 
vibration  excitation.  All  data  were  taken  at  low  speed  and  with 
atmospheric  pressure  in  the  test  chamber.  The  first  sets  of  data 
exhibited  features  which  were  not  well  understood.  Figure  62  shows 
two  sets  of  data  taken  on  the  12-bladed  fan  spinning  at  1,150  rpm 
with  no  external  excitation.  These  two  sets  of  data  show  the 
repeatability  of  the  plane-of-light  sensor  measurements.  Figure  63 
shows  a  data  run  with  the  fan  spinning  at  1,430  rpm.  These  data 
sets  show  a  pseudo-displacement  measurement  which  has  a  frequency 
of  one-fourth  of  the  fan  speed  at  1,150  rpm  and  one-fifth  of  the 
fan  speed  at  1,430  rpm.  We  believed  these  data  showed  variations 
in  the  test  shaft  speed,  within  each  revolution,  caused  by  the  16 
PPR  shaft  speed  feedback  control  system.  Note  that  amplitudes  were 
very  low  for  blades  1  and  12.  Speed-variation-caused  data  peaks 
were  minimized  on  these  blades  because  they  were  adjacent  to  the  1 
PPR  mirror,  and  all  the  data  readings  were  normalized  by  the  1  PPR 
values.  We  expected  that  if  these  were,  in  fact,  shaft  speed 
changes,  as  they  appeared  to  be,  they  would  be  greatly  reduced  when 
higher  shaft  speed  testing  was  performed. 

The  optical  simulator  was  set  up  to  provide  accurate 
information  for  the  PMT,  electronics,  and  software.  Optical 
waveforms  for  both  the  1  PPR  and  8  blades  per  revolution  were 
created  with  a  drive  circuit  and  LED  light  sources.  The  rise  time 
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Figure  60.  30  Faceted  Polygon  for  Speed  Checkout 


Figure  61.  Speed  Accuracy  Measurement  Setup 
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of  each  optical  signal  wave  form  recorded  over  100  revolutions 
period  was  better  than  1  microsecond.  The  scale  has  been  adjusted 
to  show  +  one  bit  of  sensitivity  which  showed  that  the  PMT 
electronics  and  software  were  all  operating  at  their  peak 
performance. 

From  several  sets  of  data  for  30  revolutions  of  the  polygon  at 
752,  1,129  and  1,694  rpm,  we  observed  that  the  velocity  uniformity 
appeared  to  improve  at  higher  rpm,  and  the  data  from  the  polygon 
segments  nearest  the  1  PPR  mirror  were  again  minimized  by  the 
normalization  of  the  data  by  the  1  PPR  readings. 

6.2.1  Test  Shaft  Speed  Variation  Improvement  Efforts 

After  identifying  speed  changes  in  a  single  revolution, 
we  diagnosed  the  controller  system  operation  from  the  schematics. 

We  made  several  unsuccessful  attempts  to  replace  the  controller 
timer  and  motor  torque  control  input  with  a  programmed  Apple 
Computer.  The  level  of  speed  stability  achieved  was  approximately 
ec[ual  to  that  achieved  by  the  original  control  system.  One  thing 
we  achieved  was  that  we  discovered  that  the  16  pulse  per  rev  (16 
PPR)  shaft  speed  optical  encoder  was  generating  17  pulses  per  rev 
due  to  a  stray  piece  of  tape  across  one  of  the  16  windows  in  the 
encoder  disk.  Removal  of  this  tape  strip  resulted  in  the  system 
returning  to  its  former  mode  of  operation  of  running  within  1  or  2 
rpm  of  the  speed  set  point  for  long  term  average  speed  but  did  not 
solve  the  problem  of  short  term  speed  changes. 

Inability  to  solve  this  speed  control  problem  with  what 
we  considered  to  be  a  much  Improved  control  signal  generator  led  us 
to  evaluate  the  speed  variations  of  a  synchronous  motor.  For  this 
evaluation,  we  excited  the  field  coils  of  a  commutator-less  dc 
motor  having  a  four-pole  permanent  magnet  armature  with  a 
controllable  frequency  and  controllable  voltage  ac  power  supply. 
This  motor  showed  about  the  same  rev  to  rev  speed  variations  as  the 
ISRL  induction  motor  system  and  showed  larger  speed  variations 
within  each  rev,  probably  due  to  the  four-pole  armature  and  a 
fairly  low  inertia  of  the  armature  and  shaft  system.  It  then 
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seemed  unlikely  that  we  could  improve  the  test  shaft  speed 
stability  by  replacing  the  motor-controller  system. 

6.2.2  Development  of  Speed  Correction  Algorithm 

Instead  of  replacing  the  motor  system,  we  decided  to 
record  the  16  PPR  signals  from  the  optical  shaft  position  encoder 
on  one  channel  of  our  data  system.  We  then  devised  a  method  to 
correct  the  blade  deflection  readings  through  the  use  of  this 
multipulse  per  rev  shaft  position  information.  Circuitry  to 
accomplish  the  recording  of  that  signal  was  designed,  fabricated, 
and  installed.  The  eight  test  runs  completed  are  listed  in  Table 
24.  Figures  64  and  65  show  the  speed  variations  for  seven 
consecutive  shaft  revolutions  of  Runs  1  and  2  of  non-excited  12- 
bladed  test  disk  at  1,200  and  1,800  rpm,  respectively,  as  indicated 
by  the  deviations  from  their  long  term  means  of  the  shaft  16  PPR 
and  1  PPR  data  readings.  Constant  shaft  speed  would  produce 
consecutive  identical  readings  (no  deviation  from  the  mean)  of  the 
pulse  series  from  each  shaft  position  encoder  within  the  accuracy 
of  the  encoder  and  recording  system.  These  two  plots  show  that 
shaft  speed  was  always  changing,  that  the  speed  changes  were  not 
large,  and  that  the  percentage  speed  changes  were  larger  at  1,200 
rpm  (+600  counts  out  of  500,000  or  ±0.12  percent)  than  at  1,800  rpm 
(+300  counts  out  of  333,333  or  +0.09  percent).  The  latter  would  be 
expected  because  rotary  Inertia  effects  increase  as  the  square  of 
rpm,  tending  to  reduce  accelerations  of  the  shaft  as  revolution 
speed  increases. 

Two  facts  were  readily  discernible  from  Figures  64  and 
65  when  we  realized  that  blade  tip  deflection  data  dispersions 
about  their  means  were  generally  no  greater  than  half  of  the  1  PPR 
readings  dispersion.  First,  no  measure  of  the  1  PPR  data  could  be 
adequate  for  correcting  the  blade  data  readings.  Evidently 
corrections  of  individual  blade  data  readings  must  be  made  based  on 
on  the  rotational  speed  of  the  particular  revolution  in  question. 
Secciid,  since  revolution  speed  was  constantly  changing,  as  shown  by 
curvature  in  all  the  plots  of  shaft  position  for  individual 
revolutions,  a  linear  slope  correction  for  each  revolution,  based 
on  the  1  PPR  reading  at  the  end  of  the  revolution,  would  not 
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TABLE  24 

SPEED  CORRECTION  TEST  RUN  OF  DETAILS 
16  PPR  SHAFT  ENCODER  RECORDED  IN  CHANNEL  1 
VIBRATION  MODES  STATIONARY  ON  THE  12-BLADED  DISK 
375  REVS  RECORDED  PER  RUN 


Run 

#1 

1,200  RPM 

No  Excitation 

Run 

#2 

1,800  RPM 

No  Excitation 

Run 

#3 

1,800  RPM 

2D  @  125  Hz 
Fex  =  185  Hz 

Run 

#4 

1,800 

ID  6  86  Hz 

Fex  =  116  Hz 

Run 

#5 

1,800 

3D  §  199  Hz 
Fex  =  289  Hz 

Run 

#6 

1,200 

3D  e  196  Hz 
Fex  =  256  Hz 

Run 

#7 

1,200 

1  D  @  83  Hz 
Fex  =  103  Hz 

Run 

#8 

1,200 

2D  e  121  Hz 
Fex  =  161  Hz 

151 


Figure  64.  Shaft  Speed  Variations  at  1,200  rpm. 
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Deviations  from  Mean  Levels  in  Counts 


RUM  2,  15  JAM  1988 
1800  RPM,  NO  EXCITATION 
1  PPR  MEAM  -  333,337  COUNTS 
SEVEN  SUCCESSIVE  REVS 
SHOWING  SPEED  CHANGES 
-t-  COUNTS  •  SLOWER  SPEED 
-  COUNTS  •  FASTER  SPEED 


adequately  correct  the  blade  deflection  readings  during  most 
revolutions.  Instead,  a  correction  was  required  based  on  the 
instantaneous  speed  error  (i.e.,  shaft  position  error)  when  the 
data  reading  was  made. 

Considering  the  above,  we  formulated  a  correction  to 
individual  blade  data  readings  based  on  the  location  of  the  mean  of 
each  blade  data  set  along  the  constant  speed  line  defined  by  the 
mean  of  the  1  PPR  data  set  and  the  deviation  of  the  shaft  speed 
from  that  constant  speed  line  as  defined  by  the  16  shaft  position 
encoder  and  1  PPR  reading  deviations  from  their  means. 

The  correction  was  performed  during  computer 
processing  of  the  recorded  data  set  and  is  described  below: 

^Bik^^*^Bik"“Bi^“^^SLik““sLi^'^^i  ^®SLik"*^SLi^  "  ^®SHik““sHi^ 


where 


V  =  corrected  data  value 

R  =  data  reading 

M  s  mean  of  the  data  set  readings 
Bi  =  Blade  number  i 

k  =  rev  number  k 

SLi  =  shaft  encoder  position  immediately  below  Bi 
SHi  =  shaft  encoder  position  immediately  above  Bi 
Fi  =  MSHi-MSLi  ~  interpolation  factor  for  Bi 


The  factor  Fi  provided  linear  interpolation  of  the  shaft  speed 
readings  for  the  encoder  (tachometer)  positions  last  before  and 
first  after  each  blade  mean  position  (Mg^^)  along  the  constant  speed 
line.  All  the  mean  values,  the  SLi,  SHi,  and  Fi  for  each  blade, 
and  the  standard  deviation  of  the  data  set  for  each  blade,  encoder 
position,  and  the  1  PPR  were  determined  in  a  preprocessing  run. 

Then  we  corrected  the  data  readings  and  performed  a  frequency 
analysis  of  the  data  set.  Linear  interpolation  between  the  shaft 
speed  values  at  adjacent  shaft  encoder  positions  appeared  valid  if 
the  straight  line  segments  between  the  positions  were  adequate  to 
closely  define  the  curvature  present  in  the  speed  line  for  each 
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revolution.  That  seemed  true  in  our  case,  even  at  the  low  average 
speed  of  1,200  rpms.  Note  that  we  had  17  line  segments  defined  by 
18  shaft  encoder  readings  which  were  then  zero  reset  after  the  1 
PPR,  the  16  PPR  readings,  and  the  1  PPR  reading.  Another  valuable 
feature  of  the  method  was  that  neither  the  blades  nor  shaft  encoder 
(tachometer)  positions  needed  to  be  spaced  nor  needed  to  be  assumed 
to  be  spaced  symmetrically  around  the  circle  of  rotation  because 
the  means  of  their  data  sets  located  them  accurately  in  time  of 
revolution  and  consequently  in  space. 

We  processed  a  number  of  test  run  data  sets  through 
the  preprocessing  and  data  correction  programs.  We  processed  the 
corrected  and  uncorrected  data  sets  through  a  statistical  routine 
that  computed  the  standard  deviation  (a)  values  of  each  data  set. 
The  a  values  for  Run  1  (1,200  rpm)  and  Run  2  (1,800  rpm)  are 
presented  in  Tables  25  and  26,  respectively.  Channel  1  contained 
the  16  shaft  encoder  values  and  channels  2,  3,  and  4  each  included 
plane-of-light  detector  readings  for  each  of  the  12  blades  on  the 
test  disk.  The  1  PPR  data  readings  were  included  as  the  last 
reading  of  each  revolution  for  each  channel.  The  corrected  o 
values  include  both  random  and  systematic  errors  due  to 
instrumentation  errors  and  noise,  unbalance  of  the  test  specimen, 
bearing  noise,  and  buffeting  of  the  blades  due  to  changes  in  air 
circulation  in  the  open  test  chamber,  but  errors  due  to  speed 
changes  of  the  test  shaft  were  eliminated,  or  nearly  so. 

Data  values  for  other  test  runs  3  through  8  were 
corrected  for  shaft  speed  variations.  The  subsequent  frequency 
analysis  made  on  the  corrected  data  sets  revealed  the  1,  2,  and  3 
diametral  mode  shapes  of  the  test  disk. 

6.3  STATIONARY  AND  TRAVELING  WAVE  DATA  ANALYSIS 

The  data  reduction  and  analysis  system  used  for  the  ISRL  NSMS 
data  of  this  effort  was  documented  fully  in  a  report,  "Data 
Analysis  Guide  for  the  Noncontacting  Blade  Deflection  Measurement 
System, "(7).  The  data  processing  procedures  are  described  briefly 
in  the  paragraphs  that  follow. 
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TABLE  25 

RUN  #1,  *1,200  RPM,  *NOT  EXCITED,  *2  VALUES 


Blade _ Chan  1 _ Chan  2 _ Chgin  3 _ ghan  4 


Raw 

Raw  Corrected 

Raw 

Corrected 

Raw 

Corrected 

1 

1 

3 

21 

11 

15 

11 

23 

6 

2 

2 

- 

22 

44 

6 

40 

32 

49 

8 

3 

3 

- 

42 

70 

7 

68 

34 

74 

8 

4 

4 

- 

61 

94 

7 

89 

34 

98 

10 

5 

5 

- 

80 

119 

5 

113 

34 

123 

7 

6 

6 

- 

99 

144 

7 

142 

34 

149 

14 

/ 

7 

- 

118 

169 

6 

164 

34 

174 

8 

8 

8 

- 

136 

194 

6 

187 

34 

199 

11 

9 

9 

- 

155 

219 

6 

214 

35 

223 

6 

10 

10 

- 

174 

243 

5 

237 

34 

247 

5 

11 

11 

- 

192 

268 

8 

262 

33 

273 

9 

12 

13 

- 

229 

13 

14 

- 

247 

14 

15 

- 

265 

1  PPR 

15 

16 

- 

283 

1  PPR  298 


TABLE  26 


Blade _ _ 

RUN 

_ Chan  1 

#2,  *1 

,800  RPM, 

Chan  2 

*NOT 

EXCITED,  *2 

Chan  3 

VALUES 

Chan  4 

Raw 

Raw 

Corrected 

Raw 

Corrected 

Raw  Corrected 

1 

1 

— 

1 

8 

7 

4 

4 

9 

6 

2 

2 

- 

6 

13 

5 

15 

11 

15 

6 

3 

3 

- 

12 

20 

5 

45 

34 

22 

7 

4 

4 

- 

17 

27 

6 

51 

34 

28 

7 

5 

5 

- 

23 

34 

4 

31 

31 

35 

5 

6 

6 

- 

28 

41 

5 

45 

32 

44 

10 

7 

7 

- 

33 

48 

4 

53 

33 

49 

6 

8 

8 

- 

38 

55 

4 

61 

33 

56 

8 

9 

9 

- 

44 

62 

4 

69 

34 

63 

4 

10 

10 

- 

49 

69 

5 

72 

33 

69 

4 

11 

11 

- 

54 

76 

6 

79 

34 

77 

6 

12 

12 

- 

59 

82 

4 

89 

35 

83 

5 

13 

13 

- 

65 

14 

- 

70 

15 

- 

75 

16 

- 

80 

1  PPR  84 

*  We  think  the  noisy  data  on  channel  3  was  due  to  wide  dispersion  of 
the  light  beam  from  the  transducer,  which  has  been  replaced  with 
our  spare  transducer  for  the  next  test  run  series. 
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6.3.1  NSMS  Data  Reduction  Procedures  for  Stationary  Wave 

Analysis 

UDRI  NSMS  data  reduction  software  was  designed  to  run 
on  a  VAX  computer.  To  perform  data  reduction  analysis,  the  user 
stores  test  NSMS  data  on  magnetic  tape.  The  tape  was  then  read  on 
the  UDRI  VAX,  and  the  data  from  the  recorded  runs  were  sorted  and 
cataloged  onto  disk. 

The  user  then  processed  these  data  for  analysis  and 
results-plottlng.  Terminal  access  to  the  VAX  for  these  purposes 
was  made  through  dial-up  connections  from  WPAFB  or  direct  terminal 
connections  at  UDRI. 

When  a  test  disk  was  excited  In  a  stationary  wave  mode 
on  the  disk,  the  NSMS  system  collected  data  for  blade  deflections 
at  each  of  three  plane-of-light  detectors.  These  detectors  were 
normally  connected  to  NSMS  system  channels  2,  3,  and  4.  If  each 
channel  was  working  correctly,  the  user  would  get  an  independently 
complete  set  of  data  c.i  eac^  of  the  three  channels,  and  he  needed 
to  process  only  one  of  these  to  analyze  the  specimen  response. 

For  this  type  of  test,  the  data  reduction  software 
first  would  read  in  the  channel  data  and  perform  a  statistical 
analysis  of  the  data.  The  program  then  corrected  the  channel  data 
for  shaft  speed  variations  and  repeated  the  statistical  analysis. 
The  correction  was  based  on  the  following  linear  speed  correction 
equation,  developed  by  UDRI; 


^Bik  ^  “  ^n-i  )  ”  (  ^cT  -iv  ”  ^ct  ^  ) 


Bik  “Bi 


SLik  SLi‘ 


Fi  *  [  (  R 


"  MqtM  )  -  (  I^SHik  "  **SHi  ^  ^ 


SLik  *SLi 


(88) 


Where 

Vgj^j^  =  Corrected  data  value  for  Blade  i,  at  rev  k. 

Rgj^j^  =  Recorded  reading  for  Blade  i,  at  rev  k. 

Mp.  =  Mean  value  for  Blade  i  over  all  good  recorded  revs. 
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Rg-..  =  Recorded  reading,  at  rev  k,  for  the  16  PPR  shaft 
encoder  position  immediately  preceding  Blade  i. 

Mgj,  .  =  Mean  value  over  all  good  revs  for  the  16  PPR  shaft 

encoder  position  Immediately  preceding  Blade  1. 

RgHij^  =  Recorded  reading,  at  rev  k,  for  the  16  PPR  shaft 

encoder  position  immediately  following  or  coincident 
with  Blade  i. 

MgHi  =  Mean  value  over  all  good  revs  for  the  16  PPR  shaft 

encoder  position  immediately  following  or  coincident 
with  Blade  i. 

^i  ''  =  interpolation  factor  for  Blade  i. 

”sHi  ■  “sLi 

The  shaft  encoder  data  was  stored  using  the  NSMS 
system  channel  1 . 

6.3.2  Time  and  Frequency  Analysis 

After  the  preliminary  reduction  was  completed,  a 
printable  report  of  the  statistical  analysis  was  generated.  The 
program  then  created  a  time  series  for  each  of  channels  2,  3,  and  4 
using  the  corrected  data,  and  then  performed  a  frequency 
analysis. 

The  frequency  analysis  consisted  of  a  Fourier  analysis 
of  the  time  series  generated  previously  for  one  of  the  data 
channels.  The  program  used  an  FFT  routine  from  the  ISML  library  to 
perform  the  Fourier  analysis.  The  ISML  library  was  a  commercial 
mathematics  and  statistics  software  package,  and  was  available  on 
both  the  UDRI  VAX  and  the  ASD  SEWS  ADA  VAX. 

As  a  final  step,  the  data  reduction  software  program 
writes  to  disk  any  requested  time  series  data  and  results  of  the 
frequency  analyses.  These  disk  files  can  be  used  for  plotting  time 
series  or  frequency  spectrum  data. 

An  example  of  a  typical  flow  chart  for  the  stationary 
wave  speed  correction,  time  and  spectrum  analysis  and  plotting 
software  program  is  shown  in  Figure  66. 
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CHECK.JBY 


Subdirectory  [U0REE02.NSMS.880l IS] 


SUtionery  Wove 
Aaelyait  Program 


NBDMS  data  file  from  tape: 
RS.003.880ilS.DAT:l 


Figure  66.  Standing  Wave  Software  Program  Flow  Chart 


6.3.3 


Anaiyglg 

When  a  test  disk  was  excited  in  a  traveling  wave  mode 
on  the  disk,  the  NSMS  system  again  collected  data  for  blade 
deflections  at  each  of  the  three  plane-of -light  detectors.  For 
this  test  mode,  however,  the  user  needed  results  from  each  of  the 
three  NSMS  data  channels  to  generate  an  appropriate  time  series  for 
the  vibration  analysis.  Also,  an  additional  low- response  test  run 
was  required  to  generate  a  nonvibrating  reference  data  set  for  the 
specimen's  magnet-transducer  configuration.  The  remaining  system 
channel  was  used,  as  before,  to  collect  the  16  PPR  shaft  speed 
variation  correction  data. 

For  the  traveling  wave  test,  the  nonvib  data  set  was 
first  analyzed  and  processed  as  described  previously.  The  data 
were  corrected  for  shaft  speed  variations,  and  for  channels  2,  3, 
and  4,  the  scaled  mean  of  each  blade's  passage  time  count  was 
stored.  That  is,  the  mean  of  a  blade's  passage  time  reading  over 
all  good  recorded  revolutions  was  divided  by  the  mean  1  pulse  per 
rev  count  value  (IPPR)  for  the  nonvib  run.  Since  the  nonvib  run 
and  its  corresponding  vibration  run  were  made  at  different  shaft 
speeds,  the  mean  IPPR  values  would  be  different  between  the  two 
data  sets  and  their  ratio  can  be  used  to  predict  the  nonvib 
reference  value  of  each  blade  at  the  test  condition  rpm  of  each 
rev. 

In  the  next  step,  the  traveling  wave  data  reduction 
software  read  in  each  channel's  data  from  the  vib  data  set.  The 
program  performed  a  statistical  analysis  of  the  raw  data.  The 
program  then  corrected  channel  data  for  shaft  speed  variations  and 
repeated  the  statistical  analysis.  The  corrected  data  for  each  of 
the  three  data  channels  then  was  referenced  to  the  nonvib  data  mean 
for  each  blade  as  adjusted  for  the  rev  speed  difference.  These 
corrected  differential  time  of  arrival  values  were  converted  to 
deflection  values  and  then  were  incorporated  into  time  history  data 
sequences  for  the  FFT  frequency  analysis  procedure,  as  in  the 
stationary  wave  analysis. 
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These  analysis  steps  were  implemented  in  the  traveling 
wave  software  program,  and  both  nonvib  and  vlb  data  sets  are 
Illustrated  in  the  flow  diagram,  Figure  67.  For  traveling  wave 
analysis  with  three  probes,  the  system  measured  response  at  three 
evenly  spaced  points  of  a  test  mode  wavelength.  This  arrangement 
used  three  sample  points  per  each  revolution,  and  the  Increased 
sample  rates  were  achieved  by  storing  data  from  symmetrically 
positioned  blades.  This  analysis  used  lists  of  symmetrically 
located  blades.  The  blade  number  arrangement  for  2,  3,  and  4 
diametral  modes  is  shown  in  Tables  27,  28,  and  29  respectively. 

6.3.4  Z-Grid  Data  Reduction  Procedure 

The  optical  system  traveling  wave  analysis  program  can 
also  process  the  Z-Grid  data.  The  grid  analysis  program  referred 
to  the  grid  pulses  as  "blade"  values,  but  these  were  data  sets  from 
a  single  blade  as  it  passed  across  36  grid  wires  in  each  of  the 
dual  grid  rings.  The  user  should  respond  to  the  software  prompt 
questioning  the  number  of  blades  as  36  while  analyzing  the  Z-Grid 
data. 

When  the  test  disk  was  excited  in  traveling  wave  mode, 
the  dual  Z-Grid  system  acquired  data  for  the  magnet  tipped  blade 
at  each  of  the  36  grid  wires  around  the  grid  ring.  These  data  were 
reordered  on  channels  3  ana  4  of  the  optical  data  acquisition 
system.  Also,  an  additional  low-response  data  test  run  was  made  to 
generate  the  nonvibrating  reference  data  set  as  in  the  optical 
system. 

For  the  Z-Grid  system,  the  nonvibrating  data  set  was 
first  analyzed  and  corrected  for  shaft  speed  variations.  To 
correct  for  speed  changes,  the  mean  of  a  Z-Grid  pulse  response  for 
500  revolutions  was  divided  by  the  mean  of  1  PPR  for  the  non¬ 
vibrating  run.  In  the  next  step,  the  traveling  wave  data  reduction 
software  performed  statistical  analysis  on  the  resonance  data  and 
corrected  it  for  shaft  speed  changes.  The  corrected  data  for  two 
Z-Grid  channels  takes  the  form  of  equation  (87)  in  Section  6.3.2. 
The  analysis  program  at  this  point  treated  each  Z-Grid  pulse  as  if 
it  were  a  blade  reading  on  a  test  disk  with  36  blades.  The 


161 


162 


TABLE  27 

2D  MODE  OPTICAL  PROBE  ARRAY  FOR  TIME  SERIES 

Transducers  placed  at  60-Degree  intervals  around  the  specimen. 
Symmetric  blade  ordering  uses  the  following  blade  arrays: 


Array  Set  1: 

Channel 

2 

3 

4 

Blades 

1 

7 

13 

20 

26 

32 

Array  Set  2: 

Channel 

2 

3 

4 

Blades 

2 

8 

14 

21 

27 

33 

Array  Set  3: 

Channel 

2 

3 

4 

Blades 

3 

9 

15 

22 

28 

34 

Array  set  4: 

Channel 

2 

3 

4 

Blades 

4 

10 

16 

23 

29 

35 

A  time  series 

sequence  : 

for  Set  4 

(as  < 

an  example)  would 

repeated  sequence  of 

the  data 

from  the 

following: 

Event  Number: 

1  2 

3  4 

5 

6 

Blade . : 

4  10 

16  23 

29 

35 

Channel . : 

2  3 

4  2 

3 

4 
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TABLE  28 

3D  MODE  OPTICAL  PROBE  ARRAY  FOR  TIME  SERIES 


Transducers  placed  at  40~degree  Intervals  around  the  specimen 
Symmetric  blade  ordering  uses  the  following  blade  arrays: 


Array  Set  1:  Channel  234 


Blades 

1 

5 

9 

13 

17 

21 

26 

30 

34 

Array  Set  2: 

Channel 

2 

3 

4 

Blades 

2 

6 

10 

14 

18 

22 

27 

31 

35 

Array  Set  3: 

Channel 

2 

3 

4 

Blades 

3 

7 

11 

15 

19 

23 

28 

32 

36 

Array  set  4: 

Channel 

2 

3 

4 

Blades 

4 

8 

12 

16 

20 

24 

29 

33 

37 

A  time  series 

sequence 

for  Set  4 

(as  an 

example) 

would  be  a 

repeated  sequence  of 

the  data 

from  the 

following: 

Event  Number: 

1  2 

3  4 

5  6 

7  8 

9 

Blade . : 

4  8 

12  16 

20  24 

29  33 

37 

Channel . : 

2  3 

4  2 

3  4 

2  3 

4 
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TABLE  29 

4D  MODE  OPTICAL  PROBE  ARRAY  FOR  TIME  SERIES 


Transducers  placed  at  30-degree  intervals  around  the  specimen. 
Symmetric  blade  ordering  uses  the  following  blade  arrays: 


Array  Set  1: 

Channel 

2 

3 

4 

Blades 

1 

4 

7 

10 

13 

16 

20 

23 

26 

29 

32 

35 

Array  Set  2: 

Channel 

2 

3 

4 

Blades 

2 

5 

8 

11 

14 

17 

21 

24 

27 

30 

33 

36 

Array  Set  3: 

Channel 

2 

3 

4 

Blades 

3 

6 

9 

12 

15 

18 

22 

25 

28 

31 

34 

37 

Array  set  4: 

Channel 

2 

3 

4 

Blades 

4 

7 

10 

13 

16 

19 

23 

26 

29 

32 

35 

38 

A  time  series 

sequence  for  Set 

4 

(as  an  example) 

would 

be 

repeated  sequence  of 

the  data 

from  the 

following: 

Event  Number:  1 

2  3  4 

5 

6 

7 

B  9  10 

11 

12 

Blade . :  4 

7  10  13 

16 

19 

23 

26  29  32 

35 

38 

Channel . :  2 

3  4  2 

3 

4 

2 

3  4  2 

3 

4 
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correction  and  the  analysis  formulae  for  the  Z-Grid  system  have  the 
same  form  as  in  the  optical  system. 

The  basis  for  configuring  the  dual  grid  wires  in  the 
letter  form  "Z"  was  to  enable  the  grid  wire  system  to  measure  axial 
and  tangential  components  of  the  blade  tip  deflection.  An  enhanced 
analysis  procedure  was  added  to  the  grid  wire  software  analysis  to 
resolve  blade  tip  axial  and  tangential  deflection  components  from 
the  dual  grid  system  blade  tip  time  of  arrival  data  sets.  The 
rationale  for  the  deflection  resolution  in  two  directions  is 
explained  below. 

If  the  blade  had  both  axial  and  tangential 
displacements,  their  resolution  into  independent  components 
required  the  application  of  two  grid  wires.  In  the  ISRL  grid 
system,  two  conducting  grid  wires  were  perpendicular  to  each  other 
and  at  45*  to  the  blade  tip  magnet  path.  This  arrangement  measured 
both  components  of  the  blade  tip  vibration.  If  the  blade  with  the 
source  magnet  experienced  only  axial  vibration,  then  the  time  of 
arrival  values  for  one  grid  would  be  equal  to  but  opposite  in  sign 
to  those  of  the  second  grid.  If  the  blade  experienced  only  the 
tangential  vibration,  the  time  of  arrival  values  for  both  grid 
wires  would  be  either  greater  or  smaller  by  equal  amount  depending 
on  the  phase  of  vibration  of  the  instrumented  blade.  If  the 
vibration  was  skewed  to  the  reference  axis,  the  two  orthogonal 
vibration  components  would  be  resolved. 

We  determined  the  axial  and  tangential  displacement 
components  of  the  instrumented  blade  using  our  digital  Z-Grid 
system.  To  resolve  these  components,  the  time  of  arrival  data  from 
the  two  grid  circuits  were  corrected  for  speed  changes.  The  sum 
and  difference  of  the  time  of  arrival  data  sets  were  then 
calculated  grid  wire  by  grid  wire.  When  these  two  data  sets  were 
added,  the  tangential  motion  of  the  blade  was  obtained.  Their 
difference  defined  the  axial  displacement  of  the  blade.  These 
calculation  procedures  were  incorporated  in  the  system  data 
analysis  software. 


This  method  of  analysis  was  applied  to  a  sample  data 
run,  and  the  results  from  this  analysis  matched  well  with  the 
respective  time  series  data  and  FFT  spectra  from  the  optical  system 
and  hence  confirmed  the  validity  of  the  analysis  procedure. 

The  results  from  a  comparative  study  of  the  optical  and 
Z-Grid  system  performance  for  the  integral  order  vibrations  of  UTRC 
test  disk  are  summarized  in  reference  (45) . 

6.4.  DIFFERENTIAL  DOUBLE  PULSE  INTERFEROMETRY  FRINGE  ANALYSIS 

The  double  exposure  method,  involving  two  holographic 
exposures  of  the  vibrating  disk  in  rapid  succession,  was  used  to 
record  disk  modal  patterns.  The  resulting  interference  fringe 
patterns  could  then  be  interpreted  to  provide  net  displacement  of 
the  disk  between  the  two  exposures.  This  amplitude  was  used  to 
determine  the  peak  to  peak  vibration  of  the  disk.  The  method  of 
calculation  is  explained  below. 

Consider  the  disk  vibrating  sinusoidally  in  its  normal  mode 
when  two  exposures  were  made  by  pulsing  the  ruby  laser 
(X  =  694  nm)  at  time  ti  and  tj ,  as  shown  on  the  acceleration  signal 
of  Figure  68a.  The  displacement  Z  at  time  t^  was 

Zj  =  A  sin  utj  (89) 

where  u  was  the  circular  frequency  (u  =  2tf)  and  A  was  the  half 
amplitude  vibration.  The  displacement  Z2  at  time  tj  is 

Z2  =  Asin  ut2  (90) 

The  net  displacement  is 

Z2-Z1  *  A(sin  wt2  -  sin  wti )  (91) 

but  ti  =0,  wti  =  0,  sin  wti  *  0,  and  Z2  *  0  for  integral 
order  traveling  wave  modes,  and  t2  was  the  pulse  separation 
time  (At) .  Then 
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Figure  68a.  Laser  Pulses  at  tj  and  tj  on  the 
Triggered  Accelerometer  Signal. 
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Figure  68b. 


Holographic  Plate 


Object  and  Reference  Be2un  Angles  with 
Reference  to  Hologram  Plane  in 
Double  Exposure  Method. 
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Zi  =  A  sin  (uAt) 


(92) 


Z2  - 


In  the  double  exposure  method  interference  fringes  occur  when 


Z  -  Z  -  fsn-ll  A 


(93) 


where  n  was  the  number  of  fringes  and  \  was  the  wave  length 


From  (91)  and  (92) 


A  = 


_ f2n-l)  X 

4  sin  (2irfAt) 


(94) 


The  out-of-plane  displacement  of  the  disk  must  be  corrected  for  the 
beam  separation  angle  (a  of  Figure  68b)  by  Then 


A  = 


i2nrl)  A 


4  sin  (2irfAt)  cos  a 


(95) 


and  peak  to  peak  vibration  was 


2A 


.I2n. 


2  sin  (2irfAt)  cos  a 
For  this  data  set 


(96) 


X  =  694  X  10“^  M  =  0.02732  mil 

—5 

At  =  10  seconds  -  10/isec 

cos  a  =  0.82  (a  =  34.92') 

F2P  =  135.4  Hz 

Fjp  =  182.0  Hz 

F^P  =  243.0  Hz 

and 

2A  *  (n-1/2)  - - 

sin  (6.2832xlo"^f) 

Then 

2A2JJ  =  (n-1/2)  3.9025  mils 


(97) 


(98) 
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(99) 


2A2JJ  =  (n-1/2)  2.9138  mils 

2A^jj  =  (n-1/2)  2.1823  mils  (100) 

The  equations  (98,  99,  and  100)  present  the  scaling  factors  to 
be  used  for  conversion  of  the  interferograms  to  peak-to-pealc 
deflection  in  mils  for  the  usual  cases  when  the  hologram  pulse 
interval  time  was  10  /isec.  They  were  calculated  assuming  the  sine 
equal  to  the  angle  in  radians  for  these  very  small  angles.  They  can 
be  scaled  to  other  pulse  intervals  by  dividing  the  scaling  constant 
by  [0.1  X  (pulse  interval  in  /xsec)  ]  for  any  of  the  three  modal 
vibration  patterns  used  in  the  test  series. 
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SECTION  VII 

BIADED  DISK  VIBRATION  INSTRUMENTATION  EVALUATION  TESTS 


Turbine  ar I  compressor  disks  vibrate  at  ix,  2X,  3X,  and  at 
higher  multiples  of  excitation  frequencies.  Synchronous  and 
asynchronous  excitation  conditions  of  these  disks  were  laboratory 
simulated  to  test  flexible  blades  and  bladed  disks  for  vibrations 
in  static  and  rotating  conditions.  Appropriate  hardware  and 
evaluation  tests  were  designed  to  study  the  modal  response  of  tuned 
and  mistuned  disks.  Suitable  instrumentation  schemes  for  the 
bladed  disk  modal  definition  were  implemented. 

This  section  describes  the  instrumentation  and  measurement 
systems  preparation  for  the  rotating  bladed  disk  vibration 
experiments.  Slip-ring  design  and  its  installation,  stationary  and 
traveling  wave  excitation  methods,  positioning  of  plane-of-light 
vibration  sensors  and  the  setting  up  of  the  dual  Z-grid  ring  were 
described.  Application  of  double  and  triple  pulse  laser 
interferometry  for  the  mode  shape  definition  of  the  rotating  bladed 
disk  assemblies  together  with  the  AEDC  NSMS  evaluation  test  series 
is  described. 

7.1  DISK  INSTRUMENTATION  FOR  VIBRATION  MEASUREMENTS 

7.1.1  Excitation  and  Response  Instrumentation  for  the  Test 

PjsKs 

Several  bladed  disks  were  fully  characterized  during 
this  contract  period  for  vibrations.  Static  and  rotating  disks 
were  excited  by  one  or  two  stationary  electromagnets.  Input 
voltage  and  the  frequencies  of  these  excitation  magnets  were 
controlled  by  a  combination  function  generator-power  amplifier  set. 
To  enable  higher  power  levels,  the  excitation  magnets  for  the 
rotating  disk  were  freon-cooled. 

Two  bladed  disks  -  a  12-inch-diameter,  12-bladed  disk 
and  an  11-lnch-diameter  and  38-bladed  disk  were  instrumented, 
tested  and  evaluated  for  vibration  characteristics  in  static 
conditions.  Resonances  were  identified  through  Lissajous  figures 
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on  the  oscilloscope  and  FFT  spectra  displays  on  the  spectrum 
analyzer.  Time  average  holography  recorded  several  diametral  and 
circumferential  mode  shapes  of  static  bladed  disks. 


During  the  test  runs  on  a  12-inch-diameter  disk,  the 
lead  wires  on  the  test  disk  accelerometers  and  the  leads  to  the 
four  channel  slip-ring  mounted  on  the  drive  shaft  behind  the  test 
disk  failed  when  the  spin  speed  exceeded  2,000  rpm.  Considerable 
time  was  lost  for  repair  and  reinstallation.  As  a  remedy,  we 
designed  and  installed  a  ten  channel  slip-ring  set  directly  on  the 
front  face  of  the  test  disk.  This  unit  has  ten  brass  rings  and  can 
accommodate  three  isolated  strain  gage  circuits,  two  accelerometer 
and  two  piezo  crystal  circuits.  Lead  wires  of  these  sensors 
were  laid  with  care  to  avoid  mistuning  mass  effects.  The  brush 
assembly  was  fabricated,  and  a  transfer  of  electrical  signals  from 
disk  sensors  to  the  stationary  signal  conditioning  and  analysis 
equipment  through  the  slip-ring  -  brush  assembly  was  ensured. 

The  drawing  for  the  ten-channel  slip-ring  components 
is  shown  in  Figure  69,  and  a  photograph  of  the  unit  installed  for 
use  is  shown  in  Figure  70.  Figure  71  shows  the  transducers 
installed  on  the  disk  to  use  the  new  slip-ring  circuits. 

7.2.  TRAVELING  WAVE  (INTEGRAL  ORDER)  EXCITATION  TESTS 

A  joint  effort  with  the  Unite  Technology  Research  Center 
(UTRC)  and  Arnold  Engineering  Development  Center  was  undertaken  to 
excite  and  measure  a  flexible  bladed  disk  traveling  wave  motion. 
This  excitation  condition  simulated  the  static  pressure  field 
induced  oscillations  on  a  rotating  turbine  disk. 

Discussions  with  AEDC  and  UTRC  engineers  on  the  evaluation  of 
their  NSMS  led  to  the  design  of  a  suitable  test  specimen  and  test 
setup.  The  test  specimen  had  to  meet  three  requirements:  first, 
UTRC  preferred  a  38-bladed  disk  with  modal  frequencies  in  the 
frequency  range  of  100-250  Hz  for  the  lower  order  modes;  second. 
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Figure  69.  Ten-channel  Slip-Ring  Components, 


gure  70.  Slip-Ring  iMounting  to  the  Disk  Hub. 


Figure  71.  Instrumentation  on  12-in. -Dia.  Disk 


the  disk  should  show  large  deflection  at  blade  tips  so  that  UTRC 
optical  probes  could  capture  these  deflections  with  ease;  third, 
the  size  of  the  test  specimen  was  not  to  exceed  I2~inch  diameter,  a 
size  that  was  acceptable  for  optical  holography. 

7.2.1  Test  Disk  Design 

Prior  experience  with  disk  dynamics  guided  UDRI  to 
select  the  approximate  diameter  and  thickness  for  the  proposed  test 
disk.  NASTRAN  based  finite  element  methods  were  used  to  design  the 
test  disk  to  satisfy  frequency-amplitude-size  requirements.  Four 
possible  test  disk  configurations  were  analyzed.  Their  natural 
frequencies  were  summarized  in  Table  30.  The  two  12 -inch-diameter 


TABLE  30 

NATURAL  FREQUENCIEi:  OF  FOUR  62.5-MIL-THICK  TEST  DISKS  (HZ) 

AT  ZERO  SPEED 


Disk  Dia(in) 

12.0 

12.0 

10.0 

11.0 

Blade  Length  (in) 

1.5 

3.0 

1.0 

2.0 

Mode 

ID 

59.29 

57.04 

90.82 

72. 

.56 

OC 

71.45 

63.83 

107.89 

83. 

.99 

2D 

87.16 

68.92 

133.98 

96. 

.34 

3D 

168.98 

101.86 

269.17 

166. 

.16 

4D 

— 

125.59 

434.17 

231. 

.26 

disks  produced  low  frequencies,  whereas  a  10-inch  disk  yielded 
higher  frequencies.  The  11-inch-diameter  disk  with  2-inch  blades 
was  selected  because  it  yielded  frequencies  close  to  UTRC 
specifications.  Considering  the  size  of  excitation  magnets,  blade 
length  was  later  changed  to  1-1/2  inches.  Disks  of  62.5-  and  50- 
mil  thickness  were  analyzed,  and  their  natural  frequencies  are 
presented  in  Table  31.  The  62 . 5-roil-thick  disk  had  a  4D  frequency 
of  294  Hz  at  zero  speed  and  would  have  required  much  larger 
excitation  force  to  produce  acceptable  blade  tip  deflections.  The 
disk  of  50-mil  thickness  was  selected  for  more  detailed  analysis. 
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TABLE  31 

NATURAL  FREQUENCIES  (HZ)  OF  62.5-  AND  50-MIL  THICK 
DISKS,  EACH  11-INCH  DIA.  WITH  38  BLADES  OF 
1.5-INCH  LENGTH  AT  ZERO  SPEED 


Mode 

62.5-Mil  Thickness 

50-Mil  Thickness 

ID 

72.5 

58.2 

OC 

85.8 

69.6 

2D 

103.3 

84.4 

3D 

195.0 

161.8 

4D 

294.1 

244.0 

A  reasonable  number  (38x8)  membrane  and  bending, 
constant  thickness  plate  elements  were  used  to  model  a  38-bladed, 
11-inch-diameter,  50-mil-thick  disk.  The  normal  mode  method  (Rigid 
Format  3)  of  analysis  was  used  to  determine  the  natural  frequencies 
and  modes  of  the  centrally  clamped,  nonrotating  bladed  disk.  The 
centrally  clamped  model  was  used  because  a  2-inch-dleuneter  ten- 
channel  slip-ring  assembly  was  directly  mounted  to  the  disk  center 
to  connect  disk  response  monitoring  accelerometers  and  piezo 
crystals  to  a  frequency  analyzer.  Two  cases  of  hub-clamping,  one 
without  a  slip  ring  but  centrally  clamped  over  a  1-lnch-diameter 
and  the  other  with  a  slip-ring  assembly  rigidly  fixed  at  the  2-inch 
diameter,  were  analyzed;  see  results  in  Table  32. 

To  reduce  computational  effort,  a  simple  assumption  of 
infinite  rigidity  was  assumed  for  the  1-  or  2-inch  diametral 
surface  at  the  disk  hub.  Natural  frequencies  of  the  11-inch 
diameter,  50-mil-thick  disk  with  the  slip  ring  were  within 
acceptable  limits.  Out  of  plane  normalized  vibration  contours  of 
the  disk  diametral  modes  were  determined  analytically. 

Static  bench  tests  to  establish  modal  characteristics 
were  performed  for  fabricated  50-mil  and  62.5-mil  thick  steel 
disks.  Time-average  holography  was  used  to  record  mode  shapes  of 


these  disks.  Verification  of  analytical  frequencies  and  inodes 
helped  us  to  continue  the  disk  analysis  for  the  rotational  effects 


TABLE  32 

HUB  SLIP-RING  EFFECT  ON  NATURAL  FREQUENCIES 
AT  ZERO  RPM  OF  50-MIL-THICK  DISK  OF  11-INCH  DIAMETER 
WITH  1.5 -INCH-LONG  BLADES 

Mode _ Without  Slip  Rina _ Wjt^  Slip  Rjpg 


ID 

58.2 

80.5 

OC 

69.6 

83.9 

2D 

84.4 

100.4 

3D 

161.8 

165.1 

4D 

244.0 

244.6 

7. 2. 1.1  Rotational  Effects  on  the  Test  Disk 
Frequencies 

NASTRAN's  differential  stiffness  feature, 

Rigid  Format  13,  was  used  to  analyze  the  disk  rotational  effects. 
Modal  frequency  increase  due  to  centrifugal  stiffening  during 
rotation  were  determined,  and  the  resulting  frequencies  are  shown 
in  Table  33.  Information  from  this  table  was  used  to  construct  the 
frequency-speed  diagram  to  determine  the  critical  speeds  of  the 
test  disk. 


TABLE  33 

ROTATIONAL  EFFECT  ON  DISK  NATURAL  FREQUENCIES  (HZ) 

DISK  GEOMETRY;  11-INCH-DIA. ,  50-MIL-THICK  38  EACH  OF  1.5-INCH- 

LONG  BLADES 

2-INCH-DIA.  HUB  SLIP  RING  UNIT 


SPEED 

IN  RPM 

Mode 

0 

1000 

2000 

3000 

4000 

5000 

2D 

100 

103 

109 

115 

130 

148 

3D 

155 

158 

163 

170 

185 

201 

4D 

225 

227 

234 

241 

250 

264 

178 


7. 2. 1.2  Excitation  of  a  Rotating  Disk  Due  to  a  Static 

Force  Field 

While  analyzing  the  dynamics  of  rotating 
turbine  disks,  we  need  to  discuss  important  contributions  made  to 
the  turbine  disk  dynamics  research  by  several  researchers.  Tobias 
and  Arnold  (46)  in  their  first  extensive  study  of  the  dynamics  of 
imperfect  disks,  set  up  simple  analytical  models  and  excited  a 
rotating  disk  by  a  single  stationary  magnet  for  measuring  critical 
resonances.  Ewins  (47)  extended  Tobias  and  Arnold  work  and 
developed  a  generalized  force  model  for  the  coupled  bladed  disk 
response  and  used  it  to  perform  a  numerical  study  on  a  simple 
mistuned  disk.  Ewins  (48)  also  analyzed  the  Individual  blade  and 
disk  participation  in  the  total  disk  response  of  tuned  bladed  disk 
assemblies.  Srinivasan  conducted  extensive  analytical  and 
experimental  investigation  of  the  blade-to-blade  response  of 
mistuned  turbine  disks  and  summarized  the  disk  dynamics  research  in 
a  survey  article  in  reference  (49) . 

To  understand  the  disk  vibrations  due  to  a 
static  force,  consider  a  point  force  F^  acting  on  the  disk  at  8  * 
9^.  This  point  force  can  be  represented  by  the  Fourier  series  as 

F  F 

F(0)  =  — ®  +  —  2*  cos  m(0-0  )  (101) 

2K  %  m=l  ^ 

The  generalized  force  in  each  mode  can  be  determined  from  the 
relation 


0J.,  the  position  of  static  force  was  time  variant  due  to  the  disk 
rotation  at  speed  Q.  Therefore, 

0^  =  0  -  fit  (105) 


and  letting  0=0  and  for  any  time  t  Eg.  (104)  becomes 


Qi 


cos 

sin 


nOt 


(106) 


From  this  equation,  it  was  evident  that  a  disk  rotating  past  a 
static  force  was  stimulated  for  a  specific  mode  of  vibration  at 
speed  fl . 

Considering  only  one  mode  of  vibration,  the 
generalized  excitation  force  was  Q  =  cosnOt. 

The  vibration  of  an  axisymmetric  disk  due  to  a 
stationary  concentrated  harmonic  force  cos  wt  can  be  written  as 

|g?®|  n0  cos  «t  ,  i  =  1,  2  (107) 


where  4^ ^  was  the  solution  of  the  equation  of  motion  of  disk  in 
generalized  coordinates;  was  the  amplitude  of  vibration  and  a 
function  of  radius,  excitation  force  and  magnification  factor;  0 
was  the  angular  coordinate  of  the  disk  relative  to  inertial  frame; 
cos  n0  and  sin  n0  were  the  shapes  of  vibrations  of  two  dual  modes; 
and  t  was  the  time. 

Letting  0^.  =  0  in  Eg.  (105),  we  get  the 
instantaneous  position  of  excitation  force  as 


0  =  nt 


(108) 


Substituting  Eg.  (108)  into  Eg.  (107)  and  from 
the  trigonometric  identities  cos  «t  cos  nOt  =  ^[cos(«+nfl)t  +  cos(u- 
nflt)  ]  and  cos  wt  sin  nOt  =  ^[sin(w+nfi)t  -  sin(w-nflt)  ]  Eg.  (107) 
becomes 


cos 

sin 


n0 


cos  (w+nfl)t 
sin  (w+nfljt 


(109) 
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From  Eg.  (109),  four  conditions  for  resonance  were  possible  for  an 
imperfect  disk.  These  are 

cos  (u+nQ)t  and  sin  (u+n0)t  (110a) 

and  cos  (a)-nQ)t  and  sin  ((i)-nfi)t  (110b) 

Eg.  (IlOa)  represents  two  forward  traveling  waves  whereas  Eg. 

(110b)  represents  two  backward  traveling  waves.  These  were  shown 
in  the  freguency-speed  diagram  (Figure  72),  of  the  disk. 

But  in  this  analysis,  the  disk  was  treated  as 
perfect  and  hence,  only  two  resonance  conditions  were  possible. 

By  extending  backward  traveling  wave  branches 
of  Eg.  (110)  to  intersect  the  speed  (0)  axis,  we  obtain  critical 
speeds  of  the  rotating  disk. 

Using  the  information  in  Table  33  and  the 
conditions  in  Eg.  (110) ,  we  determined  the  critical  speeds  for  our 
test  disk  from  Figure  73.  These  speeds  were  3,950  rpm  for  2D, 

3,640  rpm  for  3D,  and  3,900  for  4D.  At  these  speeds,  an  excitation 
force  with  zero  excitation  freguency  would  set  the  disk  in 
resonance . 

Assuming  nonpreferential  modes  for  the  disk, 
its  vibration  can  be  written  as 

0  =  a^^  cos  nQt  sin  wt  +  sin  n«t  sin(wt  +  c)  (111) 

and  letting  a^^  =  a2  =  a  and  €  =  ^  by  the  assumption  of  one  mode 
being  orthogonal  to  the  other.  Eg.  (Ill)  becomes 

^  =  a  sin  («-nfl)t  (112) 

This  was  a  backward  traveling  wave  at  speed 
nfl.  This  wave  appears  stationary  to  a  stationary  observer  in 
space.  The  experiments  were  designed  to  excite  the  backward 
traveling  wave  on  the  disk  due  to  a  static  force  field  in  space. 
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2s  Forward 


S  =  Sine  Mode 


C  =  Cosine  Mode 


Figure  72.  Frequency  -  Speed  Diagram  of  an  Imperfect  Disk. 
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11"  DIA.,  0.05"  THICK  DISK  WITH  38-1.5"  LONG 
BLADES  AND  2"  DIA.  HUB  SUPRING 

20  30  40 


FREQ.  (Hz) 
300 


Figure  73.  Frequency-Speed  Diagram  of  ll-in.-Dia., 
50-mil-Thickness  Tuned  Disk. 
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As  agreed  by  APL-UDRI,  AEDC-UTRC,  UDRI  designed  an 
excitation  system  to  excite  the  disk  in  its  2,  3,  and  4  diametral 
modes  and  measure  its  response. 

Stimulus  to  form  a  pure  mode  was  provided  by  placing 
permanent  magnets  at  desired  sine  wave  response  peaks.  Magnets  (4 
for  2D,  6  for  3D,  and  8  for  4D)  were  positioned  on  both  sides  of 
the  disk  at  wavelength  intervals  of  the  desired  mode-pattern  but 
offset  by  half  a  wavelength  on  the  front  and  back  sides.  Half  the 
number  of  excitation  magnets  were  positioned  in  the  front  while  the 
other  half  were  at  the  back  side  of  the  test  disk.  Aluminum  bars, 
shown  in  Figure  74  were  used  to  support  the  front  magnets  from  a 
circular  yoke.  The  back  magnets  were  held  by  a  plate  fixed  to  the 
test  chamber  floor  and  motor  support  frame  behind  the  disk.  A 
general  front  magnet  holder  mounting  configuration  is  shown  in 
Figure  75  for  the  three  test  modes.  A  specific  excitation  scheme 
for  the  2D  mode  is  illustrated  in  Figure  76. 

7.2.3  Probe  Positions 

In  these  tests,  plae-of-llght  sensors  were  used  to 
measure  blade  tip  deflections  of  the  rotating  disk.  Since  the 
speed  of  the  rotating  disk  and  the  speed  of  the  response  wave  were 
the  same,  a  stationary  probe  measured  nearly  constant  amplitude  as 
the  disk-blades  passed  the  probe.  In  effect,  one  probe,  provided 
only  a  constant  amplitude  each  time  a  blade  passed  the  probe,  so 
this  arrangement  masks  the  actual  dynamic  response  of  the  disk. 
Therefore,  more  than  one  probe  was  necessary  to  define  the  blade 
and  disk  traveling  wave  vibrations  in  integral  order  modes. 

Four  AEDC  and  three  UDRI  plane-of-light  sensors  were 
used  here  to  define  modal  deformation  of  the  test  disk.  The 
following  paragraphs  describe  some  considerations  for  probe 
locations  to  monitor  three  diametral  modes  of  the  disk. 

Positioning  of  three  UDRI  probes  to  measure  2D,  3D,  and 
4D  mode  shapes  is  shown  in  Figure  77.  A  photograph  in  Figure  78 
shows  UDRI  optical  probes  in  the  test  setup.  For  2D  mode 
excitation,  probes  were  spaced  60  degrees  apart  in  space  around  the 
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TDC 


18  Tapped  Holes, 
i“20  NC  X  1  Deep  in 


Yoke 


Figure  75.  Excitation  Magnet  Positions. 
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Probe  Locations  on  Mountinq  Rinq  (Lookinq  Aft) 
0«2D,  x*3D,  A*4D  Probe  Locations 
1/4"  Mounting  Holes  at  7"  Rad. 

Angular  Spacing  for  2D  Probes,  60^ 

Angular  Spacing  for  3D  Probes, 

Angular  Spacing  for  4D  Probes,  30 


Figure  77.  Plane -of -Light  Probe  Positions. 


188 


Figure  78.  UDRI  Optical  Probes  in  the  Test  Setup. 


perimeter  of  the  test  disk,  while  for  the  3D  and  4D  modes,  the 
probes  were  spaced  at  40  and  30  degrees,  respectively.  The  probe 
location  requirement  was  to  space  the  three  probes  120*  apart 
within  the  modal  pattern  of  each  response  so  a  standard  Fourier 
transform  routine  based  on  uniformly  spaced  time  samples  could  be 
used  to  analyze  the  optical  data. 

The  AEDC  NSMS  measuring  system  required  four  blade 
deflection  measurements  within  a  modal  wavelength.  The  transducer 
locations  need  not  be  symmetrically  spaced  within  the  wavelength, 
but  the  locations  must  be  known  relative  to  each  other.  The  four 
deflection  amplitude  readings  were  used  to  define  the  blade 
response  equation  of  the  known  mode  shape  to  yield  the  modal 
amplitude  and  frequency. 

Figure  79  shows  the  UTRC  plane-of-light  probe  mounting 
blocks.  The  four  probes  are  mounted  in  the  test  setup  (Figure  80) . 

7.2.4  Instrumentation  Rina  fYoke^  Design.  Fabrication  and 

Assembly 

A  yoke  was  used  to  mount  the  UD  and  UTRC  optical  probes 
and  the  front  excitation  magnets.  The  yoke,  shown  in  Figure  81, 
was  designed  by  UTRC  to  mount  their  sensors  and  modifications  to 
the  design  for  our  requirements  were  added  by  UDRI.  After  the  yoke 
was  fabricated  by  the  Air  Force  machine  shop,  it  was  mounted  within 
the  test  chamber  (Figure  82) .  Mounting  brackets,  shown  in  Figure 
83,  held  the  yoke  in  position.  Tolerances  provided  in  the  holes  of 
these  mounting  brackets  allowed  the  yoke  to  be  positioned 
concentrically  with  the  test  disk.  Concentricity  within  3  mils  was 
established  using  a  dial  gage  mounted  directly  on  the  test  disk  by 
slowly  rotating  the  disk  and  adjusting  the  position  of  the  yoke 
with  respect  to  the  disk. 

7.2.5  Test  Disk  Vibration  Monitoring  Instrumentation. 

Excitation  System  and  Preliminary  Tests 

The  test  disk  was  instrumented  with  three  miniature 
accelerometers  and  two  piezo  crystals.  The  Instrumented  disk  with 
its  sensing  elements  and  their  slip-ring  connections  is  shown  in 
Figure  84.  The  specimen  and  test  monitoring  system  consisted  of 
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Figure  79.  UTRC  Probe  Mounting  Blocks. 
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Through  Slots 


Figure  81.  Instrumentation  Support  Ring  (Yoke). 


Figure  82.  Assenbly  of  Instrumentation  R 
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in  the  Test  Chamber. 


Figure  83.  Mounting 
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the  instrumented  disk,  the  slip  ring  unit,  preamplifiers,  a 
spectrum  analyzer,  and  several  oscilloscopes. 

Preliminary  tests  were  conducted  to  excite  the  test 
disk  by  permanent  magnets  at  the  disk  critical  speeds.  Half-lnch- 
dlameter  Samarium-Cobalt  permanent  magnets  were  installed  to  an 
existing  mounting  plate  and  were  positioned  perpendicular  to  the 
aft  disk  surface  at  a  3.75-inch  radius  from  the  spin  axis. 

Spin  tests  were  conducted  with  four  magnets  mounted  at 
90  degree  spacing,  two  magnets  mounted  at  180-degree  spacing,  and 
with  a  single  magnet  installed.  The  accelerometer  signal  spectra 
for  three  typical  runs  are  shown  in  Figure  85.  The  magnets  were 
widely  spaced  in  comparison  to  their  diameter,  and  tended  to  excite 
all  the  lower  order  modes  of  the  bladed  disk  test  specimen.  Figure 
86  shows,  from  top  to  bottom,  vibration  spectra  generated  by  four, 
two,  and  one  magnets.  The  4-D  disk  system  mode  was  marked  by  the 
bright  spot  at  the  top  of  the  spectral  peak  in  the  range  of  260  to 
293  Hz  for  the  three  spectra  shown. 

It  can  be  seen  in  the  spectra:  that  all  magnet 
configurations  excite  the  lower  four  modes  of  the  system;  that  a 
single  magnet  excites  all  modes  of  the  system,  much  like  an  impact, 
but  the  odd-order  modes  were  excited  relatively  more  strongly  than 
the  even-order  modes;  and  that  for  excitation  by  two  or  four  evenly 
spaced  magnets,  the  even-order  modes  were  excited  more  strongly 
then  the  odd-order  modes,  except  for  the  1-D  mode.  In  addition  to 
the  characteristics  shown  in  the  spectral  plots,  holographic 
interferograms  were  produced  which  showed  that  a  particular  mode 
could  be  emphasized  by  rotating  the  specimen  in  the  magnetic  field 
configuration  at  a  speed  where  that  mode  became  dominant.  This 
test  series  showed  that  we  could  excite  backward  traveling 
vibration  waves  on  our  test  disk  with  the  magnets  on  hand.  From 
these  tests  we  observed  that  a  fine  speed  control  was  required  to 
generate  a  pure  stationary  wave.  To  circumvent  this  difficulty  and 
to  excite  a  selected  mode  with  ease,  we  decided  to  install 
excitation  magnets  on  both  sides  of  the  disk. 
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Figure  84.  Instrumentation  of  ll-ln.-Dia.  Disk. 


Accelerometer  1 


Accelerometer  2 


Crystal  1 
(Malfunction) 


Crystal  2 

(Peak  offscale  at  lOV) 


Figure  85.  Accelerometer  Output  Signals. 
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Figure  86.  Response  Spectra  Due  to  Magnetic  Excitation 


7.3.  UDRI  AND  UTRC  PLANE-OF-LIGHT  SYSTEMS  EVALUATION  TESTS 


An  AEDC  NSMS  Evaluation  test  series  was  conducted  during  the 
week  11-15  April  1988.  While  the  disk  was  running  at  constant 
speed  with  peak  response  of  each  modal  pattern,  Indicated  by 
peaking  of  the  signals  from  the  accelerometers  and  crystals  mounted 
on  the  disk  as  displayed  on  the  spectrum  analyzer  screen,  UDRI 
recorded  blade  tip  deflection  data  on  the  ISRL  NSMS  data  system; 
while  UTRC,  as  the  AEDC  NSMS  contractor,  recorded  blade  tip 
deflection  data  on  the  AEDC  NSMS  data  system.  Also,  a  triple-pulse 
holographic  interferogram  (50)  was  recorded  on  film  from  the  ISRL 
optical  instrumentation  system.  UDRI  reduced  the  ISRL  NSMS  data 
and  compared  it  to  the  AEDC  NSMS  data  and  the  holographic  data. 

7.3.1  First  Test  Series 

In  the  first  test  run,  six  permanent  magnets  spaced  at 
60  degree  intervals  on  alternating  sides  of  the  disk  at  the  3.75- 
inch  disk  radius  and  spaced  5/16  inch  from  each  surface  produced  a 
3D  mode  fixed  in  space  (a  backward  traveling  3D  mode  on  the  disk  at 
the  disk  rotational  speed)  with  peak  disk  response  occurring  at 
3,732  rpm.  The  photograph  in  Figure  87  shows  the  rotating  disk 
deflected  in  3D  standing  wave.  The  test  run  was  recorded  by  UDRI 
with  zero  response  reference  runs  recorded  first  at  1,500  and  3,000 
rpm  for  use  in  the  data  reduction  programs.  UTRC  recorded  a 
resonance  data  set  and  then  asked  to  record  spin-up  and  spin-down 
runs  through  resonance  as  additional  data  sets.  When  these  runs 
were  made,  we  found  that  a  broad  resonance  peak  occurred  more  or 
less  continuously  from  3,700  to  3,850  rpm  with  a  normal  onset  at 
the  lower  speed  but  an  abrupt  dropoff  at  the  higher  speed. 

Resonance  onset  then  could  not  be  achieved  until  the  rotational 
speed  was  reduced  to  about  the  midpoint  of  the  resonance  speed 
range . 

With  the  magnet-disk  gaps  at  5/16  inch,  a  number  of 
holograms  were  made  of  the  3D  integral  mode  at  3,720  rpm.  UTRC 
then  recorded  an  AEDC  NSMS  data  set  at  3,720  rpm  followed  by  a 
reference  data  set  at  3,720  rpm  with  the  excitation  magnets 
removed . 
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Figure  87.  Rotating  Disk  Deflected  in  3D  Mode. 


201 


For  the  next  test,  four  magnets  were  installed  at  90- 
degree  intervals  and  5/ 16-inch  spacing  from  the  test  disk  to  excite 
the  20  integral  mcJe.  We  found  that  if  the  disk  was  accelerated 
very  slowly,  the  resonance  could  be  extended  from  3,600  to  4,300 
rpm  before  it  collapsed.  UTRC  then  recorded  the  20  blade 
deflections  at  '’,920  rpm  on  the  AEOC  NSMS  and  recorded  a  reference 
data  set  at  speed  with  the  magnets  removed.  Since  the 
inter ferograms  recorded  earlier  for  30,  have  displayed  unresolvable 
number  of  fringes,  we  decided  to  reduce  the  disk  response  amplitude 
by  reducing  the  excitation  force  due  to  stationary  magnets.  The 
magnets  then  were  set  at  7/16  inch  from  the  disk  and  the  20  mode 
was  recorded  on  the  AEOC  NSMS,  on  the  ISRL  NSMS,  and  on  a  series  of 
triple-pulse  holograms  at  3,780  rpm.  The  response  amplitude  was 
decreased  nearly  half  by  this  lower  excitation  force,  and  the 
resonant  frequency  and  response  frequency  range  of  the  disk  also 
were  reduced,  but  the  resonance  response  characteristic  remained 
the  same.  Eight  magnets  then  were  installed  at  7/16  inch  from  the 
disk  to  excite  the  40  mode.  Again,  the  same  response 
characteristic  of  the  disk  occurred.  The  response  peak  seemed 
highest  at  3,630  rpm,  and  data  were  taken  there  on  the  two  NSMS 
systems  and  on  holographic  film. 

Six  magnets  were  installed  at  7/16  inch  from  the  disk 
to  repeat  the  30  mode  test  at  the  same  excitation  level  used  for 
the  20  and  40  modes  the  previous  day,  which  had  seemed  to  produce 
the  desired  nuiciber  of  fringes  on  the  holograms.  Then  the  hologram 
reference  beams  were  also  adjusted  to  remove  the  apparent  fringe 
muddling.  The  30  modal  data  then  was  recorded  by  the  three 
instrumentation  systems.  Again,  we  found  that  the  resonance 
frequency  could  be  sustained  over  a  range  of  about  300  rpm  with  a 
sharp  collapse  and  instantaneous  phase  shift  of  the  response  at  the 
upper  speed  (or  frequency)  limit  of  the  resonance,  suggesting  a 
nonlinear  disk  response. 

Ouring  the  course  of  this  testing  program,  UORI 
recorded  response  data  on  the  ISRL  NSMS  for  the  3D  mode  with  5/16- 
inch  magnet  spacing  and  for  the  2D,  4D,  and  3D  modes  with  7/16-inch 
magnet  spacing.  For  each  of  those  records,  reference  data  were 
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taken  iniined lately  before  or  lnunediately  after  the  modal  data  at 
1,500  and  3,000  rpm  where  no  disk  response  was  indicated  by  the 
monitoring  instruments  on  the  disk  or  in  the  recorded  reference 
data.  The  ISRL  NSMS  data  records  are  listed  in  Table  34. 

Similarly,  data  were  recorded  on  the  AEDC  NSMS  for  these 
conditions,  but  the  reference  data  for  these  records  were  taken  at 
the  same  rpm  as  the  resonance  data  with  the  excitation  magnets 
removed.  No  accel  or  decel  data  were  taken  on  the  AEDC  NSMS 
because  of  the  nonlinear  response  characteristic  of  the  test  disk. 
We  were  later  informed  by  the  project  monitor  that  none  of  the 
holographic  data  was  useful  except  that  which  was  recorded  for  the 
3D  mode. 


TABLE  34 

AEDC  NSMS  EVALUATION  TESTS  -  FIRST  TEST  SERIES 
ISRL  NSMS  DATA  SYSTEM  TEST  RECORDS 
APRIL  12-15,  1988 


Date 

Run  # 

Mode 

Magnet 

Gap 

(Inch) 

RPM 

Data  Type 

Remarks 

4/12 

lA 

3D 

5/16 

1500 

Reference 

Speed 

1 

No  AEDC  Data 

2A 

3D 

3000 

Reference 

Speed 

2 

ISRL  Data  OK 

3A 

3D 

3732 

Resonance 

Data 

No  Phase  Records 

4/14 

4A 

2D 

7/16 

3780 

Resonance 

Data 

ISRL  &  AEDC 

5A 

2D 

3000 

Reference 

Speed 

1 

Data  OK 

6A 

2D 

1500 

Reference 

Speed 

2 

No  Phase  Records 

4/14 

7A 

4D 

7/16 

1500 

Reference 

Speed 

1 

ISRL  &  AEDC 

8A 

4D 

3000 

Reference 

Speed 

2 

Data  OK 

9A 

4D 

3630 

Resonance 

Data 

Small  Phase 

Change 

4/15 

lOA 

3D 

7/16 

1500 

Reference 

Speed 

1 

llA 

3D 

3000 

Reference 

Speed 

2 

12A 

3D 

3560 

Resonance 

Data 

No  Real 

Resonance 

7. 3. 1.1  Nonlinear  Response  Studies  of  the  Disk 

As  mentioned  previously,  one  of  the  plans 
during  the  April  test  series  was  that  UTRC  wanted  to  accomplish  two 
types  of  data  collection  and  data  reduction  with  the  AEDC  NSMS 
instrumentation  system.  The  first  of  these  was  for  steady-state 
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data  collected  at  the  peak  modal  response  for  the  2D,  3D,  and  4D 
resonant  modes  of  the  test  disk.  The  second  was  for  data  collected 
during  spin-up  and  spin-down  through  each  resonance  response,  the 
so-called  accel  and  decel  test  runs.  However,  during  testing  we 
found  that  the  test  disk  displayed  nonlinear  response  in  that  a 
normal  onset  of  resonance  with  a  gradual  phase  shift  was  not 
displayed  during  spin-up  to  resonance,  but  that  an  abrupt  collapse 
of  response  with  an  accompanying  immediate  phase  shift  occurred  as 
speed  was  increased  through  the  resonance  response  range.  Then  as 
speed  was  decreased  back  through  the  resonance  range,  an  abrupt 
response  onset  with  an  immediate  phase  shift  occurred  at  a 
considerably  lower  spin  speed  within  the  spin-up  resonance  response 
speed  range.  This  resonance  characteristic  is  shown  in  Figures  88 
and  89  for  2D  and  3D  response  modes  for  small  magnet  gap  settings. 

We  thought  that  the  nonlinear  response  of  the 

test  disk  was  a  function  of  the  strength  of  the  force  field  of  the 

excitation  magnets.  A  number  of  test  runs  then  were  made  with 
magnet  spacing  from  the  test  disk  varied  from  5/16  inch  to  9/16 
inch,  but  the  same  nonlinear  type  of  response  was  shown  during  all 
test  runs.  UDRI  then  was  asked  to  try  air- jet  excitation  of  the 
test  disk.  This  was  accomplished  in  early  May  of  1988  with  two 

small  air  jets  supported  at  3/8-inch  spacing  from  each  side  of  the 

disk  to  excite  the  2D  modal  response.  The  results  of  this  test, 
shown  in  Figure  90,  indicated  the  same  nonlinear  response. 

Subsequently,  at  a  meeting  between  UDRI  and 
APL/POTC  personnel,  it  was  decided  to  attach  stiffener  plates  of 
2.25-inch  diameter  and  0.115-inch  thickness  to  each  side  of  the  hub 
area  of  the  instrumented  test  disk  in  the  expectation  that  this 
would  improve  the  linearity  of  the  modal  responses.  UDRI 
fabricated  the  plates,  removed  the  slip  ring  hub,  installed  the 
plates  with  high-strength  HySol  EA934  epoxy  adhesive,  and  re¬ 
installed  and  rewired  the  slip  ring  hub  set.  Then  tests  were 
conducted  to  determine  if  the  reinforced-hub  test  disk  possessed 
nonlinear  behavior  in  nonrotating  conditions.  The  test  disk  was 
hard  mounted  and  was  driven  by  an  electromagnet  exciter  at  a 
constant  level  of  excitation  voltage  but  with  different  gaps 
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Figure  90.  Disk  Nonlinear  Response  Under  Air 
Jet  Excitation. 
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between  the  electromagnet  and  the  test  disk.  The  static  test 
results  shown  In  Figures  91  and  92  indicated  that  the  disk  behaved 
like  a  soft  spring.  Variation  of  response  frequency  with  amplitude 
of  excitation  force  decreased  with  an  increase  in  the  mode  number 
of  vibration.  For  example,  leftward  shift  of  the  response 
frequency  peak  with  higher  excitation  force  was  more  pronounced  for 
ID  as  compared  to  the  3D  and  4D  modes. 

After  the  static  vibration  tests  were 
completed,  the  disk  was  reinstalled  in  the  spin  test  setup  and  was 
excited  in  the  2D  mode  at  various  magnet  gap  spacings  from  the 
disk.  The  results  of  five  successive  runs  in  which  the  magnet  gaps 
were  increased  from  3/8  inch  to  5/8  inch  in  1/ 16-inch  increments 
indicated  that  the  disk  now  displays  a  linear  response  in  the  2D 
mode  if  the  magnet  gaps  were  9/16  inch  or  greater,  as  seen  in 
Figure  93.  The  response  amplitude,  however,  was  reduced 
appreciably  at  the  larger  magnet  gaps  and  the  accelerometer  output 
amplifier  gain  was  doubled  when  the  magnet  gaps  increased  from  1/2 
inch  to  9/16  inch. 

7. 3. 1.2  Triple  Pulse  Holographic  Interferometry 

Problem  and  Its  Diagnosis 

One  of  the  requirements  introduced  by  UTRC  for 
the  triple-pulse,  double-image  holograms  was  exact  measurement  of 
the  pulse  spacing  between  the  first  and  second  pulses,  which 
produce  the  first  interferogram,  and  the  second  and  third  pulses, 
which  produce  the  second  interferogram.  The  exact  times  were 
needed  to  establish  the  frequency  and  phasing  of  the  two  images. 

The  desired  time  intervals  were  10  fisec  between  the  successive 
pulses.  UDRI  monitored  some  triple-pulse  events  with  a  high 
frequency  oscilloscope,  but  pulse  intervals  could  be  measured  with 
confidence  only  to  1  fisec  accuracy  on  the  oscilloscope  screen. 

However,  during  the  study  of  the  pulse  interval 
timing  problem,  we  found  that  the  Appollo  ruby  laser  often  produced 
spurious  pulses.  Sometimes  as  many  as  five  or  six  pulses  occurred, 
as  shown  in  the  oscillograph  screen  pictures  of  Figure  94.  Careful 
study  of  the  three  storage  circuits  of  the  triple-pulse  system 
showed  a  shorted  silicon  controlled  rectifier  (SCR) ,  a  missing 
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Figure  93. 


Figure  94.  Multi-Pulse  Operation  of  the  Ruby  Laser. 
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diode,  and  a  thyratron  that  sometiiaes  broke  down  spontaneously  at 
normal  operating  voltage,  all  in  the  energy  storage  circuit  for  the 
third  laser  pulse.  Obviously,  the  multipulses  occurred  because  of 
the  faults  in  that  circuit.  We  (a)  checked  the  power  storage,  the 
voltage  sensing,  and  pulse-timing  circuits  of  the  laser  power  head; 
(b)  identified  and  replaced  a  shorted  silicon  controlled  rectifier 
and  missing  diode;  and  (c)  reduced  the  high  voltage  level  on  the 
third  Pockels  cell  to  an  acceptable  level.  These  efforts  had 
resulted  in  suppressing  the  spurious  pulses  of  the  laser,  but  the 
system  seldom  produced  the  desired  triple-pulse  train.  Note  that 
at  every  stage  of  modification,  output  of  the  sensing  photo  diode 
and  the  input  to  the  internal  Pockels  cell  were  monitored  and 
displayed  on  a  digital  storage  oscilloscope.  We  were  convinced 
from  these  oscilloscope  records  that  the  changes  incorporated  into 
the  circuitry  of  the  laser  power  supply  unit  had  altered  the 
thresholds  of  the  three  firing  thyratrons,  and  a  new  set  of  optimum 
voltage  settings  needed  to  be  established  for  the  Pockels  cells  to 
produce  a  triple-pulse  train  with  a  specified  time  of  separation 
between  pulses. 

Phone  conferences  were  held  with  Mr.  Ralph  Page 
of  the  Appollo  Co. ,  and  he  suggested  that,  contrary  to  what  was 
prescribed  in  the  system  operations  manual,  triple-pulsing  was  more 
easily  achieved  by  first  adjusting  the  third  pulse  circuit,  then 
the  second  and  the  first,  instead  of  the  specified  ascending  order. 
He  also  said  a  success  rate  of  greater  than  30  percent  was  not  to 
be  expected  for  the  triple-pulse  system.  By  following  his 
directions,  UDRI  was  able  to  achieve  a  triple-pulse  success  rate  of 
about  25  percent  for  either  10  or  25-/isec  pulse  intervals. 

After  restoring  the  system  to  its  triple  pulse 
operation  with  pulse  separation  accuracy  of  1/is  as  displayed  in 
Figure  95,  several  triple-pulse  (sequential  double  pulse)  holograms 
were  tried  while  the  test  disk  was  in  resonance.  The  three  pulses 
were  achieved  by  sequential  switching  of  the  laser  head  three 
internal  and  the  two  external  Pockels  cells.  But  it  was  discovered 
that  even  when  no  third  pulse  was  displayed  on  the  pulse  monitoring 
scope,  the  reconstructed  image  contained  spurious  fringes.  This 
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fact  led  us  to  believe  that  the  external  Pockels  cells  in  the 
reference  path  had  a  poor  extinction  ratio  and  allowed  light 
passage  when  it  was  intended  to  completely  block  a  laser  pulse. 

This  leakage  of  laser  light  through  one  of  the  external  Pockels 
cells  combined  with  the  hologram  due  to  the  two  intended  pulses 
produced  an  effect  similar  to  Moire.  One  such  sample  inter ferogram 
was  enclosed  in  Figure  96  to  show  the  difficulty  in  interpreting 
such  holograms.  An  attempt  to  quench  the  leakage  by  placing  a 
combination  polarizer-filter  in  front  of  the  Pockels  cells  proved 
to  be  a  failure  as  the  same  effect  was  seen  in  the  interferograms 
produced  by  the  system.  It  was  concluded  at  this  stage  that  the 
trouble  was  in  the  delay  lines  and  pulse  forming  networks  of  the 
external  Pockels  cells  control  voltages,  causing  them  to  attain 
inconsistent  and  inadequate  extinction  levels,  or  in  the  Pockels 
cells  themselves. 

Meanwhile,  UTRC  had  agreed  to  utilize  double 
pulse  holograms  (51)  to  determine  the  modal  deformation  of  the 
resonating  disk  and  accordingly  requested  us  to  proceed  to  record 
double  pulse  interferograms  for  the  test  disk  while  in  resonance. 

7.3.2  Second  Test  Series;  Even  Probe  Positions 

The  second  test  series  was  required  for  the  following 
reasons:  (a)  sequential  double-pulse  or  triple-pulse  holograms 

recorded  on  the  first  test  series  were  unacceptable  due  to  Moire 
type  fringes  in  interferograms,  (b)  AEDC  NSMS  data  analysis  had 
produced  unsatisfactory  results  for  two  out  of  the  three  disk 
response  modes,  and  (c)  the  occurrence  of  abrupt  collapse  of  the 
disk  response  modes  without  gradual  phase  change  muddled  the  UTRC 
single  degree  of  freedom  analysis  method. 

For  this  second  test  series,  the  linear  dynamic 
characteristics  of  the  test  disk  had  been  improved  by  installing  a 
2 . 25-inch-dia.  x  0.11-inch-thick  stiffener  on  each  side  of  the  disk 
in  the  hub  area.  Also  wider  spacing  was  used  between  the  disk  and 
the  excitation  magnets. 

Test  matrix  objectives  and  test  conditions  for  this 
second  series  remained  the  same  -  excitation  magnets  were 
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positioned  on  alternate  sides  of  the  disk  on  circles  of  3.75  inches 
radius  to  excite  a  specific  node,  and  ISRL  optical  probes  were 
positioned  with  angular  spacings  of  60,  40,  and  30  degrees  to 
measure  the  2D,  3D,  and  4D  node  responses,  respectively. 

Test  dates  with  run  numbers  and  test  conditions  were 
provided  in  Table  35. 

The  first  test  in  this  series  was  with  magnet  gap 
spacing  of  5/8  inch  from  the  disk  and  with  25  ns  for  triple-pulse 
separation.  Interferograms  corresponding  to  Runs  #1B  and  #2B 
contained  an  unresolvable  number  of  fringes,  indicating  large 
displacements.  The  gap  was  increased  to  3/4  inch  in  Runs  #3B  and 
#4B.  AEDC  NSMS  on  site  data  analysis  by  UTRC  yielded  no  useful 
data  because  of  high  signal  to  noise  ratios  in  the  optical  probe 
data.  Threshold  levels  for  their  trigger  device  was  altered,  blade 
tips  were  refurbished,  and  the  Argon  laser  output  was  trimmed.  The 
subsequent  test  runs  (#5B,  #6B)  were  useful.  Triple-pulse 
holography  was  pursued  until  now.  Later  a  decision  was  made  to 
discontinue  triple-pulse  efforts  and  record  double-pulse  holograms 
at  resonance,  before  resonance,  and  after  resonance.  The  off- 
resonance  holograms  would  provide  qualitative  phase  information  as 
the  disk  transited  through  resonance. 

Double  pulse  holograms  were  made  for  4D  (Runs  7B 
through  lOB  and  17B  and  18B) ,  for  3D  (Runs  IIB  through  16B) ,  and 
for  2D  (Runs  19B  through  22B)  with  a  variety  of  exciter  spacings 
and  pulse  separation  times.  Two  data  sets  for  every  test 
condition,  one  at  resonance  and  the  other  at  off-resonance  speed 
(usually  at  3,000  rpm  for  ISRL  and  resonance  speed  without  magnets 
for  UTRC)  were  acquired.  Tests  were  conducted  at  various  magnet 
gap  spacings  from  9/16  to  15/16  inch  to  meet  AEDC  NSMS  data  range 
requirements . 

7 . 3 . 2 . 1  AEDC  Disk  Spin  Up  and  Spin  Down  Test 

At  the  end  of  the  standard  test  series  that 
acquired  response  of  the  disk  in  its  resonance  modes,  spin  up  and 
spin  down  tests  were  conducted  for  the  AEDC  NSMS  one  sensor  method. 
Response  amplitude  of  a  disk  accelerometer  was  monitored  while  the 
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TABLE  35 

AEDC  NSMS  EVALUATION  TESTS  -  SECOND  TEST  SERIES 
ISRL  DATA  RECORDS 
JULY  27  -  AUGUST  5,  1988 


Magnet 

Gap 


Date 

Run  # 

Mode 

(Inch) 

RPM 

7/27 

IB 

4D 

5/8 

3648 

2B 

4D 

3000 

7/28 

3B 

4D 

3/4 

3646 

4B 

4D 

3000 

7/29 

5B 

4D 

3/4 

3640 

6B 

4D 

3000 

8/1 

7B 

4D 

5/8 

3644 

8B 

4D 

3000 

8/2 

9B 

4D 

5/8 

3636 

lOB 

4D 

3648 

8/3 

IIB 

3D 

5/8 

3642 

12B 

3D 

3632 

13B 

3D 

3646 

14B 

3D 

3000 

8/4 

15B 

3D 

3/4 

3633 

16B 

3D 

3000 

8/4 

17B 

40 

9/16 

3642 

18B 

4D 

3000 

8/5 

19B 

2D 

15/16 

4063 

20B 

20 

15/16 

4058 

21B 

2D 

3000 

22B 

2D 

4058 

Data  Type 

Remarks 

Resonance  Speed 

ISRL  Data  OK 

Reference  Speed 

AEDC-No  Data 
Holograms  bad 

Resonance  Speed 

ISRL  Data  OK 

Reference  Speed 

AEDC-No  Data 
Holograms  bad 

Resonance  Speed 

ISRL  Data  OK 

Reference  Speed 

AEDC  Data  OK 
Holograms  bad 

Resonance  Speed 
Reference  Speed 

Good  Data 

Lower  Half  Power 

Point  Speed 

Upper  Half  Power 
Point  Speed 

Good  Data 

Resonance  Speed 

Lower  Half  Power 
Point  Speed 

Upper  Half  Power 
Point  Speed 
Reference  Speed 

Good  Data 

Resonance  Speed 
Reference  Speed 

Good  Data 

Resonance  Speed 
Reference  Speed 

Good  Data 

Resonance  Speed 
Near  Resonance 
Speed 

Reference  Speed 
Reference  Speed 
Without  Magnets 

Good  Data 
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disk  transited  through  resonance.  Amplitude  changes  due  to  speed 
change  In  steps  of  1  rpm  on  either  side  of  resonance  were  recorded 
for  the  2D,  3D  and  4D  modes.  Phase  change  while  the  disk  went 
through  resonant  speed  was  observed  by  monitoring  the  accelerometer 
output  signal  as  triggered  by  the  IPPR  signal.  A  smooth  phase 
change  of  nearly  90  degrees  was  observed  for  the  4D  mode  when  the 
magnet  gaps  were  3/4  Inch.  This  data  set  was  not  acquired  during 
the  first  test  series  due  to  nonlinear  behavior  of  the  disk. 

7.3.3  Third  Test  Series:  Uneven  Probe  Positions 

These  tests  were  requested  by  the  AEDC  to  verify  the 
accuracy  of  an  algorithm  used  by  UTRC  for  uneven  probe  positions 
around  the  circumference  of  the  bladed  disk.  Since  large 
displacements  were  the  prime  concern,  large  excitation  forces  were 
achieved  by  placing  the  excitation  magnets  close  to  the  disk. 

After  the  test  setup  and  data  acquisition  systems  were 
prepared,  11  test  runs  as  detailed  In  Table  36  were  made  on  the 
following  day.  For  each  test,  both  the  magnet  and  AEDC  probe 
settings  were  changed  to  the  desired  positions,  and  the  test  disk 
was  quickly  spun  up  to  Its  resonance.  Data  were  concurrently 
recorded  by  the  UDRI  NSMS  with  Its  probes  at  60,  40,  and  30  degrees 
for  the  2D,  3D,  and  4D  modes,  respectively. 

Holographic  Images  for  this  test  series  were 
not  recorded  due  to  the  film  transport  shutter  malfunctioning. 

7.4.  MODE  SHAPES  OF  ROTATING  DISKS  BY  LASER  INTERFEROMETRY 

The  double  exposure  laser  Interferometric  technique  was 
employed  to  record  the  mode  shapes  of  the  rotating  disk.  The  least 
possible  pulse  gap  of  10  microseconds  between  the  two  consecutive 
pulses  was  required  to  provide  a  minimum  countable  number  of 
fringes  on  the  surface  of  a  vibrating  bladed  disk  (52) .  Two  sets 
of  tests  were  conducted  using  the  triple  pulse  rxiby  laser  available 
ISRL.  This  technique  provided  quantitative  Interpretation  of  the 
Interferometric  fringes  due  to  resonant  vibrations  of  the  test 
object  (53)  and  served  as  a  comparative  tool  for  the  Z-grld  and 
plane~of-llght  sensor  systems. 


TABLE  36 

AEDC  NSMS  EVALUATION  TESTS  -  THIRD  TEST  SERIES 
UNEVEN  PROBE  POSITIONS 
AUGUST  24,  1988 


Run  # 

Mode 

Magnet 

Gap 

(Inch) 

RPM 

Data  Type 

Remarks 

1C 

4D 

9/16 

3648 

Resonance 

Speed 

ISRL  and  AEDC 

2C 

4D 

9/16 

3000 

Reference 

Speed 

Data  Acquisition 

3C 

4D 

9/16 

3631 

Resonance 

Speed 

OK  on  ISRL 

4C 

3D 

7/16 

3600 

Resonance 

Speed 

No  Data  on  AEDC 

5C 

3D 

7/16 

3000 

Reference 

Speed 

1 

6C 

3D 

No 

3600 

Reference 

Speed 

2 

Magnets 

7C 

4D 

No 

3631 

Reference 

Speed 

2 

Magnets 

for  Run  2C 

8C 

2D 

1/2 

3000 

Reference 

Speed 

1 

9C 

2D 

1/2 

4019 

Resonance 

Speed 

1 

IOC 

2D 

1/2 

4044 

Resonance  Speed 

2 

lie 

2D 

No 

4019 

Reference 

Speed 

Magnets 

for  Run  9C 

The  double  pulse  interferometric  method  consists  of  comparing 
two  consecutively  recorded  images  of  a  vibrating  disk  rotating  at  a 
steady  speed.  The  configuration  of  the  optical  setup  used  is 
illustrated  in  Figure  97.  The  rotational  motion  of  the  object  was 
compensated  by  an  image  derotator,  running  exactly  at  half  the 
speed  of  the  object-disk.  Circular  polarization  effects  of  the 
object  beam  were  compensated  by  the  addition  of  two  quarter  wave 
plates  before  and  after  the  image  derotator  system.  A  10- 
microsecond  interval  between  two  pulses  of  the  ruby  laser  was 
achieved  by  a  careful  selection  of  appropriate  voltage  settings  on 
the  Q-Switch  control  of  the  power  source  of  the  laser  and  the 
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occurrence  of  double  pulses  was  achieved  by  the  instrumentation 
system,  shown  in  Figure  98. 

An  important  feature  of  the  experiment  was  to  initiate  the 
first  pulse  of  the  laser  at  any  point  of  time  during  a  vibrating 
cycle  of  the  bladed  disk.  An  ideal  point  would  be  when  a 
harmonically  vibrating  disk  traverses  the  zero  velocity  point  on 
the  response  cosine  wave.  For  this  purpose,  a  filtered 
acceleration  signal  was  used  as  a  reference,  and  a  marker  signal 
from  the  charge  unit  of  the  laser  power  supply  was  made  to  override 
the  reference  sine  wave.  It  was  possible  to  shift  this  marker  to 
any  point,  and  double  pulse  holographic  images  of  the  vibrating  and 
rotating  disk  were  made  for  marker  positions  of  45,  0,  and  90 
degrees  on  the  response  sine  wave.  This  special  control  feature 
for  initiating  the  first  pulse  of  the  ruby  laser  was  achieved  by 
the  scheme  shown  in  Figure  99.  After  the  instrumentation  was 
prepared,  several  diametral  modes  of  the  test  disk  at  its  critical 
speeds  with  different  magnet  gaps  were  recorded. 

7.5  INTEGRATED  Z-GRID  -  NSMS  INSTRUMENTATION  TESTS 

The  instrumentation  ring,  which  supported  dual  Z-grid  rings 
and  three  plane-of-light  sensors,  was  mounted  concentric  with  the 
11-inch-diameter,  50-mil-thick  test  disk.  The  grid  system  source 
magnet  and  its  counter  weight  were  cemented  to  diametrically 
opposite  blade  tips,  with  a  gap  of  60  mils  from  the  grid  wires. 

A  series  of  tests  were  conducted  using  the  integrated  Z-grid 
-  NSMS  system  to  define  integral  vibrations  of  the  test  disk  in  2D, 
3D,  and  4D  modes. 

Test  runs  and  test  conditions  are  summarized  in  Table  37. 

The  sequence  for  the  integrated  system  tests  consisted  of: 

(a)  setting  up  the  excitation  magnets  to  excite  a  specific 
disk  mode; 

(b)  positioning  the  optical  probes  at  appropriate  positions 
in  the  support  ring; 
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Figure  98.  Double  Pulse  Settings  for  the  Ruby  Laser. 


SA  -  Spectrun  Analzer  KF  -  Kronhite  Filler 

CPO  -  Oscilloscope  ^1'  ^2  -  Anplifiers 

^1'  ^2  '  ■^‘^®^®i^'cmeters 


Figure  99.  Pulse  Trigger  Control  and  Display  Instrumentation. 


(c)  ensuring  the  satisfactory  functioning  of  optical  and  grid 
wire  probes  from  histogram  display;  and 

(d)  recording  a  data  set  for  each  system. 

When  the  disk  was  stabilized  in  its  resonance,  the  plane-of- 
light  sensors  (ch.  2,  3,  and  4)  collected  the  blade  deflection 
data.  Also,  an  additional  low  response  run  of  3,000  rpm  was 
recorded  as  reference  data.  Then,  optical  probes  were  disconnected 
from  the  signal  conditioning  unit  and  the  two  channels  of  the  Z- 
grid  system  were  connected  to  it  as  channels  3  and  4.  A  reference 
and  resonance  data  set  were  recorded  for  the  grid  system  under  the 
same  test  conditions.  The  data  were  analyzed  later  on  the  VAX 
computer . 

TABLE  37 

THE  INTEGRATED  Z-GRID  NSMS  SYSTEM  EVALUATION  TESTS 

Magnet 

Gap 


P^p  # 

Mode 

(inch) 

RPM 

_ Data  T^ 

tag _ 

— 

17ZG 

3D 

3/4 

3000 

Reference 

Data 

for 

Z-Grid 

18ZG 

3D 

3/4 

3657 

Resonance 

Data 

for 

Z-Grid 

19ZO 

3D 

3/4 

3657 

Resonance 

Data 

for 

NSMS 

20ZO 

3D 

3/4 

3000 

Reference 

Data 

for 

NSMS 

21ZO 

2D 

15/16 

4063 

Resonance 

Data 

for 

NSMS 

22ZO 

2D 

15/16 

3000 

Reference 

Data 

for 

NSMS 

23ZG 

2D 

16/16 

3000 

Reference 

Data 

for 

Z-Grid 

24ZG 

2D 

25/16 

4063 

Resonance 

Data 

for 

Z-Grid 

Runs  21  through  24  were 
were  not  retrievable. 

repeated 

as  the  magnetic  tape  data 

29ZO 

4D 

9/16 

3646 

Resonance 

Date 

for 

NSMS 

30ZO 

4D 

9/16 

3000 

Reference 

Data 

for 

NSMS 

31ZG 

4D 

9/16 

3000 

Reference 

Data 

for 

Z-Grid 

32ZG 

4D 

9/16 

3646 

Resonance 

Data 

for 

Z-Grid 

SECTION  VIII 
RESULTS  AND  DISCUSSION 


This  research  and  development  program  developed  and 
implemented  advanced  finite  element  based  deterministic  and 
probabilistic  methods  to  analytically  predict  the  normal  modes  of 
flexible  bladed  disks.  An  instrumentation  system  using  either  a 
four  plane-of-light  and  dual  Z-grid  measurement  head  was  developed, 
and  the  blade  tip  deflections  of  a  rotating  bladed  disk  in  backward 
traveling  wave  resonance  were  measured.  Full-field  double  pulse 
interferometric  method  served  as  reference  for  comparing  the  blade- 
tip  deflections. 

8.1  COMPARISON  OF  NATURAL  FREQUENCIES  OF  STATIC  AND  ROTATING  DISKS 

Two  flexible  disk  natural  frequencies  were  analytically 
predicted.  SAP,  MAGNA,  PROTEC,  and  NASTRAN  programs  (the  first 
three  for  static  disks  and  the  last  for  both  static  and  rotating 
disk  analysis)  were  used. 

The  zero  speed  normal  mode  data  was  first  derived  from  static 
bench  tests  and  compared  with  normal  mode  analysis  results  from  the 
NASTRAN  and  the  other  finite  element  methods.  All  analytical 
methods  had  predicted  double  frequencies  for  perfectly  axisymmetric 
or  tuned  circular  bladed  disks.  Small  errors  in  twin  mode 
frequencies  were  essentially  due  to  numerical  computational  errors 
in  the  eigenvalue  extraction  procedures.  Of  the  existing 
eigenvalue  extraction  methods  in  NASTRAN,  only  the  Fast  Eigenvalue 
Extraction  method  (FEEM)  yielded  quick  (less  CPU  time/mode  for  the 
same  size  problem)  and  reliable  double  modes  compared  to  other 
methods.  SAP  and  MAGNA  analysis  methods  were  slow  compared  to  any 
of  the  NASTRAN  methods  for  the  same  accuracy  level.  Note,  however, 
that  the  Mindlin  plate  theory  based  PROTEC  analysis  produced 
analytical  frequencies  which  were  closer  to  the  experimental  values 
for  static  disks.  A  12-bladed  12-inch-diameter,  62.5-mil-thick 
steel  disk  and  a  38-bladed  11-in. -diameter,  50-mil-thick  steel  disk 
were  analyzed,  and  their  static  frequencies  are  compared  in  Table 
38.  The  agreement  was  better  in  the  former  than  the  latter  because 
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the  clamping  area  for  the  11-inch  diameter  disk  was  over  a  larger 
surface  area  due  to  a  hub  slip  ring.  This  resulted  in  a  smaller 
active  annulus  area  for  the  disk.  The  frequency  discrepancy  varied 
from  4  to  8  percent,  and  the  errors  in  modal  frequencies  were 
essentially  due  to  the  simplistic  assumption  of  infinite  rigidity 
assumed  for  the  multichannel  slipring  assembly  for  the  11-inch  disk 
over  a  2-inch  diameter  hub  area.  A  large  concentrated  mass  due  to 
the  slipring,  near  the  disk  rotational  axis  also  lowered  the 
natural  frequencies  of  the  static  disk. 

TABLE  38 

COMPARISON  OF  STATIC  DISKS  RESONANCE  FREQUENCIES 


38-Bladed  11”  Dia,  50-mil 

12  Bladed  12"  Dia.  5/64”  Thickness  Disk  with  a  2"  Dia. 


Mode 

Thickness  Disk 

Analytical  Experimental 

Hub  Slip-Ring  and 
Analytical 

0.1”  Hub  Stiffness 

Experimental 

ID 

77 

69 

80.5 

— 

OC 

92 

85 

83.9 

— 

2D 

110 

105 

100.4 

95 

3D 

200 

194 

165.1 

150 

4D 

277 

270 

244.0 
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Resonances  for  rotating  disks  occur  differently  for  different 
load  types.  Under  harmonic  excitation,  standing  waves  as  a  result 
of  the  superposition  of  forward  and  backward  traveling  waves  build 
up  on  spinning  disk,  and  these  waves  travel  around  the  disk  with 
respect  to  the  disk.  For  example,  if  a  rotating  disk  was  excited 
harmonically  at  the  frequency,  p  =  w  +  nQ,  where  u  was  the  disk 
resonance  frequency,  n  was  the  nodal  diameter  number,  Q  was  the 
rotational  speed,  the  +  sign  for  the  forward  traveling  wave  and  the 
-  sign  for  the  backward  traveling  wave,  the  disk  will  resonate  at 
the  frequency  w.  This  type  of  disk  excitation  was  achieved  by  our 
freon-cooled  electromagnets  in  vacuum.  These  resonances  for  the 


227 


12-inch  disk  were  compared  with  the  NASTRAN  (Rigid  Format  13) 
predicted  frequencies  and  the  comparison  was  good.  Table  39. 


TABLE  39 

COMPARISON  OF  STATIC  DISKS  RESONANCE  FREQUENCIES 


Mode 

12-Bladed  12"  Dia. 

at  4,000 

Analytical 

5/64" 

rpm 

Thickreco  Disk 

Experimental 

ID 

104 

108 

OC 

110 

117 

2D 

143 

153 

3D 

226 

227 

4D 

298 

298 

A  38-bladed  disk  was  designed  for  the  integral  order  vibration 
tests.  The  disk,  in  these  tests  was  excited  by  a  static  force 
field.  From  the  equation  p  =  w  +  nfl,  when  the  excitation  frequency 
p  =  0,  i.e.,  a  static  force,  w  =  nfl  and  only  the  backward  traveling 
waves  would  be  present  in  the  disk.  These  resonant  waves  travel  at 
the  same  speed  as  the  rotating  disk.  The  speed  Q  is  also  called 
the  critical  speed.  Integral  order  waves  were  excited,  measured, 
and  compared.  The  resonances  were  monitored  by  accelerometers, 
piezocrystals  and  strain  gages.  Spin  tests  were  conducted  to 
establish  the  disk  critical  speeds  and  were  observed  to  be 
dependent  on  the  excitation  force  levels.  A  comparison  of  the 
critical  speeds  was  provided  only  when  the  disk  exhibited  linear 
response;  i.e.,  when  the  magnet  gap  was  9/16  inch  or  greater  from 
the  rotating  disk  surface.  The  analytical  and  experimental  results 
are  compared  in  Table  40.  The  measured  critical  frequencies  for 
the  second  and  third  diametral  modes  provided  good  agreement  with 
the  analytical  prediction.  An  error  of  1.5  percent  for  the  2D,  3D, 
and  5.5  percent  for  the  4D  mode  was  present.  Low  measured 
frequencies  were  possibly  due  to  overestimation  of  centrifugal 
force  effects  for  the  disk  with  a  slipring  mass  at  the  hub.  The  ID 
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critical  speed  mode  was  not  measured  because  the  lower  branch  of 
the  traveling  wave  would  intersect  the  speed  axis  at  a  speed  much 
higher  than  the  design  speed  of  the  drive  shaft  system.  We  decided 
not  to  run  the  tests  at  these  high  rotational  speeds. 


TABLE  40 

NATURAL  FREQUENCIES  (HZ)  OF  38-BLADED  DISK  WITH  2 . 0-INCH-DIA. 
HUB  SLIP  RING  AT  ZERO  RPM  AND  AT  CRITICAL  SPEED 


Zero 

Soeed 

Critical  Frequencies 
_ (Speeds)  in  RPM _ 

Mode 

Analvtical 

Exoerimental 

Analvtical 

Experimental 

2D 

100.4 

95 

132  (3950) 

135  (4058) 

3D 

165.1 

150 

185  (3640) 

184  (3670) 

4D 

244.6 
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260  (3900) 

251  (3760) 

8.2  MODE  SHAPE  ANALYSIS 

Static  disk  mode  shapes  were  first  determined  by  the  time 
average  holographic  method  in  which  image  exposure  time  exceeded 
the  period  of  oscillation  by  a  minimum  factor  of  two.  Good 
correlation  between  the  predicted  and  measured  modes  was  observed. 
Spatically  orthogonal  dual  modes  were  identified  for  both  test 
disks,  but  only  one  set  of  the  analytical  and  experimental  dual 
mode  shapes  of  the  38-bladed  disk  are  enclosed  in  Figures  100a  and 
100b. 

For  the  mode  shape  definition  of  rotating  disks,  double  pulse 
interferograms  were  recorded  when  the  disk  was  in  resonance. 
Interferograms  of  the  disk  in  2,  3,  and  4D  modes  are  shown  in 
Figure  101.  Figure  102  shows  typical  displacement  patterns  for  the 
3  and  4  diametral  modes  for  two  excitation  cases,  recorded  in  lO^s 
pulse  separation  time.  The  peak-to-peak  vibration  amplitude  of  the 
harmonic  wave  was  derived  by  sine  extrapolation  of  fringes.  The 
increase  in  fringe  number  for  a  larger  exposure  time  was  evident 
and  typical  holograms  for  2D  mode  for  10  and  15  /is  pulse  separation 
were  shown  in  Figure  103. 
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3D  Mode  -  162  Hz  4D  Mode  -  244  Hz 

Lgure  100a.  Analytical  Mode  Shapes  of  the  Static  Test  Disk. 


23 


Figure  100b.  Experimental  Mode  Shapes  of  the  Static 
Test  Disk. 
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C  2  ^ 


3D  MODE 


4D  MODE 


Inter ferograms  for  the  3D  and  4D  Modes  from  Two 
Excitation  Force  Levels. 


Figure  102. 
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Figure  103.  Inter ferograms  for  the  2D  Mode  for  10-  and  15-ys  Pulse  Separation. 


8.3  BLADE  TIP  DEFLECTIONS 


Quantitative  blade  tip  deflections  of  resonating  disks  were 
determined  by  AEDC  and  UDRI  NSMS  optical  methods  and  compared  with 
those  from  interferometric  methods.  Results  from  a  limited  Z-Grid 
test  data  were  also  analyzed. 

8.3.1  The  NSMS  and  Interferometric  Methods 

Maximum  deflections  at  the  blade  tips  were  determined 
from  (a)  point-wise  measurements  from  AEDC  and  UDRI  NSMS-Systems, 
and  from  (b)  wholefield  interferogram  data  recorded  when  the  disk 
had  stabilized  in  resonance. 

UDRI  data  analysis  programs  had  determined  the  maximum 
blade  tip  deflections  from  the  software  generated  time  series  and 
their  spectral  components  for  the  2,3,  and  4D  modes.  Although  a 
large  number  of  time  series  and  FFT  plots  (an  approximate  total  of 
110)  were  generated,  only  a  few  plots  are  included  in  this  report 
for  illustration.  Figures  104a  and  104b  show  the  NSMS  signals  in 
time  and  frequency  domains  for  the  3D  mode  in  the  large  amplitude- 
first  test  series.  Figures  105  through  107  show  the  time  series 
and  FFT  components  for  the  2,  3,  and  4D  modes  in  second  test  series 
in  which  the  test  disk  exhibited  linear  response.  The  plots  for 
the  third  test  series  results  are  not  presented  in  this  report  for 
want  of  space. 

A  summary  of  results  analyzed  by  the  UDRI  NSMS  data 
analysis  software  for  the  first,  second,  and  third  test  series  was 
provided  in  Tables  41,  42,  and  43,  respectively.  In  our  analysis 
procedure,  symmetrical  arrays  of  blades  were  considered  for  3  probe 
data.  Software  constructed  time  series  data  and  FFT  components  of 
the  data  series  yielded  nearly  identical  peak-to-peak  amplitude  of 
vibrations  for  blade  tips  of  the  resonating  test  disk.  The  2D,  3D, 
and  4D  modes  for  all  excitation  test  conditions  were  analyzed. 

AEDC  analysis  yielded  results  for  the  first  and  second  test 
series  and  a  summary  of  test  conditions  and  results  are  provided  in 
Tables  44  and  45.  No  results  were  obtained  for  the  third  test 


235 


Rub  12A/11A,  lS-Apr-88,  3-D  Bi  3560  rpa 
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Run  19B/aiB«  e5-Rug-88>  2D  ni  4663  rpm 


Figure  105a.  Time-Series  Data  for  the  2D  Mode  -  Second  Test  Series 


Run  19B/S1B*  OS-Aug-88*  2D  ni  4B63  rpa 


Figure  105b.  Fourier  Components  for  the  2D  Mode  -  Second  Test  Series 


Figure  106a.  Time-Series  Data  for  the  3D  Mode  -  Second  Test  Series. 


series  because  of  a  malfunction  in  the  AEDC  data  acquisition 
system.  The  first  test  series  yielded  large  vibration  amplitudes 
because  the  excitation  magnets  were  located  close  to  the  disk 
surface.  A  comparison  of  the  results  from  this  test  series  is 
provided  in  Table  46. 

TABLE  41 

BLADE  TIP  DEFLECTIONS  (P-P) 

CALCULATED  FROM  ISRL  NSMS  DATA 
FOR  THE  FIRST  TEST  SERIES 


Mode 

Gap 
(in. ) 

Speed 

(rpm) 

UDRI 

(p-p  in  mils) 

3D 

7/16 

3720 

110±2 

3D 

7/16 

3560 

52±2 

3D 

5/16 

3732 

160±4 

2D 

5/16 

3970 

No  Data 

2D 

7/16 

3780 

176±5 

4D 

7/16 

3630 

62±2 
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TABLE  42 


BLADE  TIP  DEFLECTIONS  (P-P) 
CALCULATED  FROM  ISRL  NSMS  DATA 
FOR  THE  SECOND  TEST  SERIES 


Res/Ref. 

Data  Resonance 

Record  RPM 

Magnet 

Gap 

rinchl 

Mode 

Mils  Defl.  -  P-P* 

1B/2B 

3648 

5/8 

4D 

25.5  ±  3.1 

3B/4B 

3646 

3/4 

4D 

15.1  ±  2.7 

5B/6B 

3640 

3/4 

4D 

19.5  ±  1.3 

7B/8B 

3644 

5/8 

4D 

22.4  ±  2.1 

9B/8B 

3636 

5/8 

4D 

18.0  ±  3.7 

10B/8B 

3648 

5/8 

4D 

25.5  ±  2.0 

11B/14B 

3642 

5/8 

3D 

44.5  ±  2.6 

12B/14B 

3632 

5/8 

3D 

34.7  ±  2.5 

13B/14B 

3646 

5/8 

3D 

45.9  ±  2.5 

15B/16B 

3633 

3/4 

3D 

27.8  ±  1.8 

18B/17B 

3642 

9/16 

4D 

108.3  ±2.3 

19B/21B 

4063 

15/16 

2D 

53.1  ±  3.1 

19B/22B 

4063 

15/16 

2D 

51.7  ±  1.6 

20B/21B 

4058 

15/16 

2D 

33.3  ±  2.0 

20B/22B 

4058 

15/16 

2D 

30.2  ±  1.9 

♦Deflection 

Values  are 

the  Mean 

and  the 

Dispersion  of  Results 

from  Four  Data  Reduction  Runs  Using  Different  Blade  Array  Sets 
for  the  Input  Time  History  Data. 
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TABLE  43 


BLADE  TIP  DEFLECTIONS  CALCULATED  FROM  ISRL  NSMS  DATA 
FOR  THE  THIRD  TEST  SERIES 


Mode 

Resonance 
Speed (rpm) 

Gap 

(inch) 

Peak-Peak 

Amplitude 

(mils) 

4D 

3631 

9/16 

30  +  2 

3D 

3600 

7/16 

50  ±  3.0 

2D 

4019 

1/2 

38.3  +  1.5 

2D 

4044 

1/2 

41.5  ±  0.5 

The  wholefield  deflection  patterns  for  the  test  disk  were 
recorded  by  double  pulse  inter ferograms  and  the  inter ferogram 
fringes  were  interpreted  to  determine  maximum  displacements  at 
blade  tips  by  the  calculation  procedure  described  in  SECTION  IV. 
Peak-to-peak  amplitudes  were  evaluated  based  on  pulse  width 
separation  time  compared  to  the  disk  response  frequency.  AEDC  and 
UDRI  independently  calculated  blade  tip  displacements  from  the 
recorded  interferograms,  and  these  results  are  included  in  Table 
47.  A  comparison  of  the  optical  and  interferometric  results  in 
this  table  for  selected  test  runs  presents  a  range  for  blade 
deflections  defined  by  counting  the  number  of  fringes  in  different 
sectors  of  the  2,  3,  and  4D  integral  order  mode  shapes.  Reasonable 
agreement  was  found  between  the  AEDC  and  UDRI  interferometric 
fringe  analysis  and  between  the  interferometric  and  NSMS  deflection 
data. 

8.3.2  The  Z-Grid  System 

The  data  acquired  by  the  integrated  Z-Grid  and  NSMS 
Systems  described  in  Section  VII. 5  was  analyzed  and  the  results 
from  tne  2,  3,  and  4D  modal  deflection  tests  for  selected 
excitation  magnet  spacings  are  summarized  in  Table  48.  Similarity 
in  the  NSMS  and  Z-Grid  time  series  data  could  be  seen  from  Figures 
108  and  109.  Comparing  the  time  series  and  the  frequency 
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TABLE  44 


UTRC  NSMS  FOUR-SENSOR  ANALYSIS  RESULTS  -  FIRST  TEST  SERIES 


Run  Mux 
No.  No. 


TEST  CONFIGURATION 


RESULTS 


Rotor 

Speed 

(rpm) 

Excita¬ 

tion 

(E/O) 

Magnet 

Gap 

(in.) 

Sensor 

Spacing 

(cleg) 

Mode 

Order 

Frequency 

(Hz) 

Amplitude 
(mils  p-p) 

Phase  Lag 
(deg) 

3720 

3 

5/16 

36.36.36 

3 

186.0 

110 

31 

3 

186.0 

108 

32 

3 

186.0 

110 

30 

3720 

3 

5/16 

36.36.36 

3 

186.0 

110 

31 

3 

186.0 

112 

3 

186.0 

112 

3970 

2 

5/16 

36.36.36 

2 

132.4 

198 

28 

2 

132.3 

196 

2 

132.4 

200 

3970 

2 

5/16 

36.36.36 

2 

132.3 

196 

29 

2 

132.4 

196 

2 

132.3 

198 

3780 

2 

7/16 

54.54.54 

2 

126.0 

140 

41 

2 

126.0 

140 

2 

126.0 

160 

3780 

2 

7/16 

54.54.54 

2 

126.0 

130 

41 

2 

126.0 

138 

2 

126.0 

140 

3630 

4 

7/16 

36.18.18 

n 

242.0 

50 

27 

242.0 

44 

242.0 

50 

3630 

4 

7/16 

38.18.18 

242.0 

50 

27 

242.0 

50 

■■ 

242.0 

50 
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RUN  13Z/14Z  ZGRID  3D  7/161a.  gap  f  367e/3eM  rpm 


D«i«  SaaplM  •  ll/Hs 


Figure  109.  Z-Grid  Time  Series  Data  for  the  same  Sample  Rate. 
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TABLE  46 


NSMS  BLADE  TIP  DEFLECTION  COMPARISON  -  FIRST  TEST  SERIES 
(PEAK-TO-PEAK  DEFLECTIONS  IN  MILS) 


Mode 

Gap 

(in) 

Speed 

(rpm) 

UDRI 

UTRC 

2D 

5/16 

3790 

No  Data 

198-200 

2D 

7/16 

3780 

176±5 

130-160 

3D 

5/16 

3732 

160±4 

No  Data 

3D 

7/16 

3720 

110±2 

110-112 

4D 

7/16 

3630 

62±2 

44-50 

components  of  a  specific  mode,  the  Z-Grid  measurement  method  often 
yielded  lower  blade  tip  deflections  than  those  from  the  NSMS 
method.  This  would  be  possible  if  the  source  magnet  was  located 
near  the  nodal  diameter.  Also  the  nonsymmetry  of  time  series 
indicated  noise  in  the  grid  electronic  circuits.  Another 
explanation  for  a  lower  disk  response  measurement  by  the  Z-Grid 
system,  in  general  was  the  centrifugal  stiffening  effect  of  the 
source  magnet  and  its  counterweight.  Additional  masses  would 
increase  the  disk  stiffness  from  the  relationship,  the  centrifugal 

force,  =  ("nagnet  “counterweight)  “ 

was  the  angular  speed,  and  r  was  the  radial  distance  of  the  extra 

mass  from  the  spin  axis. 

UDRI  had  determined  the  axial  and  tangential  displacement 
components  of  the  instrumented  blade,  using  our  digital  Z-Grid 
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system.  The  time  of  arrival  data  from  the  two  grid  circuits  were 
first  corrected  for  speed  changes  and  then  summed  and  differenced 
to  obtain  tangential  and  axial  components,  respectively.  This 
method  of  analysis  was  applied  to  the  data  run  13Z/14Z.  The 
results  from  this  analysis  are  shown  in  Figures  110  through  117. 
Figures  110  and  111  show  segments  of  Z-grid  time  series  data 
samples  of  channels  3  and  4  respectively.  Figure  112  shows  the  sum 
of  the  channel  3  and  4  time  series  signals  and  represents  the 
tangential  motion  of  the  instrumented  blade.  Nearly  zero 
displacement  was  seen,  and  this  correlated  with  the  expected  nearly 
zero  tangential  vibration  of  the  flat-plate  test  disk.  Figure  113 
shows  the  time  series  of  difference  data  from  the  two  channels. 
Peak-to-peak  vibration  amplitude  of  50  mil  was  shown  for  the  axial 
vibration  component.  The  Fourier  components  lent  support  to  the 
analysis  procedure.  Figures  114  and  115  show  the  Fourier 
components  of  the  channel  3  and  4  time  series  data,  whereas  Figures 
116  and  117  show  the  spectral  components  of  tangential  and  axial 
vibration  from  the  sum  and  difference  time  series  data, 
respectively.  A  zero  displacement  in  Figure  116  and  a  half 
amplitude  value  of  25  mil  seen  in  Figure  117  matched  well  with  the 
respective  time  signals  and  also  matched  well  with  the  NSMS 
measurements  taken  at  the  same  test  condition. 

Further  research  will  be  required  to  understand  the  complete 
performance  of  the  Z-grid  system. 

8.4  TRAVELING  WAVE  MODE  SHAPE  CHARACTERISTICS 

For  mode  shape  definition,  holograms  of  the  rotating  disk  in 
resonance  were  made.  Two  additional  holograms,  one  at  a  lower 
speed  and  the  other  at  a  higher  than  critical  speed  were  recorded 
to  show  the  variation  of  the  traveling  wave  position  and  amplitude 
with  speed.  Such  changes  in  a  4D  wave  are  shown  in  Figure  118. 
Similar  changes  in  wave-travel  are  also  noticed  in  Figure  119. 

From  the  of f-resonance  records,  a  stationary  displacement  wave 
of  the  shape  A  cos  n0  could  easily  be  seen.  The  location  of  the 
nodal  diameter  with  respect  to  the  static  magnets  at  speeds 
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Run  14Z/i3Z»  10-Nov-88#  3D  at  3670/3000  rpm 


Run  14Z/132#  ie-Nov-88#  3D  at  3670/3000  rpm 


Figure  113.  Axial  Vibration  Component  Due  to  the  Difference  of 
Ch.  3  and  Ch.  4  Time  Series  Dc.ta. 


Run  14Z/13Z*  10-Nosr’88*  3D  at  3670/3000  rpm 


Figure  115.  Fourier  Components  of  the  Ch.  4  Z-Grid  Data. 


Run  i4Z/13Z#  ie-Nov-88#  3D  at  3678/3086  rpm 


i4Z/i3Z,  le-Nov-SS*  3D  at  3670/3000  rpm 


Figure  117,  Fourier  Components  of  the  Axial  Vibration. 


Magnet  - ;  Magnet  - -  Magnet 


Figure  118.  Traveling  Wave  Positions  Before,  At,  and  After  Resonance  for  4D  Mode. 


Nodal  line 


Figure  119.  Traveling  Wave  Positions  Before,  At,  and  After  Resonance  for  3D  Mode 


TABLE  47 


NSMS  AND  LASER  INTERFEROMETRY  BLADE  TIP  DEFLECTION  COMPARISONS: 

SECOND  TEST  SERIES 

(Peak-to-Peak  Deflections  in  mils) 


UTRC _  _ UDRI 


Gap 

Speed 

Interfe- 

NSMS 

Interfe- 

NSMS 

Mode 

(in) 

(rpm) 

rometry 

rometry 

2D 

15/16 

4063 

23-28 

30-32 

No  Data 

24.8 

51.7 

15/16 

26-30 

31.3 

30.2 

3D 

5/8 

3642 

48-50 

30-50 

5/8 

31-34 

20-36 

50.4 

44.5 

5/8 

48-52 

34-54 

3D 

3/4 

3633 

24-26 

32 

29-31 

18-34 

29.9 

27.8 

29-33 

20-36 

4D 

5/8 

3644 

30-34 

36-54 

22.4 

25-29 

No  Data 

33.2 

13-25 

No  Data 

4D 

9/16 

3642 

29-36 

25-31 

34-38 

34.9 

27.6 

below  and  above  resonance  indicated  the  lead  and  lag  phase  angles 
of  the  force-wave  relative  to  the  response  wave.  When  the  disk  was 
in  resonance,  the  peak  of  the  response  wave  was  at  ff/2  relative  to 
the  excitation  force. 

Near  resonance  amplitude  changes  of  the  traveling  wave  with 
speed  were  illustrated  in  Figure  120.  The  wave  build-up  and 
decline  effects  were  evident  from  the  varying  number  of  fringes  in 
holographic  mode  shapes  when  the  disk  was  below,  at,  and  above  the 
resonant  speed.  Fewer  fringes  were  present  in  holograms  made  at 
the  lower  and  upper  half-power  points  of  the  resonance  curve. 
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p 


4,058  rpm,  Amp  =  14  fringes 


4,060  rpm.  Amp  =  9  fringes 

Figure  120.  Resonant  Wave  Build-up  with  Speed  -  2D  Mode. 


Maximum  number  of  fringes  at  the  resonant  speed  of  4,058  rpm  for 
the  2D  indicated  maximum  dynamic  displacement  for  the  disk.  Also, 
maximum  deflection  occurred  when  the  traveling  wave  node  coincided 
with  the  excitation  force  location.  The  holographic  mode  shape 
records  in  Figures  118  through  120  seemed  to  indicate  the  opposite, 
but  that  was  because  the  double  pulse  holographic  technique 
recorded  the  nodes  of  the  actual  mode  as  antinodes  in  the  image  and 
vice  versa.  The  force-response  relationship  of  resonating  bladed 
disks  at  critical  speeds  is  further  discussed  in  (54) . 

8.5  BLADE-TO-BLADE  MISTUNING  STUDIES 


A  perfect  axisymmetric  disk  under  symmetric  load  conditions 
will  produce  symmetrical  modal  response,  but  our  test  disk  with 
several  direct  mounted  motion  sensing  devices  and  their  leadwires 
altered  the  disk  stiffness  locally  and  produced  different  blade 
responses. 

The  UDRI  NSMS  software  was  used  to  analyze  the  test  run  of 
August  3,  1988  and  extract  the  blade-to-blade  deflections.  The 
details  of  this  test  and  results  were  given  in  Table  49. 

The  routine  NSMS  data  analysis  procedure  for  the  above  test 
conditions  yielded  average  blade  deflections  35,  45,  and  45  mils 
when  the  disk  transited  through  resonance.  The  presence  of 
harmonics  other  than  the  primary  component  was  evident  for  the  off- 
resonance  conditions  of  the  disk  due  to  the  participation  from 
other  modes.  Blade  deflections  for  the  upper  half  power  point  were 
higher  than  those  of  resonance  condition  for  the  3D  mode  but  the 
phase  change  expected  was  not  observed. 
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TABLE  48 

BLADE  TIP  AMPLITUDES  (MILS)  COMPARISON 
BY  NSMS  AND  Z-GRID  METHODS 


Mode 

Magnet 
Spacing 
^in. ) 

Critical 

Speed 

rRPM) 

NSMS 

fD-D.mils) 

Z-Grid 

fo-D.mils) 

2D 

1/2 

4044 

41.5 

34.0  ±  3 

15/16 

4058 

30.2  ±  1.9 

— 

3D 

7/16 

3670 

60  ±  1.0 

58  +  1 

3/4 

3633 

27.8  ±  1.8 

— 

4D 

9/16 

3631 

30.1  ±  2 

— 

5/8 

3644 

22.4 

11  ±  1 

To  arrive  at  the  average  deflection  of  an  array  of  blades,  the 
NSMS  data  was  rearranged  and  data  values  from  several  blades  were 
used  to  construct  a  single  time  series,  which  was  used  in  turn  as 
an  input  for  the  FFT  software  analysis.  The  data  for  the  four  sets 
of  symmetrically  arranged  blades  was  analyzed  to  produce  the 
average  deflection  spectrums.  Three  such  blade  ordering  schemes 
were  already  presented  in  Tables  27,  28,  and  29  in  SECTION  VI. 
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TABLE  49 


TEST  DISK  RESONANCE  TEST  CONDITIONS 


Run 

# 

Mode 

Speed 

(rpm) 

Gap 

(in.) 

Data  Acquired  at 

Average 

Deflections 

(p-p,  mils) 

12B/14B 

3D 

3632 

5/8 

Lower  half-power  point 

Array 

1  32 

Array 

2  37 

Array 

3  37 

Array 

4  31 

11B/14B 

3D 

3642 

5/8 

Resonance 

Array 

1  44 

Array 

2  47 

Array 

3  47 

Array 

4  46 

13B/14B 

3D 

3646 

5/8 

Upper  half-power  point 

Array 

1  43 

Array 

2  48 

Array 

3  49 

Array 

4  43 

In  the  present  study,  data  readings  for  each  blade  when  the 
disk  was  in  resonance  weic^  extracted.  Each  data  value  was 
corrected  for  speed  variation  and  the  statistical  mean  over  500 
rotations  calculated.  The  nonvibratory  data  at  3,000  rpm  was 
similarly  treated.  The  resonance  data  were  then  referenced  with 
nonvibration  data  to  get  the  number  of  counts  for  each  of  the 
blades.  These  counts  were  a  direct  measure  of  blade  tip 
deflections.  Histograms,  the  plot  of  number  of  counts  versus  blade 
numbers  were  then  plotted.  Figure  121,  shows  the  blade-to-blade 
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130 


Figure  \2l •  NSMS  Determined  Blade 


deflections  for  the  optical  probe  2,  when  the  disk  was  stabilized 
at  the  lower  half-power  point  speed,  3,632  rpm.  Similarly  the 
Figure  122  shows  resonant  blade  deflections  at  3,642  rpm. 

Fringe  number  differences  in  different  sectors  of  a  specific 
modal  pattern  were  attributed  to  imperfection  in  the  disk,  that  is, 
mistuning.  Such  asymmetric  displacement  patterns  for  the  2,  3,  and 
4D  modes  were  recorded  in  our  experiments.  A  typical  modal 
asymmetric  fringe  pattern  for  3D  before  and  at  resonance  are  shown 
in  Figure  123.  This  behavior  was  evident  in  3D  and  4D  resonant 
mode  shapes  of  the  disk  shown  in  Figure  124.  The  offset  effect 
seen  in  recorded  mode  shapes  was  due  to  the  back  and  forth 
switching  between  two  closely  spaced  resonant  frequencies  for 
different  sectors  of  the  disk.  This  was  likely  due  to  the  variable 
stiffness  and  mass  distributions  of  the  instrumentation  system 
installed  to  this  very  light  test  disk.  Fringe  muddling  due  to 
bias  fringes  was  minimum  at  resonance  speeds  and  a  sample  bias 
fringe  hologram  is  shown  in  Figure  125. 

8 . 6  NONLINEAR  DYNAMICS 

Tests  were  also  conducted  with  various  magnet  spacings  ranging 
from  1/2  to  1  inch  to  study  the  effect  of  excitation  force  on  the 
dynamic  response  of  the  bladed  disk.  Amplitude  changes  due  to 
speed  increments  in  steps  of  l  rpm  on  either  side  of  resonance  were 
recorded  for  the  2D,  3D,  and  4D  modes.  A  phase  change  of  nearly  90 
degrees  as  the  disk  passed  through  resonance  was  evident  from 
Figures  126  and  127  for  the  2D  and  3D  modes.  A  smooth  phase  change 
did  not  occur  when  the  magnets  were  positioned  closer  than  9/16 
inch.  This  was  attributed  to  the  nonlinearity  of  the  disk.  The 
nonlinear,  hard  spring  response  characteristic  of  the  rotating 
bladed  disk  is  presented  in  Figure  128. 
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AT  3,642  RPM 
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®19'®24'®31 
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-  Blade  Number  Deflection  in  mils 

Figure  123 .  Modal  Asymmetry  due  to  Disk  Imperfection 
at  Two  Speeds  for  3D. 


39  mils 


~  33.5  mil 


50.5  mils 


-  46  mils 


271 


.tude 


3D  ReSPONSC-CflP-5^0  IN. 


2 


128.  Hard  Spring  Characteristics  of  Rotating 
Flexible  Disk. 


SECTION  IX 

CONCLUSIONS  AND  RECOMMENDATIONS 


UDRI  has  concentrated  its  efforts  to  pursue  the  objectives  of 
the  proposed  program  for  developing  advanced  analytical, 
ex^jerimental  and  instrumentation  methods  to  characterize  the 
dynamics  of  rotating  turbine  engine  components.  This  research  and 
development  work  performed  during  the  period  1  September  1985 
through  30  June  1989,  has  resulted  in  many  accomplishments.  Two 
nonintrusive  instrumentation  systems,  one  based  on  plane-of-light 
optical  sensors  and  the  other  based  on  digital  demodulation  of  twin 
gridwires  were  developed  and  demonstrated  for  rotating  blade 
vibration  measurements.  Finite  element  based  deterministic  and 
probabilistic  analysis  methods  were  developed  and  used  to  analyze 
turbine  engine  components  for  vibrations.  This  section  provides 
the  conclusions  from  the  program  and  offers  some  recommendations 
for  future  research. 

9.1  CONCLUSIONS 


The  following  conclusl — ^  are  based  upon  the  analytical  and 
experimental  investigation  of  turbine  engine  component  dynamics. 

9.1.1  NSMS  Instrumentation 

Four  nonintrusive  plane-of-light  fiber  optic  sensors 
for  blade  tip  deflection  measurements  of  rotating  bladed  disks  and 
the  associated  data  acquisition,  analysis,  and  display  systems  were 
developed,  demonstrated  and  delivered  to  ISRL. 

■  This  system  is  well  suited  for  laboratory  use  and  has  an 
inexpensive  white  light  source  for  illumination. 

■  Its  fast  turn  around  time  allows  the  user  to  calculate  test 
results,  plan  the  next  appropriate  test  condition  and 
conduct  the  next  test  quickly. 

■  Stationary  wave  measurements  on  a  12-lnch  diameter  disk 
using  a  single  sensor  provided  accurate  results. 
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■  In  a  joint  effort  with  AEDC  (UTRC)  on  critical  speed 
response  measurements,  the  UDRI  NSMS  results  compared  well 
with  those  from  the  UTRC  system.  Our  system  provided 
results  faster  than  the  UTRC  system. 

■  Fast  rise  time  of  1  PPR  trigger  signal  was  achieved  by  the 
UDRI  devised  folded  beam  path. 

■  Noise  in  the  digital  data  due  to  speed  changes  during  each 
revolution  was  eliminated  by  a  special  UDRI  developed 
algorithm  which  took  these  speed  changes  into  account. 

This  feature  was  not  available  with  the  AEDC  System. 

■  The  Data  processing  technique  was  suitably  modified  to 
capture  the  blade-to-blade  vibratory  response.  This 
capability  has  potential  for  application  in  real  turbine 
systems . 

■  Besides  histograms,  the  UDRI  system  displays  real-time  data 
for  five  blades.  This  feature  was  unique  to  the  UDRI 
NSMS. 

9.1.2  2-Grid  Instrumentation 

A  novel  2 -channel  grid  wire  system  was  designed  and 
tested  for  monitoring  an  instrumented  blade  tip  vibrations  of  a 
rotating  disk.  This  system  was  delivered  to  ISRL. 

■  The  design  of  this  system  is  unique.  Demodulation  was 
achieved  by  digital  circuits  unlike  the  analog  methods 
used  by  Rolls  Royce  and  G.E. 

■  The  Z-grid  signal  conditioning  units  were  integrated  into 
the  NSMS  data  acquisition  system. 

■  The  common  data  acquisition  system  for  both  Z-grid  and  NSMS 
was  novel.  However,  it  affected  the  system  speed  and 
stability. 

■  The  Z-grid  system  was  demonstrated  in  conjunction  with  the 
optical  NSMS  for  rotating  disk  critical  speed  experiments. 
This  system  consistently  gave  lower  blade  tip  deflections 
than  the  NSMS  and  interferometric  methods. 
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■  We  recognized  the  need  to  vary  the  trigger  signal  delay 
time  to  compensate  the  speed  and  signal  notch  effects. 

This  will  hopefully  increase  the  accuracy  of  measurements. 

■  A  UDRI  developed  algorithm  resolved  successfully  the  grid 
data  into  tangential  and  axial  deflection  components. 

■  Comparing  the  NSMS  and  Z-Grid  systems,  the  NSMS  System 
provided  data  samples  on  all  blades  of  the  disk  whereas  the 
Z-Grid  system  provided  data  samples  at  a  much  higher  rate 
on  a  single  blade.  The  data  sampling  rate  ratios  for  our 
optical  and  grid  systems  are  1:6,  1:4,  and  1:3  for  the  2D, 
3D,  and  4D  modes  respectively. 

■  Higher  sampling  rates  are  possible  by  increasing  the  number 
of  probes  in  NSMS  and  decreasing  the  grid  spacing  in  Z-Grid 
systems . 

9.1.3  General  Instrumentation  Development 

Some  important  developments  in  the  mechanical,  optical 
and  electronic  instrumentation  which  have  furthered  the 
experimental  research  at  ISRL  are: 

■  The  UDRI  designed  10-channel  slip  ring  unit  allowed  test 
disks  to  rotate  at  higher  spin  rates  than  before.  This 
direct  mount  design  also  eliminated  downtime. 

■  For  AEDC  critical  speed  tests,  UDRI  provided  a  unique 
design  for  the  multimagnet  configurations  to  excite  pure 
backward  traveling  waves  in  a  rotating  disk. 

■  The  sine  extrapolation  equation  for  interpreting 
interferogram  fringes  for  the  full-field  modal  deflection 
definition  has  been  explicitly  derived.  The  user  simply 
Inputs  the  number  of  fringes  and  exposure  time  to  determine 
the  disk  deflection  topology. 

■  The  Ruby  laser  has  been  repaired  several  times,  and  checked 
to  establish  double  and  triple  pulse  operation.  Pulse 
width  separations  have  been  accurately  determined  by  a 
digital  storage  oscilloscope. 
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■  Pockels  cells  extinction  ratios  could  not  be  determined  for 
lack  of  suitable  equipment  which  curtailed  triple  pulse 
interferometry  implementation. 

9.1.4  Analytical  Research 

Analytical  techniques  have  been  developed  and  used  to 
understand  the  dynamic  behavior  of  simple  to  complex  turbine 
components.  Some  important  developments  were: 

■  Time  and  frequency  domain  forced  response  of  rotating 
bladed  disk  assemblies  with  mass  and  stiffness  asymmetry 
were  explored.  Amplitude  and  location  of  anti-nodal  lines 
were  significantly  affected  by  mass  distribution  around  the 
disk. 

■  Stiffness  matrix  singularities  resulting  from  highly  curved 
and  twisted  geometries  of  advanced  blade  profiles  were 
observed  in  the  analysis.  A  procedure  to  suppress  these 
singularities  was  developed  and  applied  to  low  aspect  ratio 
fan  (LARF)  blades  and  the  results  were  compared  with  those 
of  PW  and  G.E.  A  good  correlation  for  all  modes  except  one 
was  observed. 

■  An  integral  bladed  (Wennerstrom)  disk  was  analytically 
modeled  by  PATRAN  and  analyzed  by  the  NASTRAN  cyclic 
symmetry  feature  for  vibrations.  Approximate  equations  to 
couple  the  interfacial  displacements  of  the  disk  and  blade 
was  developed.  Experimental  and  analytical  results  for 
non-rotating  conditions,  agreed  well. 

■  Probabilistic  finite  element  code  (PROTEC)  was  developed  to 
account  for  uncertainties  in  the  element  mass,  geometry  and 
material  properties.  This  code  has  been  successfully 
applied  to  analyze  stress,  frequency  and  deflection 
variances  of  the  AEDC  disk.  The  need  to  upgrade  this  code 
to  account  for  centrifugal  effects  is  recognized. 
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9 . 2  RECOMMENDATIONS 

The  following  recommendations  are  offered  to  the  Air  Force. 
These  are  directly  related  to  the  current  program  and  if 
implemented,  will  enhance  the  engine  structural  design  technology. 

■  In  the  present  investigation,  efforts  were  concentrated  on 
the  determination  of  one  mode  (cosine  or  sine  wave) 
response  of  a  bladed  disk.  It  is  recommended  that  the  twin 
mode  critical  speed  response  investigation  of  a  rotating 
disk  be  investigated  so  that  the  disk  modal  interpretation 
could  be  applicable  to  practical  turbine  disks. 

■  The  experiments  conducted  on  the  38-bladed  disk  required 
testing  within  a  specific  frequency  range  which  constrained 
the  test  specimen  to  a  very  flexible  disk  where  nonlinear 
dynamic  behavior  was  easily  triggered.  A  stiff  disk  would 
be  more  stable  over  its  entire  critical  speed  range. 

■  Stability  characteristics  of  twin  modes  of  vibrating  disk 
and  their  interaction  under  a  stationary  system  of  forces 
are  worth  pursuing  to  define  disk  critical  speeds. 

■  Modal  definition  when  flexible  disks  are  near  their 
resonance-speeds  should  be  investigated  in  detail,  since 
future  designs  suggest  the  use  of  lightweight  flexible 
bladed  disks  for  turbine  engines. 

■  Investigation  on  a  realistic  disk  having  both  axial  and 
tangential  blade  motions  is  needed  to  evaluate  the 
measurement  capability  of  NSMS  methods. 

■  Extension  of  the  Z-Grid  system  application  to  monitor 
several  blade  tip  deflections  by  sectored  grids  should  be 
pursued . 

■  The  application  of  the  NSMS  test  methods  to  actual 
compressor  and  turbine  disks  would  be  very  desirable. 

■  Upgrading  the  NSMS  systems  with  several  probes  to 
completely  define  the  deflection  patterns  of  vibrating 
blades  is  required. 
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■  The  present  systems  monitor  blade  tip  deflections,  but  a 
combined  analytical  and  experimental  approach  to  derive 
stress  patterns  for  blade  critical  locations  needs  to  be 
developed. 

■  Multiprobe  and  multigrid  designs  at  leading  and  trailing 
edges  will  allow  complete  definition  of  flexural  and 
torsional  motions  of  turbine  and  compressor  airfoils. 
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